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REFERENCE TO RELATED APPLICATION 

This application claims the benefit of U.S. Provisional Application 
No. 60/080,457, filed on April 2, 1998. 

BACKGROUND OF THE INVENTION 



1. Field of the Invention 

This invention relates generally to power generation and processing systems 
1 0 and in particular to distributed generation power systems. 

2. Description of the Prior Art 

Conventional power generation and distribution systems are configured to 
1 5 maximize the specific hardware used. In the case of a turbine power motor, for 

example, the output or bus voltage in a conventional power distribution system varies 
with the speed of the turbine. In such systems, the turbine speed must be regulated to 
control the output or bus voltage. Consequently, the engine cannot be run too low in 
speed else the bus voltage would not be high enough to generate some of the voltages 
20 that are needed. As a result, the turbine would have to be run at higher speeds and 
lower temperatures, making it less efficient. 

What is needed therefore is a power generation and distribution system where 
the bus voltage is regulated by a bi-directional controller independent of turbine speed. 

25 SUMMARY OF THE INVENTION 

The present invention provides in a first aspect, a power controller which 
provides a distributed generation power networking system in which bi-directional 
power converters are used with a common DC bus for permitting compatibility 
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between various energy components. Each power converter operates essentially as a 
customized bi-directional switching converter configured, under the control of the 
power controller, to provide an interface for a specific energy component to the DC 
bus. The power controller controls the way in which each energy component, at any 
5 moment, will sink or source power, and the manner in which the DC bus is regulated. 
In this way, various energy components can be used to supply, store and/or use power 
in an efficient manner. The various energy components include energy sources, loads, 
storage devices and combinations thereof. 

In another aspect, the present invention provides a turbine system including a 

10 turbine engine, a load, a power controller, an energy reservoir for providing transient 
power to the DC bus and an energy reservoir controller, in communication with the 
power controller for providing control to the energy reservoir. The power controller 
includes an engine power conversion in communication with the turbine engine, an 
utility power conversion in communication with the load and a DC bus. 

15 These and other features and advantages of this invention will become further 

apparent from the detailed description and accompanying figures that follow. In the 
figures and description, numerals indicate the various features of the invention, like 
numerals referring to like features throughout both the drawing figures and the written 
description. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a block diagram of a power controller according to the present 
invention. 

FIG. 2 is a detailed block diagram of a power converter in the power controller 
25 illustrated in FIG. 1. 

FIG. 3 is a simplified block diagram of a turbine system including the power 
architecture of the power controller illustrated in FIG. 1. 

FIG. 4 is a block diagram of the power architecture of a typical implementation 
of the power controller illustrated in FIG. 1 . 
30 FIG. 5 is a schematic diagram of the internal power architecture of the power 

controller illustrated in FIG. 1. 
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FIG. 6 is a functional block diagram of an interface between load/utility grid 
and turbine generator using the power controller according to the present invention. 

FIG. 7 is a functional block diagram of an interface between load/utility grid 
and turbine generator using the power controller for a stand-alone application 
according to the present invention. 

FIG. 8 is a schematic diagram of an interface between a load/utility grid and 
turbine generator using the power controller according to the present invention. 

FIG. 9 is a block diagram of the software architecture for the power controller 
including external interfaces. 

FIG. 10 is a block diagram of an EGT control mode loop for regulating the 
temperature of the turbine. 

FIG. 1 1 is a block diagram of a speed control mode loop for regulating the 
rotating speed of the turbine. 

FIG. 12 is a block diagram of a power control mode loop for regulating the 
power producing potential of the turbine. 

FIG. 13 is a state diagram showing various operating states of the power 
controller. 

FIG. 14 is a block diagram of the power controller interfacing with a turbine 
and fuel device. 

FIG. 15 is a block diagram of the power controller in multi-pack configuration. 

FIG. 16 is a block diagram of a utility grid analysis system for the power 
controller according to the present invention. 

FIG. 17 is a graph of voltage against time for the utility grid analysis system 
illustrated in FIG. 16. 

FIG. 1 8 is a diagram of the power controller shown in FIG. 1 6, including brake 
resistor. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS) 
Referring to FIG.l, power controller 10 provides a distributed generation power 
networking system in which bi-directional (i.e. reconfigurable) power converters are 
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used with a common DC bus for permitting compatibility between one or more energy 
components. Each power converter operates essentially as a customized bi-directional 
switching converter configured, under the control of power controller 10, to provide an 
interface for a specific energy component to DC bus 24. Power controller 10 controls 
the way in which each energy component, at any moment, will sink or source power, 
and the manner in which DC bus 24 is regulated. In this way, various energy 
components can be used to supply, store and/or use power in an efficient manner. 

One skilled in the art will recognize that the particular configurations shown 
herein are for illustrative purposes only. In particular, the present invention is not 
limited to the use of three bi-directional converters as shown in FIG. 1. Rather, the 
number of power converters is dependent on various factors, including but not limited 
to, the number of energy components and the particular power distribution 
configuration desired. For example, as illustrated in FIGS. 5 and 6, power controller 
10 can provide a distributed generation power system with as few as two power 
converters. 

The energy components, as shown in FIG. 1, include energy source 12, 
utility/load 18 and storage device 20. The present invention is not limited to the 
distribution of power between energy source 12, energy storage device 20 and 
utility/load 18, but rather may be adapted to provide power distribution in an efficient 
manner for any combination of energy components. 

Energy source 12 may be a gas turbine, photovoltaics, wind turbine or any other 
conventional or newly developed source. Energy storage device 20 may be a flywheel, 
battery, ultracap or any other conventional or newly energy storage device. Load 1 8 
may be an utility grid, dc load, drive motor or any other conventional or newly 
developed utility/load 18. 

Referring now to FIG. 2, a detailed block diagram of power converter 14 in 
power controller 10, shown in FIG. 1, is illustrated. Energy source 12 is connected to 
DC bus 24 via power converter 14. Energy source 12 may be, for example, a gas 
turbine driving an AC generator to produce AC which is applied to power converter 14. 
DC bus 24 connects power converter 14 to utility/load 18 and additional energy 
components 36. Power converter 14 includes input filter 26, power switching system 
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28, output filter 34, signal processor 30 and main CPU 32. In operation, energy source 
12 applies AC to input filter 26 in power converter 14. The filtered AC is then applied 
to power switching system 28 which may conveniently be a series of insulated gate 
bipolar transistor (IGBT) switches operating under the control of signal processor (SP) 
5 30 which is controlled by main CPU 32. One skilled in the art will recognize that other 
conventional or newly developed switches may be utilized as well. The output of the 
power switching system 28 is applied to output filter 34 which then applies the filtered 
DC to DC bus 24. 

In accordance with the present invention, each power converter 14, 16 and 22 
10 operates essentially as a customized, bi-directional switching converter under the 
control of main CPU 32, which uses SP 30 to perform its operations. Main CPU 32 
provides both local control and sufficient intelligence to form a distributed processing 
system. Each power converter 14, 16 and 22 is tailored to provide an interface for a 
yQ specific energy component to DC bus 24. Main CPU 32 controls the way in which 

j7j 1 5 each energy component 12,18 and 20 sinks or sources power, and DC bus 24 is 
W regulated at any time. In particular, main CPU 32 reconfigures the power converters 

m 14, 16 and 22 into different configurations for different modes of operation. In this 

^ way, various energy components 12, 1 8 and 20 can be used to supply, store and/or use 

O power in an efficient manner. In the case of a turbine power generator, for example, a 

ry 20 conventional system regulates turbine speed to control the output or bus voltage. In the 
jjf power controller, the bi-directional controller independently of turbine speed regulates 

N« the bus voltage. 

Operating Modes 

FIG. 1 shows the system topography in which DC bus 24, regulated at 800 v 
25 DC for example, is at the center of a star pattern network. In general, energy source 12 
provides power to DC bus 24 via power converter 14 during normal power generation 
mode. Similarly, during the power generation mode, power converter 16 converts the 
power on DC bus 24 to the form required by utility /load 18, which may be any type of 
load including a utility web. During other modes of operation, such as utility start up, 
30 power converters 14 and 16 are controlled by the main processor to operate in different 
manners. 
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For example, energy is needed to start the turbine. This energy may come from 
load/utility grid 18 (utility start) or from energy storage 20 (battery start), such as a 
battery, flywheel or ultra-cap. During a utility start up, power converter 16 is required 
to apply power from load 1 8 to DC bus 24 for conversion by power converter 14 into 
the power required by energy source 12 to startup. During utility start, energy source 
or turbine 12 is controlled in a local feedback loop to maintain the turbine revolutions 
per minute (RPM). Energy storage or battery 20 is disconnected from DC bus 24 while 
load/utility grid 10 regulates V DC on DC bus 24. 

Similarly, in the battery start mode, the power applied to DC bus 24 from which 
energy source 12 is started may be provided by energy storage 20 which may be a 
flywheel, battery or similar device. Energy storage 20 has its own power conversion 
circuit in power converter 22, which limits the surge current into DC bus 24 capacitors, 
and allows enough power to flow to DC Bus 24 to start energy source 12. In particular, 
power converter 16 isolates DC bus 24 so that power converter 14 can provide the 
required starting power from DC bus 24 to energy source 12. 
Electronics Architecture 

Referring to FIG. 3, a simplified block diagram of a turbine system 50 using the 
power controller electronics architecture of the present invention is illustrated. The 
turbine system 50 includes a fuel metering system 42, turbine engine 58, power 
controller 52, energy reservoir conversion 62, energy /reservoir 64 and load/utility grid 
60. The fuel metering system 42 is matched to the available fuel and pressure. The 
power controller 52 converts the electricity from turbine engine 58 into regulated DC 
then converts it to utility grade AC electricity. By separating the engine control from 
the converter that creates the utility grade power, greater control of both processes is 
realized. All of the interconnections are comprised of a communications bus and a 
power connection. 

The power controller 52 includes an engine power conversion 54 and utility 
power conversion 56 which provides for the two power conversions that take place 
between the turbine 58 and the load/utility grid 60. One skilled in the art will recognize 
that the power controller 52 can provide a distributed generation power system with as 
few as two power converters 54 and 56. The bi-directional (i.e. reconfigurable) power 
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converters 54 and 56 are used with a common regulated DC bus 66 for permitting 
compatibility between the turbine 58 and load/utility grid 60. Each power converter 54 
and 56 operates essentially as a customized bi-directional switching converter 
configured, under the control of the power controller 10, to provide an interface for a 
5 specific energy component 58 or 60 to the DC bus 66. The power controller 10 

controls the way in which each energy component, at any moment, will sink or source 
power, and the manner in which the DC bus 66 is regulated. Both of these power 
conversions 54 and 56 are capable of operating in a forward or reverse direction. This 
allows starting the turbine 58 from either the energy reservoir 64 or the load/utility grid 

10 60. The regulated DC bus 66 allows a standardized interface to energy reservoirs such 
as batteries, flywheels, and ultra-caps. The architecture of the present invention 
permits the use of virtually any technology that can convert its energy to/from 
electricity. Since the energy may flow in either direction to or from the energy 
reservoir 64, transients may be handled by supplying energy or absorbing energy. Not 

1 5 all systems will need the energy reservoir 64. The energy reservoir 64 and its energy 
reservoir conversion 62 are not contained inside the power controller 52. 

Referring to FIG. 4, the power architecture 68 of a typical implementation of 
the power controller 70 is shown. The power controller 70 includes a generator 
converter 72 and output converter 74 which provides for the two power conversions 

20 that take place between the turbine 76 and the load/utility grid 78. In particular, the 
generator converter 72 provides for AC to DC power conversion and the output 
converter 74 provides for DC to AC power conversion. Both of these power converters 
72 and 74 are capable of operating in a forward or reverse direction. This allows 
starting the turbine 76 from either the energy storage device 86 or the load/utility grid 

25 78. Since the energy may flow in either direction to or from the energy storage device 
86, transients may be handled by supplying energy or absorbing energy. The energy 
storage device 86 and its DC converter 84 are not contained inside the power controller 
70. The DC converter 84 provides for DC to DC power conversion. 

Referring to FIG. 5, a schematic 90 of a typical internal power architecture, 

30 such as that shown in FIG. 4, is shown. The turbine has an integral PMG that can be 
used as either a generator (for starting) or a generator (normal mode of operation). 
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Because all of the controls can be performed in the digital domain and all switching 
(except for one output contactor) is done with solid state switches, it is easy to shift the 
direction of the power flow as needed. This permits very tight control of the turbine 
during starting and stopping. In a typical configuration, the power output is a 480 
5 VAC, 3-phase output. One skilled in the art will recognize that the present invention 
may be adapted to provide for other power output requirements such as a 3 -phase, 400 
VAC, and single-phase, 480 VAC. 

Power controller 92 includes generator converter 94 and output converter 96. 
Generator converter 94 includes IGBT switches 94, such as a seven-pack IGBT module 

1 0 94, driven by control logic 98, providing a variable voltage, variable frequency 3-phase 
drive to the PMG 100. Inductors 102 are utilized to minimize any current surges 
associated with the high frequency switching components which may affect the PMG 
100 to increase operating efficiency. 

IGBT module 94 is part of the electronics that controls the engine of the 

1 5 turbine. IGBT module 94 incorporates gate driver and fault sensing circuitry as well as 
a seventh IGBT used to dump power into a resistor. The gate drive inputs and fault 
outputs require external isolation. Four external, isolated power supplies are required 
to power the internal gate drivers. IGBT module 94 is typically used in a turbine 
system that generates 480 VAC at its output terminals delivering up to 30 kWatts to a 

20 freestanding or utility-connected load. During startup and cool down (and occasionally 
during normal operation), the direction of power flow through the seven-pack reverses. 
When the turbine is being started, power is supplied to the DC bus 1 12 from either a 
battery (not shown) or from the utility grid 108. The DC is converted to a variable 
frequency AC voltage to generator the turbine. 

25 For utility grid connect operation, control logic 1 10 sequentially drives the solid 

state IGBT switches, typically configured in a six-pack IGBT module 96, associated 
with load converter 96 to boost the utility voltage to provide start power to the 
generator converter 94. The IGBT switches in load converter 96 are preferably 
operated at a high (15 kHz) frequency, and modulated in a pulse width modulation 

30 manner to provide four quadrant converter operation. Inductors 104 and AC filter 
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capacitors 106 are utilized to minimize any current surges associated with the high 
frequency switching components which may affect load 108. 

Six-pack IGBT module 96 is part of the electronics that controls the converter 
of the turbine. IGBT module 96 incorporates gate driver and fault sensing circuitry. 
The gate drive inputs and fault outputs require external isolation. Four external, 
isolated power supplies are required to power the internal gate drivers. IGBT module 
96 is typically used in a turbine system that generates 480 VAC at its output terminals 
delivering up to approximately 30 k Watts to a free-standing or utility-connected load. 
After the turbine is running, six-pack IGBT module 96 is used to convert the regulated 
DC bus voltage to the approximately 50 or 60 hertz utility grade power. When there is 
no battery (or other energy reservoir), the energy to run the engine during startup and 
cool down must come from utility grid 108. Under this condition, the direction of 
power flow through the six-pack IGBT module 96 reverses. DC bus 1 12 receives its 
energy from utility grid 108, using six-pack IGBT module 96 as a boost converter (the 
power diodes act as a rectifier). The DC is converted to a variable frequency AC 
voltage to generator the turbine. To accelerate the engine as rapidly as possible at first, 
current flows at the maximum rate through seven-pack IGBT module 94 and also six- 
pack IGBT module 96. 

Dual IGBT module 1 14, driven by control logic 1 16, is used to provide an 
optional neutral to supply 3 phase, 4 wire loads. 
Startup 

Energy is needed to start the turbine. Referring to FIGS. 3 and 4, this energy 
may come from utility grid 60 or from energy reservoir 64, such as a battery, flywheel 
or ultra-cap. When utility grid 60 supplies the energy, utility grid 60 is connected to 
power controller 52 through two circuits. First is an output contactor that handles the 
full power (30 kWatts). Second is a "soft-start" or "pre-charge" circuit that supplies 
limited power (it is current limited to prevent very large surge currents) from utility 
grid 66 to DC bus 62 through a simple rectifier. The amount of power supplied 
through the soft-start circuit is enough to start the housekeeping power supply, power 
the control board, and run the power supplies for the IGBTs, and close the output 
contactor. When the contactor closes, the IGBTs are configured to create DC from the 
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AC waveform. Enough power is created to run the fuel metering circuit 42, start the 
engine, and close the various solenoids (including the dump valve on the engine). 

When energy reservoir 64 supplies the energy, energy reservoir 64 has its own 
power conversion circuit 62 that limits the surge circuit into DC bus capacitors. 
Energy reservoir 64 allows enough power to flow to DC bus 62 to run fuel-metering 
circuit 42, start the engine, and close the various solenoids (including the dump valve 
on the engine). After the engine becomes self-sustaining, the energy reservoir starts to 
replace the energy used to start the engine, by drawing power from DC bus 62. In 
addition to the sequences described above, power controller senses the presence of 
other controllers during the initial power up phase. If another controller is detected, the 
controller must be part of a multi-pack, and proceeds to automatically configure itself 
for operation as part of a multi-pack. 
System Level Operation 

Referring to FIG. 6, a functional block diagram 130 of an interface between 
utility grid 132 and turbine generator 148 using power controller 136 of the present 
invention is shown. In this example, power controller 136 includes two bi-directional 
converters 138 and 140. Permanent magnet generator converter 140 starts turbine 148 
(using the generator as a generator) from utility or battery power. Load converter 138 
then produces AC power using an output from generator converter 140 to draw power 
from high-speed turbine generator 148. Power controller 136 also regulates fuel to 
turbine 148 and provides communications between units (in paralleled systems) and to 
external entities. 

During a utility startup sequence, utility 132 supplies starting power to turbine 
148 by "actively" rectifying the line via load converter 138, and then converting the 
DC to variable voltage, variable frequency 3 -phase power in generator converter 136. 
As is illustrated in FIG. 7, for stand-alone applications 150, the start sequence is the 
same as the utility start sequence shown in FIG. 6 with the exception that the start 
power comes from battery 170 under the control of an external battery controller. Load 
1 52 is then fed from the output terminals of load converter 158. 

Referring to FIG. 8, a schematic illustration 180 of an interface between utility 
grid 132 and turbine generator 148 using the power controller is illustrated. Control 
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logic 184 also provides power to fuel cutoff solenoids 198, fuel control valve 200 and 
igniter 202. An external battery controller (not shown), if used, connects directly to 
DC bus 1 90. In accordance with an alternative embodiment of the invention, a fuel 
system (not shown) involving a compressor (not shown) operated from a separate 
variable speed drive can also derive its power directly from DC bus 190. 

In operation, control and start power comes from either the external battery 
controller (for battery start applications) or from the utility, which is connected to a 
rectifier using inrush limiting techniques to slowly charge internal bus capacitor 190. 
For utility grid connect operation, control logic 184 sequentially drives solid state 
IGBT switches 214 associated with load converter 192 to boost the utility voltage to 
provide start power to generator converter 186. Switches 214 are preferably operated 
at a high (1 5 kHz) frequency, and modulated in a pulse width modulation manner to 
provide four quadrant converter operation. In accordance with the present invention, 
load converter 192 either sources power from DC bus 190 to utility grid 222 or from 
utility grid 222 to DC bus 190. A current regulator (not shown) may achieve this 
control. Optionally, two of the switches 214 serve to create an artificial neutral for 
stand-alone applications (for stand-alone applications, start power from an external DC 
supply (not shown) associated with external DC converter 220 is applied directly to DC 
bus 190). 

Solid state (IGBT) switches 214 associated with generator converter 186 are 
also driven from control logic 184, providing a variable voltage, variable frequency 3- 
phase drive to generator 218 to start turbine 208. Control logic 184 receives feedback 
via current sensors I sen s as turbine 206 is ramped up in speed to complete the start 
sequence. When turbine 206 achieves a self sustaining speed of, for example, approx. 
40,000 RPM, generator converter 186 changes its mode of operation to boost the 
generator output voltage and provide a regulated DC bus voltage. 

PMG filter 188 associated with generator converter 186 includes three 
inductors to remove the high frequency switching component from permanent magnet 
generator 208 to increase operating efficiency. Output AC filter 194 associated with 
load converter 192 includes three or optionally four inductors (not shown) and AC 
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filter capacitors (not shown) to remove the high frequency switching component. 
Output contactor 210 disengages load converter 192 in the event of a unit fault. 

During a start sequence, control logic 184 opens fuel cutoff solenoid 198 and 
maintains it open until the system is commanded off. Fuel control 200 may be a 
5 variable flow valve providing a dynamic regulating range, allowing minimum fuel 
during start and maximum fuel at full load. A variety of fuel controllers, including but 
not limited to, liquid and gas fuel controllers, may be utilized. One skilled in the art 
will recognize that the fuel control can be by various configurations, including but not 
limited to a single or dual stage gas compressor accepting fuel pressures as low as 

1 0 approximately l A psig. Igniter 202, a spark type device similar to a spark plug for an 
internal combustion engine, is operated only during the start sequence. 

For stand-alone operation, turbine 206 is started using external DC converter 
220 which boosts voltage from a battery (not shown), and connects directly to the DC 
bus 190. Load converter 192 is then configured as a constant voltage, constant 

1 5 frequency (for example, approximately 50 or 60 Hz) source. One skilled in the art will 
recognize that the output is not limited to a constant voltage, constant frequency 
source, but rather may be a variable voltage, variable frequency source. For rapid 
increases in output demand, external DC converter 220 supplies energy temporarily to 
DC bus 190 and to the output. The energy is restored after a new operating point is 

20 achieved. 

For utility grid connect operation, the utility grid power is used for starting as 
described above. When turbine 206 has reached a desired operating speed, converter 
192 is operated at utility grid frequency, synchronized with utility grid 222, and 
essentially operates as a current source converter, requiring utility grid voltage for 

25 excitation. If utility grid 222 collapses, the loss of utility grid 222 is sensed, the unit 
output goes to zero (0) and disconnects. The unit can receive external control signals 
to control the desired output power, such as to offset the power drawn by a facility, but 
ensure that the load is not backfed from the system. 
Power Controller Software 

30 Referring to FIG. 9, power controller 230 includes main CPU 232, generator SP 

234 and converter SP 236. Main CPU software program sequences events which occur 
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inside power controller 230 and arbitrates communications to externally connected 
devices. Main CPU 232 is preferably a MC68332 microprocessor, available from 
Motorola Semiconductor, Inc. of Phoenix, Arizona. Other suitable commercially 
available microprocessors may be used as well. The software performs the algorithms 
that control engine operation, determine power output and detect system faults. 

Commanded operating modes are used to determine how power is switched 
through the major converts in the controller. The software is responsible for turbine 
engine control and issuing commands to other SP processors enabling them to perform 
the generator converter output converter power switching. The controls also interface 
with externally connected energy storage devices (not shown) that provide black start 
and transient capabilities. 

Generator SP 234 and converter SP 236 are connected to power controller 230 
via serial peripheral interface (SPI) bus 238 to perform generator and converter control 
functions. Generator SP 234 is responsible for any switching which occurs between 
DC bus 258 and the output to generator. Converter SP 236 is responsible for any 
switching which occurs between DC bus 258 and output to load. As illustrated in FIG. 
5, generator SP 234 and converter SP 236 operate IGBT modules. 

Local devices, such as a smart display 242, smart battery 244 and smart fuel 
control 246, are connected to main CPU 232 in power controller 230 via intracontroller 
bus 240, which may be a RS485 communications link. Smart display 242, smart 
battery 244 and smart fuel control 246 performs dedicated controller functions, 
including but not limited to display, energy storage management, and fuel control 
functions. 

Main CPU 232 in power controller 230 is coupled to user port 248 for 
connection to a computer, workstation, modem or other data terminal equipment which 
allows for data acquisition and/or remote control. User port 248 may be implemented 
using a RS232 interface or other compatible interface. 

Main CPU 232 in power controller 230 is also coupled to maintenance port 250 
for connection to a computer, workstation, modem or other data terminal equipment 
which allows for remote development, troubleshooting and field upgrades. 
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Maintenance port 250 may be implemented using a RS232 interface or other 
compatible interface. 

The main CPU processor software communicates data through a TCP/IP stack 
over intercontroller bus 252, typically an Ethernet 10 Base 2 interface, to gather data 
5 and send commands between power controllers (as shown and discussed in detail with 
respect to FIG. 1 5). In accordance with the present invention, the main CPU processor 
software provides seamless operation of multiple paralleled units as a single larger 
generator system. One unit, the master, arbitrates the bus and sends commands to all 
units. 

1 0 Intercontroller bus 254, which may be a RS485 communications link, provides 

high-speed synchronization of power output signals directly between converter SPs, 
such as converter SP 236. Although the main CPU software is not responsible for 
communicating on the intercontroller bus 254, it informs converter SPs, including 
converter SP 236, when main CPU 232 is selected as the master. 

1 5 External option port bus 256, which may be a RS485 communications link, 

allows external devices, including but not limited to power meter equipment and auto 
disconnect switches, to be connected to generator SP 234. 

In operation, main CPU 232 begins execution with a power on self-test when 
power is applied to the control board. External devices are detected providing 

20 information to determine operating modes the system is configured to handle. Power 
controller 230 waits for a start command by making queries to external devices. Once 
received, power controller 230 sequences up to begin producing power. As a 
minimum, main CPU 232 sends commands to external smart devices 242, 244 and 246 
to assist with bringing power controller 230 online. If selected as the master, the 

25 software may also send commands to initiate the sequencing of other power controllers 
(FIG. 15) connected in parallel. A stop command will shutdown the system bringing it 
offline. 

System I/O 

The main CPU 232 software interfaces with several electronic circuits (not 
30 shown) on the control board to operate devices that are universal to all power 

controllers 230. Interface to system I/O begins with initialization of registers within 
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power controller 230 to configure internal modes and select external pin control. Once 
initialized, the software has access to various circuits including discrete inputs/outputs, 
analog inputs/outputs, and communication ports. These external devices may also have 
registers within them that require initialization before the device is operational. 

Each of the following sub-sections provides a brief overview that defines the 
peripheral device the software must interface with. The contents of these sub-sections 
do not define the precise hardware register initialization required. 

Communications 

Referring to FIG. 9, main CPU 232 is responsible for all communication 
systems in power controller 230. Data transmission between a plurality of power 
controllers 230 is accomplished through intercontroller bus 252. Main CPU 232 
initializes the communications hardware attached to power controller 230 for 
intercontroller bus 252. 

Main CPU 232 provides control for external devices, including smart devices 
242, 244 and 246, which share information to operate. Data transmission to external 
devices, including smart display 242, smart battery 244 and smart fuel control 246 
devices, is accomplished through intracontroller communications bus 240. Main CPU 
232 initializes any communications hardware attached to power controller 230 for 
intracontroller communications bus 240 and implements features defined for the bus 
master on intracontroller communications bus 240. 

Communications between devices such as switch gear and power meters used 
for master control functions exchange data across external equipment bus 246. Main 
CPU 232 initializes any communications hardware attached to power controller 230 for 
external equipment port 246 and implements features defined for the bus master on 
external equipment bus 246. 

Communications with a user computer is accomplished through user interface 
port 248. Main CPU 232 initializes any communications hardware attached to power 
controller 230 for user interface port 248. In a typical configuration, at power up, the 
initial baud rate will be selected to 19200 baud, 8 data bits, 1 stop, and no parity. The 
user has the ability to adjust and save the communications rate setting via user interface 
port 248 or optional smart external display 242. The saved communications rate is 
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used the next time power controller 230 is powered on. Main CPU 232 communicates 
with a modem (not shown), such as a Hayes compatible modem, through user interface 
port 248. Once communications are established, main CPU 232 operates as if were 
connected to a local computer and operates as a slave on user interface port 248 (it only 
responds to commands issued). 

Communications to service engineers, maintenance centers, and so forth are 
accomplished through maintenance interface port 250. Main CPU 232 initializes the 
communications to any hardware attached to power controller 230 for maintenance 
interface port 250. In a typical implementation, at power up, the initial baud rate will 
be selected to 1 9200 baud, 8 data bits, 1 stop, and no parity. The user has the ability to 
adjust and save the communications rate setting via user port 248 or optional smart 
external display 242. The saved communications rate is used the next time power 
controller 230 is powered on. Main CPU 232 communicates with a modem, such as a 
Hayes compatible modem, through maintenance interface port 250. Once 
communications are established, main CPU 232 operates as if it were connected to a 
local computer and operates as a slave on maintenance interface port 250 (it only 
responds to commands issued). 

Controls 

Referring to FIG. 9, main CPU 232 orchestrates operation for motor, converter, 
and engine controls for power controller 230. The main CPU 232 does not directly 
perform motor and converter controls. Rather, generator and converter SP processors 
234 and 236 perform the specific control algorithms based on data communicated from 
main CPU 232. Engine controls are performed directly by main CPU 232 (see FIG. 
14). 

Main CPU 232 issues commands via SPI communications bus 238 to generator 
SP 234 to execute the required motor control functions. Generator SP 234 will operate 
the motor (not shown) in either a DC bus mode or a RPM mode as selected by main 
CPU 232. In the DC bus voltage mode, generator SP 234 uses power from the motor 
to maintain the DC bus at the setpoint In the RPM mode, generator SP 234 uses 
power from the motor to maintain the engine speed at the setpoint. Main CPU 232 
provides Setpoint values. 
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Main CPU 232 issues commands via SPI communications bus 238 to converter 
SP 236 to execute required converter control functions. Converter SP 236 will operate 
the converter (not shown) in a DC bus mode, output current mode, or output voltage 
mode as selected by main CPU 232. In the DC bus voltage mode, converter SP 236 
regulates the utility power provided by power controller 230 to maintain the internal 
bus voltage at the setpoint. In the output current mode, converter SP 236 uses power 
from the DC bus to provide commanded current out of the converter. In the output 
voltage mode, converter SP 236 uses power from the DC bus to provide commanded 
voltage out of the converter. Main CPU 232 provides Setpoint values. 

Referring to FIGS. 10-12, control loops 260, 282 and 300 are used to regulate 
engine controls. These loops include exhaust gas temperature (EGT) control (FIG. 10), 
speed control (FIG. 1 1) and power control (FIG. 12). All three of the control loops 
260, 282 and 300 are used individually and collectively by main CPU 232 to provide 
the dynamic control and performance required of power controller 230. These loops 
are joined together for different modes of operation. 

The open-loop light off control algorithm is a programmed command of the 
fuel device used to inject fuel until combustion begins. In a typical configuration, main 
CPU 232 takes a snap shot of the engine EGT and begins commanding the fuel device 
from about 0% to 25% of full command over about 5 seconds. Engine light is declared 
when the engine EGT rises about 28° C (50° F) from the initial snap shot. 

Referring to FIG. 10, EGT control mode loop 260 provides various fuel output 
commands to regulate the temperature of the turbine. Engine speed signal 262 is used 
to determine the maximum EGT setpoint temperature 266 in accordance with 
predetermined setpoint temperature values. EGT setpoint temperature 266 is compared 
by comparator 268 against feedback EGT signal 270 to determine error signal 272, 
which is then applied to a proportional-integral (PI) algorithm 274 for determining the 
fuel command required to regulate EGT at the setpoint. Maximum/minimum fuel 
limits 278 are used to limit EGT control algorithm fuel command output 276 to protect 
from integrator windup. Resultant output signal 280 is regulated EGT signal fuel flow 
command. In operation, EGT control mode loop 260 operates at about a 100 ms rate. 
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Referring to FIG. 11, speed control mode loop 282 provides various fuel output 
commands to regulate the rotating speed of the turbine. Feedback speed signal 288 is 
read and compared by comparator 286 against setpoint speed signal 284 to determine 
error signal 290, which is then applied to PI algorithm 292 to determine the fuel 
command required to regulate engine speed at the setpoint. EGT control (FIG. 10) and 
maximum/minimum fuel limits are used in conjunction with the speed control 
algorithm 282 to protect output signal 294 from surge and flame out conditions. 
Resultant output signal 298 is regulated turbine speed fuel flow command. In a typical 
implementation, speed control mode loop 282 operates at about a 20ms rate. 

Referring to FIG. 12, power control mode loop 300 regulates the power 
producing potential of the turbine. Feedback power signal 306 is read and compared 
by comparator 304 against setpoint power signal 302 to determine error signal 308, 
which is then applied to PI algorithm 3 1 0 to determine the speed command required to 
regulate output power at the setpoint. Maximum/minimum speed limits are used to 
limit the power control algorithm speed command output to protect output signal 3 12 
from running into over speed and under speed conditions. Resultant output signal 316 
is regulated power signal turbine speed command. In a typical implementation, the 
maximum operating speed of the turbine is generally 96,000 RPM and the minimum 
operating speed of the turbine is generally 45,000 RPM. The loop operates generally at 
about a 500 ms rate. 

Start Only Battery 

Referring to FIG. 14, energy storage device 470 may be a start only battery. In 
the DC bus voltage control mode, start only battery 470 provides energy to regulate 
voltage to the setpoint command. Main CPU 472 commands the bus voltage to control 
at different values depending on the configuration of power controller 478. In the state 
of charge (SOC) control mode, the start only battery system provides a recharging 
power demand when requested. Available recharging power is generally equivalent to 
maximum engine power less power being supplied to the output load and system 
parasitic loads. Main CPU 472 transmits a recharging power level that is the minimum 
of the original power demand and available recharging power. 
Transient Battery 
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The transient battery provides the DC bus voltage control as described below as 
well as the state of charge (SOC) control mode described for the start only battery. The 
transient battery contains a larger energy storage device than the start only battery. 
DC Bus Voltage Control 

DC bus 462 supplies power for logic power, external components and system 
power output. TABLE 1 defines the setpoint the bus voltage is to be controlled at 
based on the output power configuration of power controller 478: 



TABLE 1 

10 B3 POWER OUTPUT SETPOINT 

480/400 VAC Output 800 Vdc 

240/208 VAC Output 400 Vdc 

In the various operating modes, power controller 478 will have different control 
1 5 algorithms responsible for managing the DC bus voltage level. Any of the battery 

options 470 as well as SPs 456 and 458 have modes that control power flow to regulate 
the voltage level of DC bus 462. Under any operating circumstances, only one device 
is commanded to a mode that regulates DC bus 462. Multiple algorithms would 
require sharing logic that would inevitably make system response slower and software 
20 more difficult to comprehend. 
System States 

Referring to FIG. 13, state diagram 320 showing various operating states of 
power controller 478 is illustrated. Sequencing the system through the entire operating 
procedure requires power controller to transition through the operating states defined in 
25 TABLE 2. 



30 



STATE 
# 



0 
1 



SYSTEM 
STATE 



TABLE 2 

DESCRIPTION 



Power Up Performs activities of initializing and testing the system. 
Stand By Closer power to bus and continues system monitoring 
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10 6 



15 



10 



20 11 



25 



30 



12 



13 



14 



15 



while waiting for a start command. 
Prepare to Initializes any external devices preparing for the start 
Start procedure. 

Bearing Configures the system and commands the engine to be 
Lift Off rotated to a predetermined RPM, such as 25,000 RPM. 
Open Loop Turns on ignition system and commands fuel open loop 
Light Off to light the engine. 

Closed Loop Continues motoring and closed fuel control until the 

Acceleration system reaches the no load state. 

Run Engine operates in a no load self-sustaining state 

producing power only to operate the controller. 
Load Converter output contactor is closed and system is 

producing power. 
Re-Charge System operates off of fuel only and produces power 

for recharging energy storage device if installed. 
Cooldown System is motoring engine to reduce EGT before 

shutting down. 

Re-Start Reduces engine speed to begin open loop light when a 
start command is received in the cooldown state. 

Re-Light Performs a turbine re-light in transition from the 

cooldown to warmdown state. Allows continued engine 
cooling when motoring is no longer possible. 

Warmdown Sustains turbine operation with fuel at a predetermined 
RPM, such as 50,000 RPM, to cool when engine 
motoring is not possible. 

Reconfigures the system after a cooldown to enter the 
stand by state. 

Turns off all outputs when presence of fault which 
disables power conversion exists. Logic power is still 
available for interrogating system faults. 
Disable Fault has occurred where processing may no longer 



Shutdown 



Fault 
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f 



be possible. All system operation is disabled. 
Main CPU 472 begins execution in the "power up" state 322 after power is 
applied. Transition to the "stand by" state 324 is performed upon successfully 
completing the tasks of the "power up" state 322. Initiating a start cycle transitions the 
5 system to the "prepare to start" state 326 where all system components are initialized 
for an engine start. The engine then sequences through start states and onto the 
"run/load" state 328. To shutdown the system, a stop command which sends the 
system into either "warm down" or "cool down" state 332 is initiated. Systems that 
have a battery may enter the "re-charge" state 334 prior to entering the "warm down" 

10 or "cool down" state 332. When the system has finally completed the "warm down" or 
"cool down" process 332, a transition through the "shut down" state 330 will be made 
before the system re-enters the "standby" state 324 awaiting the next start cycle. 
During any state, detection of a fault with a system severity level indicating the system 
should not be operated will transition the system state to "fault" state 334. Detection of 

15 faults that indicate a processor failure has occurred will transition the system to the 
"disable" state 336. 

One skilled in the art will recognize that in order to accommodate each mode of 
operation, the state diagram is multidimensional to provide a unique state for each 
operating mode. For example, in the "prepare to start" state 326, control requirements 

20 will vary depending on the selected operating mode. Therefore, the presence of a 
stand-alone "prepare to start" state 326, stand-alone transient "prepare to start" state 
326, utility grid connect "prepare to start" state 326 and utility grid connect transient 
"prepare to start" state 326 will be required. Each combination is known as a system 
configuration (SYSCON) sequence. Main CPU 472 identifies each of the different 

25 system configuration sequences in a 1 6-bit word known as a SYSCON word, which is 
a bit- wise construction of an operating mode and system state number. In a typical 
configuration, the system state number is packed in bits 0 through 1 1 . The operating 
mode number is packed in bits 12 through 15. This packing method provides the 
system with the capability of sequence through 4096 different system states in 16 

30 different operating modes. 
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Separate "power up" 322, "re-light" 338, "warm down" 332, "fault" 334 and 
"disable" 336 states are not required for each mode of operation. The contents of these 
states are mode independent. 

" Power Up" State 

Operation of the system begins in the "power up" state 322 once application of 
power activates main CPU 472. Once power is applied to power controller 478, all the 
hardware components will be automatically reset by hardware circuitry. Main CPU 
472 is responsible for ensuring the hardware is functioning correctly and configure the 
components for operation. Main CPU 472 also initializes its own internal data 
structures and begins execution by starting the Real-Time Operating System (RTOS). 
Successful completion of these tasks directs transition of the software to the "stand by" 
state 324. Main CPU 472 performs these procedures in the following order: 



1. 


Initialize main CPU 472 


2. 


Perform RAM Test 


3. 


Perform FLASH Checksum 


4. 


Start RTOS 


5. 


Run Remaining POST 


6. 


Initialize SPI Communications 


7. 


Verify Generator SP Checksum 


8. 


Verify Converter SP Checksum 


9. 


Initialize IntraController Communications 


10. 


Resolve External Device Addresses 


11. 


Look at Input Line Voltage 


12. 


Determine Mode 


13. 


Initialize Maintenance Port 


14. 


Initialize User Port 


15. 


Initialize External Option Port 


16. 


Initialize InterController 


17. 


Chose Master/Co-Master 


18. 


Resolve Addressing 


19. 


Transition to Stand By State (depends on operating mode) 
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" Stand By" State 

Main CPU 472 continues to perform normal system monitoring in the "stand 
by" state 324 while it waits for a start command signal. Main CPU 472 commands 
either energy storage device 470 or utility 468 to provide continuous power supply. In 
5 operation, main CPU 472 will often be left powered on waiting to be start or for 

troubleshooting purposes. While main CPU 472 is powered up, the software continues 
to monitor the system and perform diagnostics in case any failures should occur. All 
communications will continue to operate providing interface to external sources. A 
start command will transition the system to the "prepare to start" state 326. 
10 " Prepare to Start" State 

Main CPU 472 prepares the control system components for the engine start 
process. Many external devices may require additional time for hardware initialization 
before the actual start procedure can commence. The "prepare to start" state 326 
provides those devices the necessary time to perform initialization and send 
1 5 acknowledgment to the main CPU 472 that the start process can begin. Once also 

systems are ready to go, the software shall transition to the "bearing lift off state 328. 
" Bearing Lift Off ' State 

Main CPU 472 commands generator SP 456 to motor the engine 454 from 
typically about 0 to 25,000 RPM to accomplish the bearing lift off procedure. A check 
20 is performed to ensure the shaft is rotating before transition to the next state occurs. 
" Open Loop Light Off ' State 

Once the motor 452 reaches its liftoff speed, the software commences and 
ensures combustion is occurring in the turbine. In a typical configuration, main CPU 
472 commands generator SP 456 to motor the engine 454 to a dwell speed of about 
25 25,000 RPM. Execution of the open loop light off state 340 starts combustion. Main 
CPU 472 then verifies that the engine 454 has not met the "fail to light" criteria before 
transition to the "closed loop accel" state 342. 
" Closed Loop Accel" State 

Main CPU 472 sequences engine 454 through a combustion heating process to 
30 bring the engine 454 to a self-sustaining operating point. In a typical configuration, 
commands are provided to generator SP 456 commanding an increase in engine speed 
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to about 45,000 RPM at a rate of about 4000 RPM/sec. Fuel controls are executed to 
provide combustion and engine heating. When engine 454 reaches "no load" (requires 
no electrical power to motor), the software transitions to "run" state 344. 
" Run" State 

5 Main CPU 472 continues operation of control algorithms to operate the engine 

at no load. Power may be produced from engine 454 for operating control electronics 
and recharging any energy storage device 470 for starting. No power is output from 
load converter 458. A power enable signal transitions the software into "load" state 
346. A stop command transitions the system to begin shutdown procedures (may vary 
1 0 depending on operating mode), 
" Load" State 

Main CPU 472 continues operation of control algorithms to operate the engine 
454 at the desired load. Load commands are issued through the communications ports, 
display or system loads. A stop command transitions main CPU 472 to begin 
1 5 shutdown procedures (may vary depending on operating mode). A power disable 
signal can transition main CPU 472 back to "run" state 344. 
" Re-charge" State 

Systems that have an energy storage option may be required to charge energy 
storage device 470 to maximum capacity before entering the "warmdown" 348 or 

20 "cooldown" 332 states. During the "re-charge" state 334 of operation, main CPU 472 
continues operation of the turbine producing power for battery charging and controller 
supply. No out power is provided. When the energy storage device 470 has charged, 
the system transitions to either the "cooldown" 332 or "warmdown" 348 state 
depending on system fault conditions. 

25 " Cool Down" State 

"Cool down" state 332 provides the ability to cool the turbine after operation 
and a means of purging fuel from the combustor. After normal operation, software 
sequences the system into "cool down" state 332. In a typical configuration, engine 
454 is motored to a cool down speed of about 45,000 RPM. Airflow continues to 

30 move through engine 454 preventing hot air from migrating to mechanical components 
in the cold section. This motoring process continues until the engine EGT falls below 
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a cool down temperature of about 193°C (380°F). Cool down may be entered at much 
lower than the final cool down temperature when engine 454 fails to light. The 
engine's combustor requires purging of excess fuel which may remain. The software 
always operates the cool down cycle for a minimum purge time of 60 seconds. This 
purge time ensures remaining fuel is evacuated from the combustor. Completion of 
this process transitions the system into the "shutdown" state 330. For user 
convenience, the system does not require a completion of the enter "cooldown" state 
332 before being able to attempt a restart. Issuing a start command transitions the 
system into the "restart" state 350. 
" Restart" State 

Engine 454 is configured from the "cool down" state 332 before engine 454 can 
be restart. In a typical configuration, the software lowers the engine speed to about 
25,000 RPM at a rate of 4,000 RPM/sec. Once the engine speed has reached this level, 
the software transitions the system into the "open loop light off state to perform the 
actual engine start. 

" Shutdown" State 

During the "shutdown" state 330, the engine rotor is brought to rest and system 
outputs are configured for idle operation. In a typical configuration, the software 
commands the rotor to rest by lowering the engine speed at a rate of 2,000 RPM/sec or 
no load condition, whichever is faster. Once the speed reaches about 14,000 RPM, the 
generator SP is commanded to reduce the shaft speed to about 0 RPM in less than 1 
second. 

" Re-light" State 

When a system fault occurs where no power is provided from the utility or 
energy storage device 470, the software re-ignites combustion to perform a warm 
down. The generator SP is configured to regulate voltage (power) for the internal DC 
bus. Fuel is added as defined in the open loop light off fuel control algorithm to ensure 
combustion occurs. Detection of engine light will transition the system to "warm 
down" state 348. 

"Warm Down" State 
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Fuel is provided when no electric power is available to operate engine 454 at a 
no load condition to lower the operating temperature in "warm down" state 348. In a 
typical configuration, engine speed is operated at about 50,000 RPM by supplying fuel 
through the speed control algorithm. Engine temperatures less than about 343°C 
5 (650°F) causes the system to transition to "shutdown" state 330. 
" Fault" State 

The present invention disables all outputs placing the system in a safe 
configuration when faults that prohibit safe operation of the turbine system are present. 
Operation of system monitoring and communications will continue if the energy is 
10 available. 

" Disable" State 

The system disables all outputs placing the system in a safe configuration when 
faults that prohibit safe operation of the turbine system are present. System monitoring 
and communications will most likely not continue. 

15 Modes of Operation 

The turbine works in two major modes - utility grid-connect and stand-alone. 
In the utility grid-connect mode, the electric power distribution system i.e., the utility 
grid, supplies a reference voltage and phase, and the turbine supplies power in 
synchronism with the utility grid. In the stand-alone mode, the turbine supplies its own 

20 reference voltage and phase, and supplies power directly to the load. The power 
controller switches automatically between the modes. 

Within the two major modes of operation are sub-modes. These modes include 
stand-alone black start, stand-alone transient, utility grid connect and utility grid 
connect transient. The criteria for selecting an operating mode is based on numerous 

25 factors, including but not limited to, the presence of voltage on the output terminals, 
the black start battery option, and the transient battery option. 

Referring to FIG. 14, generator converter 456 and load converter 458 provide 
an interface for energy source 460 and utility 468, respectively, to DC bus 462. For 
illustrative purposes, energy source 460 is a turbine including engine 454 and generator 

30 452. Fuel device 474 provides fuel via fuel line 476 to engine 454. Generator 

converter 456 and load converter 458 operate as customized bi-directional switching 
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converters under the control of controller 472. In particular, controller 472 
reconfigures the generator converter 456 and load converter 458 into different 
configurations to provide for the various modes of operation. These modes include 
stand-alone black start, stand-alone transient, utility grid connect and utility grid 
connect transient as discussed in detail below. Controller 472 controls the way in 
which generator 452 and utility 468 sinks or sources power, and DC bus 462 is 
regulated at any time. In this way, energy source 460, utility/load 468 and energy 
storage device 470 can be used to supply, store and/or use power in an efficient 
manner. Controller 472 provides command signals via line 479 to engine 454 to 
determine the speed of turbine 460. The speed of turbine 460 is maintained through 
generator 452. Controller 472 also provides command signals via control line 480 to 
fuel device 474 to maintain the EGT of the engine 454 at its maximum efficiency 
point. Generator SP 456 is responsible for maintaining the speed of the turbine 460, 
but putting current into generator 452 or pulling current out of generator 452. 
Stand-alone Black Start 

Referring to FIG. 14, in the stand-alone black start mode, energy storage device 
470, such as battery, is provided for starting purposes while energy source 460, such as 
turbine including engine 454 and generator 452, supplies all transient and steady state 
energy. Referring to TABLE 3, controls for a typical stand-alone black start mode are 
shown. 

TABLE 3 



SYSTEM 
STATE 



ENGINE MOTOR CONVERTER ENERGY STORAGE 

CONTROLS CONTROLS CONTROLS CONTROLS 



Power Up 
Stand By 
Prepare to Start 

Bearing Lift Off - RPM 

Open Loop Light Off Open Loop RPM 

Light 

Closed Loop Accel EGT RPM 

Run Speed DC Bus 

Load Speed DC Bus 



Voltage 



DC Bus 
DC Bus 
DC Bus 
DC Bus 

DC Bus 

SOC 

SOC 
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Recharge 
Cool Down 
Restart 
Shutdown 
Re-light 
Warm Down 
Fault 
Disable 



Speed 



Speed 
Speed 



DC Bus 
RPM 
RPM 
RPM 
DC Bus 
DC Bus 



SOC 
DC Bus 
DC Bus 
DC Bus 



10 



15 



20 



25 



30 



Stand-alone Transient 

In the stand-alone transient mode, storage device 479 is provided for the 
purpose of starting and assisting the energy source 460, in this example the turbine, to 
supply maximum rated output power during transient conditions. Storage device 479, 
typically a battery, is always attached to DC bus 462 during operation, supplying 
energy in the form of current to maintain the voltage on DC bus 462. Converter/SP 
458 provides a constant voltage source when producing output power. As a result, 
load 468 is always supplied the proper AC voltage value that it requires. Referring to 
TABLE 4, controls for a typical stand-alone transient mode are shown. 

TABLE 4 



SYSTEM 
STATE 



ENGINE MOTOR CONVERTER ENERGY STORAGE 

CONTROLS CONTROLS CONTROLS CONTROLS 



Power Up 
Stand By 
Prepare to Start 
Bearing Lift Off 
Open Loop Light Off 

Closed Loop Accel 

Run 

Load 

Recharge 

Cool Down 

Restart 



RPM 

Open Loop RPM 
Light 

EGT RPM 
Power &EGT RPM 
Power & EGT RPM 
Power & EGT RPM 
RPM 
RPM 



Voltage 



DC Bus 
DC Bus 
DC Bus 
DC Bus 

DC Bus 
DC Bus 
DC Bus 
DC Bus 
DC Bus 
DC Bus 
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Shutdown - RPM - DC Bus 

Re-light Speed DC Bus 

Warm Down Speed DC Bus 

Fault - 

5 Disable - 



Utility Grid Connect 

Referring to FIG. 14, in the utility grid connect mode, the energy source 460, in 
this example the turbine is connected to the utility grid 468 providing load leveling and 
1 0 management where transients are handled by the utility grid 468. The system operates 
as a current source, pumping current into utility 468. Referring to TABLE 5, controls 
for a typical utility grid connect mode are shown. 



TABLE 5 



15 



20 



25 



30 



SYSTEM 
STATE 


ENGINE 
CONTROLS 


MOTOR 
CONTROLS 


CONTROLS 


CON 


Power Up 








N/A 


Stand By 








N/A 


Prepare to Start 






DC Bus 


N/A 


Bearing Lift Off 




RPM 


DC Bus 


N/A 


Open Loop Light Off 


Open Loop 
Light 


RPM 


DC Bus 


N/A 


Closed Loop Accel 


EGT 


RPM 


DC Bus 


N/A 


Run 


Power & EGT 


RPM 


DC Bus . 


N/A 


Load 


Power & EGT 


RPM 


DC Bus 


N/A 


Recharge 


N/A 


N/A 


N/A 


N/A 


Cool Down 




RPM 


DC Bus 


N/A 


Restart 




RPM 


DC Bus 


N/A 


Shutdown 




RPM 


DC Bus 


N/A 


Re-light 


Speed 


DC Bus 




N/A 


Warm Down 


Speed 


DC Bus 




N/A 


Fault 








N/A 


Disable 








N/A 
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Utility Grid Connect Transient 

In the utility grid connect transient mode, the energy source 460, in this 
example the turbine, is connected to the utility grid 468 providing load leveling and 
5 management. The turbine that is assisted by energy storage device 470, typically a 
battery, handles transients. The system operates as a current source, pumping current 
into utility 468 with the assistance of energy storage device 470. Referring to TABLE 
6, controls for a typical utility grid connect transient mode are shown. 



10 



TABLE 6 



SYSTEM 
STATE 



ENGINE MOTOR CONVERTER ENERGY STORAGE 

CONTROLS CONTROLS CONTROLS CONTROLS 



Power Up 
15 Stand By 

Prepare to Start 
Bearing Lift Off 
Open Loop Light Off 

20 Closed Loop Accel 

Run 

Load 

Recharge 

Cool Down 
25 Restart 

Shutdown 

Re-light 

Warm Down 

Fault 
30 Disable 



RPM 

Open Loop RPM 
Light 

EGT RPM 
Power &EGT RPM 
Power & EGT RPM 
Power & EGT RPM 
RPM 
RPM 
RPM 

Speed DC Bus 

Speed DC Bus 
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Multi-pack Operation 

In accordance with the present invention, the power controller can operate in a 
single or multi-pack configuration. In particular, power controller, in addition to being 



96691.01 01 



-30- 



153501-0053 

a controller for a single turbogenerator, is capable of sequencing multiple systems as 
well. Referring to FIG. 15, for illustrative purposes, multi-pack system 510 including 
three power controllers 518, 520 and 522 is shown. The ability to control multiple 
controllers 518, 520 and 522 is made possible through digital communications 
interface and control logic contained in each controllers main CPU (not shown). 

Two communications busses 530 and 534 are used to create the intercontroller 
digital communications interface for multi-pack operation. One bus 534 is used for 
slower data exchange while the other bus 530 generates synchronization packets at a 
faster rate. In a typical implementation, for example, an IEEE-502.3 bus links each of 
the controllers 518, 520 and 522 together for slower communications including data 
acquisition, start, stop, power demand and mode selection functionality. An RS485 
bus links each of the systems together providing synchronization of the output power 
waveforms. 

One skilled in the art will recognize that the number of power controllers that 
can be connected together is not limited to three, but rather any number of controllers 
can be connected together in a multi-pack configuration. Each power controller 5 1 8, 
520 and 522 includes its own energy storage device 524, 526 and 528, respectively, 
such as a battery. In accordance with another embodiment of the invention, power 
controllers 518, 520 and 522 can all be connected to the same single energy storage 
device (not shown), typically a very large energy storage device which would be rated 
too big for an individual turbine. Distribution panel, typically comprised of circuit 
breakers, provides for distribution of energy. 

Multi-pack control logic determines at power up that one controller is the 
master and the other controllers become slave devices. The master is in charge of 
handling all user-input commands, initiating all inter-system communications 
transactions, and dispatching units. While all controllers 518, 520 and 522 contain the 
functionality to be a master, to alleviate control and bus contention, one controller is 
designated as the master. 

At power up, the individual controllers 518, 520 and 522 determine what 
external input devices they have connected. When a controller contains a minimum 
number of input devices it sends a transmission on intercontroller bus 530 claiming to 
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be master. All controllers 5 1 8, 520 and 522 claiming to be a master begin resolving 
who should be master. Once a master is chosen, an address resolution protocol is 
executed to assign addresses to each slave system. After choosing the master and 
assigning slave addresses, multi-pack system 510 can begin operating. 
5 A co-master is also selected during the master and address resolution cycle. 

The job of the co-master is to act like a slave during normal operations. The co-master 
should receive a constant transmission packet from the master indicating that the 
master is still operating correctly. When this packet is not received within a safe time 
period, 20 ms for example, the co-master may immediately become the master and take 

1 0 over master control responsibilities. 

Logic in the master configures all slave turbogenerator systems. Slaves are 
selected to be either utility grid-connect (current source) or standalone (voltage source). 
A master controller, when selected, will communicate with its output converter logic 
(converter SP) that this system is a master. The converter SP is then responsible for 

1 5 transmitting packets over the intercontroller bus 530, synchronizing the output 

waveforms with all slave systems. Transmitted packets will include at least the angle 
of the output waveform and error-checking information with transmission expected 
every quarter cycle to one cycle. 

Master control logic will dispatch units based on one of three modes of 

20 operation: (1) peak shaving, (2) load following, or (3) base load. Peak shaving 
measures the total power consumption in a building or application using a power 
meter, and the multi-pack system 510 reduces the utility consumption of a fixed load, 
thereby' reducing the utility rate schedule and increasing the overall economic return of 
the turbogenerator. Load following is a subset of peak shaving where a power meter 

25 measures the total power consumption in a building or application and the multi-pack 
system 10 reduces the utility consumption to zero load. In base load, the multi-pack 
system 10 provides a fixed load and the utility supplements the load in a building or 
application. Each of these control modes require different control strategies to 
optimize the total operating efficiency. 

30 A minimum number of input devices are typically desired for a system 5 10 to 

claim it is a master during the master resolution process. Input devices that are looked 
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for include a display panel, an active RS232 connection and a power meter connected 
to the option port. Multi-pack system 5 1 0 typically requires a display panel or RS232 
connection for receiving user-input commands and power meter for load following or 
peak shaving. 

In accordance with the present invention, the master control logic dispatches 
controllers based on operating time. This would involve turning off controllers that 
have been operating for long periods of time and turning on controllers with less 
operating time, thereby reducing wear on specific systems. 
Utility Grid Analysis and Transient Ride Through 

Referring to FIGS. 16-18, transient handling system 580 for power controller 
620 is illustrated. Transient handling system 580 allows power controller 620 to ride 
through transients which are associated with switching of correction capacitors on 
utility grid 616 which causes voltage spikes followed by ringing. Transient handling 
system 580 also allows ride through of other faults, including but not limited to, short 
circuit faults on utility grid 616, which cleared successfully, cause voltage sags. 
Transient handling system 580 is particularly effective towards handling transients 
associated with digital controllers, which generally have a slower current response rate 
due to A/D conversion sampling. During a transient, a large change in the current can 
occur in between A/D conversions. The high voltage impulse caused by transients 
typically causes an over current in digital power controllers. 

As is illustrated in FIG. 17, a graph 590 showing transients typically present on 
utility grid 616 is shown. The duration of a voltage transient, measured in seconds, is 
shown on the x-axis and its magnitude, measured in volts, is shown on the y-axis. A 
capacitor switching transient, such as shown at 592, which is relatively high in 
magnitude (up to about 200%) and short in duration (somewhere between 1 and 20 
milliseconds) could be problematic to operation of a power controller. 

Referring to FIGS. 16-18, changes on utility grid 616 are reflected as changes 
in the magnitude of the voltage. In particular, the type and seriousness of any fault or 
event on utility grid 616 can be determined by magnitude estimator 584, which 
monitors the magnitude and duration of any change on utility grid 616. 
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In accordance with the present invention, the effect of voltage transients can be 
minimized by monitoring the current such that when it exceeds a predetermined level , 
switching is stopped so that the current can decay, thereby preventing the current from 
exceeding its predetermined level. The present invention thus takes advantage of 
5 analog over current detection circuits that have a faster response than transient 

detection based on digital sampling of current and voltage. Longer duration transients 
indicate abnormal utility grid conditions. These must be detected so power controller 
. 620 can shut down in a safe manner. In accordance with the present invention, 
algorithms used to operate power controller 620 provide protection against islanding of 
10 power controller 620 in the absence of utility-supplied grid voltage. Near short or near 
open islands are detected within milliseconds through loss of current control. Islands 
whose load is more closely matched to the power controller output will be detected 
through abnormal voltage magnitudes and frequencies as detected by magnitude 
estimator 584. 

1 5 In particular, referring to FIG. 1 8, power controller 620 includes brake resistor 

612 connected across DC bus 622. Brake resistor 612 acts as a resistive load, 
absorbing energy when converter SP 608 is turned off In operation, when converter 
SP 608 is turned off, power is no longer exchanged with utility grid 616, but power is 
still being received from the turbine, which is absorbed by brake resistor 612. The 

20 present invention detects the DC voltage between generator and converter SPs 606 and 
608. When the voltage starts to rise, brake resistor 612 is turned on to allow it to 
absorb energy. 

In a typical configuration, AC motor 618 produces three phases of AC at 
variable frequencies. AC/DC converter 602 under the control of motor SP 606 

25 converts the AC to DC which is then applied to DC bus 622 (regulated for example at 
800 vDC) which is supported by capacitor 610 (for example, at 800 microfarads with 
two milliseconds of energy storage). AC/DC converter 604, under the control of 
converter SP 608, converts the DC into three-phase AC, and applies it to utility grid 
616. In accordance with the present invention, current from DC bus 622 can by 

30 dissipated in brake resistor 612 via modulation of switch 614 operating under the 

control of motor SP 606. Switch 614 may be an IGBT switch, although one skilled in 
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the art will recognize that other conventional or newly developed switches may be 
utilized as well. 

Motor SP 606 controls switch 614 in accordance to the magnitude of the 
voltage on DC bus 622. The bus voltage of DC bus 622 is typically maintained by 
5 converter SP 608, which shuttles power in and out of utility grid 616 to keep DC bus 
622 regulated at, for example, 800 vDC. When converter SP 608 is turned off, it no 
longer is able to maintain the voltage of DC bus 622, so power coming in from the 
motor causes bus voltage of DC bus 622 to rise quickly. The rise in voltage is detected 
by motor SP 606, which turns on brake resistor 612 and modulates it on and off until 
10 the bus voltage is restored to its desired voltage, for example, 800 vDC. Converter SP 
608 detects when the utility grid transient has dissipated, i.e., AC current has decayed 
to zero 

and restarts the converter side of power controller 620. Brake resistor 612 is sized so 
that it can ride through the transient and the time taken to restart converter. 

1 5 Referring to FIGS. 1 8 and 20, in accordance with the present invention, both 

the voltage and zero crossings (to determine where the AC waveform of utility grid 
616 crosses zero) are monitored to provide an accurate model of utility grid 616. 
Utility grid analysis system includes angle estimator 582, magnitude estimator 584 and 
phase locked loop 586. The present invention continuously monitors utility grid 

20 voltage and based on these measurements, estimates the utility grid angle, thus 

facilitating recognition of under/over voltages and sudden transients. Current limits 
are set to disable DC/AC converter 604 when current exceeds a maximum and wait 
until current decays to an acceptable level. The result of measuring the current and 
cutting it off is to allow DC/AC converter 604 to ride through transients better. Thus 

25 when DC/ AC converter 608 is no longer exchanging power with utility grid 616, 
power is dissipated in brake resistor 612, 

In accordance with the present invention, converter SP 608 is capable of 
monitoring the voltage and current at utility grid 616 simultaneously. In particular, 
power controller 620 includes a utility grid analysis algorithm. One skilled in the art 

30 will recognize that estimates of the utility grid angle and magnitude may be derived via 
conventional algorithms or means. The true utility grid angle $xc> which is the angle 
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of the generating source, cycles through from 0 to 2n and back to 0 at a rate of 60 
hertz. The voltage magnitude estimates of the three phases are designated Vi mag; V 2 
mag and V 3 mag and the voltage measurement of the three phases are designated Vj, V 2 
and V3. 

5 A waveform, constructed based upon the estimates of the magnitude and angle 

for each phase, indicates what a correct measurement would look like. For example, 
using the first of the three phase voltages, the cosine of the true utility grid angle 6 A c is 
multiplied by the voltage magnitude estimate Vj mag , with the product being a cosine- 
like waveform. Ideally, the product would be voltage measurement Vi. 

10 Feedback loop 588 uses the difference between the absolute magnitude of the 

measurement of Vi and of the constructed waveform to adjusts the magnitude of the 
magnitude estimate Vi mag . One skilled in the art will recognize that the other two 
phases of three-phase signal can be adjusted similarly, with different angle templates 
corresponding to different phases of the signal. Thus, magnitude estimate Vi mag and 

1 5 angle estimate 0est are used to update magnitude estimate Y\ mag . Voltage magnitude 
estimates Vi mag , V 2 mag and V 3 mag are steady state values used in a feedback 
configuration to track the magnitude of voltage measurements Vi, V 2 and V 3 . By 
dividing the measured voltages Vi by the estimates of the magnitude Vi mag , the cosine 
of the angle for the first phase can be determined (similarly, the cosine of the angles of 

20 the other signals will be similarly determined). 

In accordance with the present invention, the most advantageous estimate for 
the cosine of the angle, generally the one that is changing the most rapidly, is chosen to 
determine the instantaneous measured angle. In most cases, the phase that has an 
estimate for the cosine of an angle closest to zero is selected since it yields the greatest 

25 accuracy. Utility grid analysis system 580 thus includes logic to select which one of 
the cosines to use. The angle chosen is applied to angle estimator 582, from which an 
estimate of the instantaneous angle 9 E st of utility grid 616 is calculated and applied to 
phase locked loop 586 to produce a filtered frequency. The angle is thus differentiated 
to form a frequency that is then passed through a low pass filter (not shown). Phase 

30 locked loop 586 integrates the frequency and also locks the phase of the estimated 
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instantaneous angle e ES T , which may have changed in phase due to differentiation and 
integration, to the phase of true utility grid angle 6 A c- 

In a typical operation, when the phase changes suddenly on measured voltage 
Vi, the algorithm of the present invention compares the product of the magnitude 
estimate Vi mag and the cosine of true utility grid angle 9 A c against the real magnitude 
multiplied by the cosine of a different angle. A sudden jump in magnitude would be 
realized. 

Thus, three reasonably constant DC voltage magnitude estimates are generated. 
A change in one of those voltages indicates whether the transient present on utility grid 
616 is substantial or not. One skilled in the art will recognize that there are a number 
of ways to determine whether a transient is substantial or not, i.e. whether abnormal 
conditions exist on the utility grid system, which require power controller 620 to shut 
down. A transient can be deemed substantial based upon the size of the voltage 
magnitude and duration. Examples of the criteria for shutting down power controller 
620 are shown in FIG. 1 7. Detection of abnormal utility grid behavior can also be 
determined by examining the frequency estimate. 

On detecting abnormal utility grid behavior, a utility grid fault shutdown is 
initiated. When system controller 620 initiates a utility grid fault shutdown, output 
contactor is opened within a predetermined period of time, for example, 100 msec, and 
the main fuel trip solenoid (not shown) is closed, removing fuel from the 
turbogenerator. A warm shutdown ensues during which control power is supplied from 
generator 618 as it slows down. In a typical configuration, the warm-down lasts about 
1-2 minutes before the rotor (not shown) is stopped. The control software does not 
allow a restart until utility grid voltage and frequency are within permitted limits. 

Having now described the invention in accordance with the requirements of the 
patent statutes, those skilled in this art will understand how to make changes and 
modifications in the present invention to meet their specific requirements or conditions. 
For example, the power controller, while described generally, may be implemented in 
an analog or digital configuration. In the preferred digital configuration, one skilled in 
the art will recognize that various terms utilized in the invention are generic to both 
analog and digital configurations of power controller. For example, converters 
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referenced in the present application is a general term which includes inverters, signal 
processors referenced in the present application is a general term which includes digital 
signal processors, and so forth. Correspondingly, in a digital implementation of the 
present invention, inverters and digital signal processors would be utilized. Such 
changes and modifications may be made without departing from the scope and spirit of 
the invention as set forth in the following claims. 
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ABSTRACT OF THE DISCLOSURE 
POWER CONTROLLER 

A power controller provides a distributed generation power networking system in 
which bi-directional power converters are used with a common DC bus for permitting 
compatibility between various energy components. Each power converter operates 
essentially as a customized bi-directional switching converter configured, under the 
5 control of the power controller, to provide an interface for a specific energy component to 
the DC bus. The power controller controls the way in which each energy component, at 
any moment, will sink or source power, and the manner in which the DC bus is regulated. 
In this way, various energy components can be used to supply, store and/or use power in 
an efficient manner. The various energy components include energy sources, loads, 
1 0 storage devices and combinations thereof. 
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WHAT IS CLAIMED IS: 



1 1 . A power controller for distributing power among a plurality of energy 

2 components, comprising: 

3 a DC bus; and 

4 a plurality of power converters, each of which is connected between one of said 

5 energy components and said DC bus and is responsive to said power controller, wherein 

6 said power controller provides a distributed generation power system by controlling the 

7 way each energy component sinks or sources power and said DC bus is regulated. 

1 2. A power controller claimed in claim 1 , wherein each of said power converters 

2 operates as a customized bi-directional switching converter configured, under the control 

3 of said power controller, to provide an interface for said energy component to said DC 

4 bus. 

1 3. The power controller claimed in claim 1 , wherein each of said power converters 

2 comprises: 

3 a power switching system; and 

4 a processing system for providing control to said power switching system. 

1 4. The power controller claimed in claim 3, wherein said processing system further 

2 comprises: 

3 a signal processor; and 

4 a central processing unit for providing control to said signal processor. 

1 5. The power controller claimed in claim 3, wherein said central processing unit 

2 reconfigures said power converter into different configurations for different modes of 

3 operation. 

1 6. The power controller claimed in claim 3, wherein said power switching system 

2 comprises a plurality of insulated gate bipolar transistor switches. 

1 7. The power controller claimed in claim 1 , wherein said plurality of energy 

2 components includes an energy source. 
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1 8. The power controller claimed in claim 7, wherein said energy source comprises a 

2 gas turbine. 

1 9. The power controller claimed in claim 8, wherein said gas turbine drives an AC 

2 generator to produce AC which is applied to said power converter. 

1 1 0. The power controller claimed in claim 9, wherein said controller regulates DC bus 

2 voltage independently of turbine speed. 

1 11. The power controller claimed in claim 1 , wherein said plurality of energy 

2 components includes an energy storage device. 

1 12. The power controller claimed in claim 1 1 , wherein said energy storage device 

2 comprises a flywheel, 

1 13. The power controller claimed in claim 1 1, wherein said energy storage device 

2 comprises a battery, 

1 14. The power controller claimed in claim 1 1 , wherein said energy storage device 

2 comprises an ultracap. 

1 15, The power controller claimed in claim 1 , wherein said plurality of energy 

2 components includes a load. 

1 1 6. The power controller claimed in claim 1 5, wherein said load comprises an AC 

2 utility. 

1 17. The power controller claimed in claim 16, wherein said load comprises a DC 



1 1 8. The power controller claimed in claim 16, wherein said load comprises a drive 

2 motor. 

1 19. The power controller claimed in claim 1 , wherein said plurality of energy 

2 components includes an energy source, a load and a storage device. 



2 



load. 
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1 20. The power controller claimed in claim 1 9, wherein during a utility start up mode 

2 of operation, one of said power converters applies power from said load to said DC bus 

3 for conversion by another of said power converters into power required by said energy 

4 source to startup and said storage device is disconnected from said DC bus while said 

5 load regulates DC voltage on said DC bus. 

1 21 . The power controller claimed in claim 20, wherein said energy source comprises 

2 a turbine and is controlled in a local feedback loop to maintain said turbine revolutions 

3 per minute (RPM). 

1 22. The power controller claimed in claim 19, wherein during a utility start up mode 

2 of operation, one of said power converters isolates said DC bus so that another of said 

3 power converters provides the required starting power from said DC bus to said energy 

4 source. 

1 23 . The power controller claimed in claim 1 , wherein said power converters comprise 

2 a generator converter and output converter. 

1 24. The power controller claimed in claim 23, wherein^said power controller is 

2 digital, said generator converter is a generator inverter and said output converter is an 

3 output inverter. 

1 25. The power controller claimed in claim 4, wherein said central processing unit 

2 sequences events which occur inside said power controller and arbitrates communications 

3 to externally connected devices. 

1 26. The power controller claimed in claim 4, wherein said central processing unit is a 

2 Motorola MC68332 microprocessor. 

1 27. The power controller claimed in claim 19, wherein said energy source comprises 

2 a turbogenerator, said energy storage device comprises a battery and said load comprises 

3 a utility grid. 

1 28. The power controller claimed in claim 19, wherein said energy source comprises 

2 a turbogenerator, said energy storage device comprises a capacitive device and said load 

3 comprises a utility grid. 

1 29. The power controller claimed in claim 19 wherein said energy source comprises a 

2 turbogenerator, and said load comprises a utility grid. 
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30. The power controller claimed in claim 1, wherein said plurality of energy 
components comprises a turbogenerator and a load. 

31. A turbine system, comprising: 
a turbine engine; 

a load; and 

a power controller for converting electricity from said turbine engine into 
regulated DC and then to AC electricity, wherein said power controller includes an 
engine power conversion in communication with said turbine engine, an utility power 
conversion in communication with said load and a DC bus. 

32. The turbine system claimed in claim 31, further comprising: 

a fuel metering system in communication with an energy reservoir controller and 
said power controller. 

33. The turbine system claimed in claim 3 1 , wherein said power controller provides a 
distributed generation power system utilizing said engine power conversion and said 
utility power conversion. 

34. The turbine system claimed in claim 3 1 , wherein said engine power conversion 
and said utility power conversion operate as a customized bi-directional switching 
converters, under control of said power controller, to provide an interface for said turbine 
engine and said load to said DC bus. 

35. A method for controlling the distribution of power among a plurality of energy 
components, comprising the steps of: 

connecting a power converter between a DC bus and each of said energy 
components; 

controlling the way each of said energy components sinks or sources power; and 
controlling the way said DC bus is regulated responsive to operation of each of 
said energy components. 

36. The method claimed in claim 35, wherein said step of controlling the way each of 
said energy components sinks or sources power, further comprises the step of: 

controlling the way each of said power converters is configured. 

37. The method claimed in claim 35, wherein said step of controlling the way each of 
said energy components sinks or sources power, further comprises the step of: 
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3 controlling said power converters such that they operate as customized, bi- 

4 directional switching converters configured to provide an interface for said energy 

5 components to said DC bus. 

1 38. The method claimed in claim 36, wherein said each of said power converters 

2 comprises a power switching system and said step of controlling the way each of said 

3 power converters is configured, further comprises the step of: 

4 providing control to said power switching system. 

1 39. The method claimed in claim 35, wherein said step of controlling the way each of 

2 said power converters is configured, further comprises the step of: 

3 reconfiguring said each of said power converters into different configurations for 

4 different modes of operation. 

1 40. The method claimed in 35, wherein said plurality of energy components include 

2 an energy source, a load and a storage device. 

1 41 . The method claimed in claim 40, wherein said energy source comprises a 

2 turbogenerator, said energy storage device comprises a battery and said load comprises a 

3 utility grid. 

1 42. The method claimed in claim 40, wherein said energy source comprises a 

2 turbogenerator, said energy storage device comprises a capacitive device and said load 

3 comprises a utility grid. 

1 43. The method claimed in claim 40, further comprising the steps of: 

2 during a utility start up mode of operation, utilizing one of said power converters 

3 to apply power from said load to said DC bus for conversion by another of said power 

4 converters into power required by said energy source to startup; 

5 disconnecting said storage device from said DC bus while said load regulates DC 

6 voltage on said DC bus. 

1 44. The method claimed in claim 43, wherein said energy source comprises a turbine 

2 and further comprising the step of: 

3 utilizing a local feedback loop to maintain said revolutions per minute (RPM) of 

4 said turbine. 

1 45. The method claimed in claim 40 further comprising the steps of: 
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2 during a utility start up mode of operation, isolating one of said power converters 

3 from said DC bus so that another of said power converters provides the required starting 

4 power from said DC bus to said energy source. 

1 46. The method claimed in claim 35, wherein said power converters comprise a 

2 generator converter and output converter. 

1 47. The method claimed in claim 46, wherein said distribution of power is controlled 

2 digitally, said generator converter is a generator inverter and said output converter is an 

3 output inverter. 

1 48. The method claimed in claim 35, further comprising the steps of: 

2 sequencing events which occur; and 

3 arbitrating communications to externally connected devices. 

1 49. A power control system, comprising: 

2 a plurality of power controllers for distributing power among a plurality of energy 

3 components; and 

4 a first bus for providing communication between said plurality of power 

5 controllers. 

1 50. The power control system claimed in claim 49, further comprising: 

2 a second bus for providing communication between said plurality of power 

3 controllers. 

1 51. The power control system claimed in claim 50, wherein said first bus provides for 

2 data exchange at a first rate and said second bus provides for data exchange at a second 

3 rate. 

1 52. The power control system claimed in claim 5 1 , wherein said first bus provides for 

2 sychronization of output power waveforms. 

1 53. The power control system claimed in claim 52, wherein said second bus provides 

2 for data acquisition. 

1 54. The power control system claimed in claim 52, wherein said second bus provides 

2 for start and stop operation of said system. 

1 55. The power control system claimed in claim 52, wherein said second bus provides 

2 mode selection. 
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1 56. The power control system claimed in claim 49, wherein each of said plurality of 

2 power controllers comprises: 

3 a DC bus; and 

4 a plurality of power converters, one of said plurality of power converters 

5 connected between one of said plurality of energy components and said DC bus and is 

6 responsive to said power controller, wherein said power controller provides a distributed 

7 generation power system by controlling the way each energy component sinks or sources 

8 power and said DC bus is regulated. 

1 57. The power control system claimed in claim 52, wherein each of said power 

2 converters operates as a customized bi-directional switching converter configured, under 

3 the control of said power controller, to provide an interface for said energy component to 

4 said DC bus. 

1 58. The power control system claimed in claim 49, wherein said plurality of power 

2 controllers includes a master controller with remainder of said plurality of power 

3 controllers being slave controllers-. 

1 59. The power control system claimed in claim 58, wherein said master controller 

2 handles user-input commands. 

1 60. The power control system claimed in claim 58, wherein said master controller 

2 initiates inter-system communication transactions. 

1 61 . The power control system claimed in claim 58, wherein said master controller 

2 dispatches said slave controllers. 

1 62. The power control system claimed 58, wherein said system determines a master 

2 controller and assigns addresses to said slave controllers. 

1 63. The power control system claimed in claim 58, wherein said slave controllers 

2 includes a co-master controller for monitoring operation of said master controller. 

1 64. The power control system claimed in claim 63, wherein said co-master controller 

2 receives a signal from said master controller indicating that said master controller is 

3 operating correctly. 

1 65. The power control system claimed in claim 58, wherein said slave controllers are 

2 selected to be either a current source or standalon 
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1 66. The power control system claimed in 58, wherein said master controller further 

2 comprises: 

3 a signal processor; and 

4 a central processing unit for providing control to said signal processor. 

1 67. The power control system claimed in claim 58, wherein said signal processor is 

2 digital 

1 68. The power control system claimed in claim 66, wherein said master controller 

2 communicates with said signal processor that it is a master controller. 

1 69. The power control system claimed in claim 66, wherein said signal processor 

2 transmits packets over said first bus for synchronizing output waveforms with said slave 

3 controllers. 

1 70. The power control system claimed in claim 69, wherein said packets include an 

2 angle of said output waveforms. 

1 71 . The power control system claimed in claim 69, wherein said packets include 

2 error-checking information. 

1 72. The power control system claimed in claim 58, wherein said master controller 

2 measures total power consumption and said system, in response to said measurement, 

3 adjusts utility consumption to a fixed load. 

1 73. The power control system claimed in claim 59, wherein said master controller 

2 adjusts utility consumption to a zero load. 

1 74. The power control system claimed in claim 58, wherein said system provides a 

2 fixed load and utility supplements said load in an application. 

1 75. The power control system claimed in claim 58, wherein said master controller 

2 dispatches slave controllers based on operating time. 

1 76. The power control system claimed in claim 75, wherein said master controller 

2 turns off slave controllers that have been operating for long periods of time and turning 

3 on slave controllers with less operating time by reducing wear on specific slave 

4 controllers. 

1 77. The power control system claimed in claim 49, wherein system provides for 

2 multi-turbogenerator control. 

1 78. The power control system claimed in claim 77, further comprising: 
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a plurality of turbine engines; 
a plurality of loads; 

wherein said plurality of power controllers convert electricity from said plurality 
of turbine engines into regulated DC and then to AC electricity, wherein said plurality of 
power controllers include an engine power conversion in communication with said 
plurality of turbine enginee, an utility power conversion in communication with said 
plurality of loads and a DC bus. 

79. A digital power controller for distributing power among a plurality of energy 
components, comprising: 

a DC bus; and 

a plurality of power inverters, each of which is connected between said energy 
component and said DC bus and is responsive to said power controller, wherein said 
power controller provides a distributed generation power system by controlling the way 
each energy component sinks or sources power and said DC bus is regulated, 

wherein said energy components include an energy source, energy storage device 
and load. 

80. The digital power controller claimed in claim 79, wherein said controller operates 
in a first and second mode. 

81 . The digital power controller claimed in claim 80, wherein in said first mode, said 
load supplies a reference voltage and phase, and said energy source supplies power in 
synchronism with said load. 

82. The digital power controller claimed in claim 80, wherein in said second mode, 
said energy source supplies its own reference voltage and phase, and supplies power 
directly to said load. 

83. The digital power controller claimed in claim 79, wherein said energy source 
comprises a turbine and said load comprises a grid. 

84. The digital power controller claimed in claim 80, wherein said controller switches 
automatically between said first and second modes. 

85. The digital power controller claimed in claim 80, wherein each of said first and 
second modes comprises a plurality of sub-modes. 
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1 86. The digital power controller claimed in claim 85, wherein in a first sub-mode, 

2 said energy storage device provides for starting and said energy source supplies transient 

3 and steady state energy. 

1 87. The digital power controller claimed in claim 85, wherein in a second sub-mode, 

2 said storage device starts and assists said energy source to supply maximum output power 

3 during transient conditions. 

1 88. The digital power controller claimed in claim 87, wherein said storage device is 

2 always attached to said DC bus during operation, supplying energy to maintain voltage 

3 on said DC bus. 

1 89. The digital power controller claimed in claim 85, wherein in a third sub-mode, 

2 said energy source is connected to said load providing load leveling and management and 

3 said load handles transients. 

1 90. The digital power controller claimed in claim 85, wherein in a fourth sub-mode, 

2 said energy source is connected to said load providing load leveling and management and 

3 said energy storage device handles transients. 

1 91 . A method of controlling the distribution of power among a plurality of energy 

2 components using a computer including a digital signal processor comprising the steps 

3 of: 

4 interfacing a plurality of power inverters between a DC bus and each of said 

5 energy components; 

6 controlling the way each of said energy components sinks or sources power; and 

7 controlling the way said DC bus is regulated responsive to operation of each of 

8 said energy components, wherein said energy components include an energy source, a 

9 load and a storage device. 

1 92. The method claimed in claim 9 1 , further comprising the steps of: 

2 applying power to said power inverters for start up; 

3 error checking said power inverters; 

4 initializing internal data structures; and 

5 starting operating system. 

1 93. The method claimed in claim 92, further comprising the steps of: 

2 monitoring said system and performing diagnostics should failures occur; and 



96472.01 03 



153501-0053 



3 commanding either energy storage device or load to provide continuous power 

4 supply. 

1 94. The method claimed in claim 93, further comprising the steps of: 

2 initializing external devices; and 

3 acknowledging start process can begin. 

1 95. The method claimed in claim 94, further comprising the steps of: 

2 commanding signal processor to motor said energy source; and 

3 ensuring said energy source is rotating. 

1 96. The method claimed in claim 95, further comprising the step of: 

2 once said energy source reaches a predetermined speed, ensuring combustion is 

3 occurring. 

1 97. The method claimed in claim 96, further comprising the step of: 

2 sequencing said energy source through a heating process to bring said energy 

3 source to a self-sustaining operating point. 

1 98. The method claimed in claim 97, further comprising the step of: 

2 continuing operation of control algorithms to operate said energy source at no 

3 load. 

1 99. The method claimed in claim 98, further comprising the step of: 

2 continuing operation of control algorithms to operate said energy source at a 

3 desired load. 

1 100. The method claimed in claim 99, further comprising the step of: 

2 charging energy storage device to maximum capacity. 

1 101. The method claimed in claim 1 00, further comprising the steps of: 

2 cooling said energy source after operation; and 

3 purging fuel. 

1 1 02. The method claimed in claim 101, further comprising the step of: 

2 configuring said energy source before said energy source is restarted. 

1 103. The method claimed in claim 1 02, further comprising the steps of: 

2 resting said energy source; and 

3 configuring system outputs for idle operation. 
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1 1 04. The method claimed in claim 102, further comprising the step of: 

2 re-igniting combustion to perform a warm down when a system fault occurs 

3 where no power is provided from load or energy storage device. 

1 1 05. The method claimed in claim 1 04, further comprising the step of: 

2 providing fuel when no electric power is available to operate said energy source at 

3 a no load condition to lower operating temperature in warm down state. 

1 106. The method claimed in claim 105, further comprising the step of: 

2 monitoring said system for faults. 

1 1 07. The method claimed in claim 106, further comprising the step of: 

2 disabling all outputs so that said system is placed in a safe configuration when 

3 faults that prohibit safe operation occur. 

1 1 08. A method of controlling the distribution of power in a system including a turbine, 

_ 2 among a plurality of energy components, using a computer including a digital signal 

yQ 3 processor comprising the steps of: 

[H 4 interfacing a plurality of power inverters between a DC bus and each of said 

^ 5 energy components; 

pi 6 controlling the way each of said energy components sinks or sources power; and 

f " 7 controlling the way said DC bus is regulated responsive to operation of each of 

D 8 said energy components, wherein said energy components include an energy source, a 

pj 9 load and a energy storage device. 

!Jf 1 1 09. The method claimed in claim 108, further comprising the step of: 
H ; 2 varying a speed command to regulate power of said system. 

1 110. The method claimed in claim 108, further comprising the step of: 

2 varying a fuel flow command to regulate speed of said turbine. 

1 111. The method claimed in claim 108, further comprising the step of: 

2 varying a fuel flow command to regulate exhaust gas temperature of said turbine. 

1 112. The method claimed in claim 108, wherein said power inverters include first and 

2 second power inverters under the control of first and second signal processors, 

3 respectively. 

1 113. The method claimed in claim 1 1 2, further comprising the step of: 

7 1 
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2 varying a current command associated with said first signal processor to regulate 

3 a speed of said turbine. 

1 114. The method claimed in claim 1 12, further comprising the step of: 

2 varying a current command associated with said second signal processor to 

3 regulate voltage of said DC bus. 

1 115. The method claimed in claim 1 12, further comprising the step of: 

2 varying a current command associated with said first signal processor to regulate 

3 voltage of said DC bus. 

1 116. The method claimed in claim 1 12, further comprising the step of: 

2 providing power from said DC bus in accordance with said second signal 

3 processor to provide a constant AC voltage output. 

1 117. The method claimed in claim 1 08, further comprising the step of: 

2 providing power bi-directionally from said energy storage device to regulate 

3 voltage of said DC bus. 

1 118. The method claimed in claim 108, further comprising the step of: 

2 providing power from said DC bus in accordance with said second signal 

3 processor to provide a constant AC current output. 

1 119. The method claimed in claim 108, further comprising the step of: 

2 varying an AC current command to said second signal processor to regulate a 

3 constant turbine EGT. 

1 120. The method claimed in claim 1 08, further comprising the step of: 

2 providing power bi-directionally from said energy storage device to regulate a device 

3 state of charge. 

1 121. A method of controlling the distribution of power in a system including a turbine, 

2 among a plurality of energy components, using a computer including a digital signal 

3 processor comprising the steps of: 

4 interfacing a plurality of power inverters between a DC bus and each of said 

5 energy components; 

6 controlling the way each of said energy components sinks or sources power; and 
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controlling the way said DC bus is regulated responsive to operation of each of 
said energy components, wherein said energy components include an energy source, a 
load and a energy storage device. 

1 22. The method claimed in claim 121, further comprising the step of: 
varying a speed command to regulate power of said system. 

123. The method claimed in claim 122, further comprising the step of: 
varying a fuel flow command to regulate speed of said turbine. 

124. The method claimed in claim 123, further comprising the step of: 

varying a fuel flow command to regulate exhaust gas temperature of said turbine. 

125. The method claimed in claim 124, wherein said power inverters include first and 
second power inverters under the control of first and second signal processors, 
respectively. 

126. The method claimed in claim 125, further comprising the step of: 

varying a current command associated with said first signal processor to regulate 
a speed of said turbine. 

127. The method claimed in claim 126, further comprising the step of: 
varying a current command associated with said second signal processor to 

regulate voltage of said DC bus. 

128. The method claimed in claim 127, further comprising the step of: 

varying a current command associated with said first signal processor to regulate 
voltage of said DC bus. 

129. The method claimed in claim 128, further comprising the step of: 
providing power from said DC bus in accordance with said second signal 

processor to provide a constant AC voltage output. 

1 30. The method claimed in claim 1 29, further comprising the step of: 
providing power bi-directionally from said energy storage device to regulate 

voltage of said DC bus. 

131. The method claimed in claim 130, further comprising the step of: 
providing power from said DC bus in accordance with said second signal 

processor to provide a constant AC current output. 

132. The method claimed in claim 131, further comprising the step of: 



96472.01 03 




153501-0053 



2 



varying an AC current command to said second signal processor to regulate a 



3 constant turbine EGT. 

1 133. The method claimed in claim 132, further comprising the step of: 

2 providing power bi-directionally from said energy storage device to regulate a device 

3 state of charge. 

1 1 34. A power controller for distributing power among a plurality of energy 

2 components, comprising: 

3 a DC bus; 

4 a plurality of power converters, each of which is connected between one of said 

5 energy components and said DC bus and is responsive to said power controller, wherein 

6 said power controller provides a distributed generation power system by controlling the 

7 way each energy component sinks or sources power and said DC bus is regulated; and 

8 means for detecting transients associated with one of said energy components. 

1 135. The power controller claimed in claim 134, further comprising: 

2 means for suspending power transfer between one of said energy components and 

3 one of said power converters. 

1 136. The power controller claimed in claim 135, further comprising: 

2 means for resuming power transfer between one of said energy components and 

3 one of said power converters once current in said one of said energy components has 

4 decayed to near zero, 

1 137. The power controller claimed in claim 136, further comprising: 

2 means for dissipating via a resistive load said power fed into said DC bus by said 

3 other one of said energy components. 

1 138. The power controller claimed in claim 1 34, further comprising: 

2 means for estimating phase voltage magnitudes and grid phase angle in a 

3 feedback process, in conjunction with measurements of actual phase voltages to improve 

4 estimated peak voltage magnitudes; 

5 means for estimating an instantaneous angle of each phase of an utility grid based 

6 on said estimated peak voltage magnitudes and measured phase voltages; 

7 means for utilizing most accurate angle estimate to calculate an estimate of an 

8 instantaneous phase angle of said grid; 



96472 01 03 




9 
10 
11 
12 
13 

1 

2 
3 
4 
5 
6 
7 
8 
1 

2 
3 
1 
2 
3 
4 
1 
2 
3 
1 

2 
3 
4 
5 
6 
7 
8 



153501-0053 

means for differentiating and filtering to form an estimate of grid frequency; 
means for integrating said grid frequency to produce an estimated grid phase 
angle; and 

means for correcting said estimated grid phase angle to converge in phase with an 
estimate of an instantaneous phase angle of said grid. 

139. A method for controlling the distribution of power among a plurality of energy 
components, comprising the steps of: 

connecting a power converter between a DC bus and each of said energy 
components; 

controlling the way each of said energy components sinks or sources power; 
controlling the way said DC bus is regulated responsive to operation of each of 
said energy components; and 

detecting transients associated with one of said energy components. 

140. The method claimed in claim 139, further comprising the step of: 
suspending power transfer between one of said energy components and one of 

said power converters. 

141. The method claimed in claim 140, further comprising the step of: 

resuming power transfer between one of said energy components and one of said 
power converters once current in said one of said energy components has decayed to near 
zero. 

142. The method claimed in claim 141 , further comprising the step of: 
dissipating via a resistive load said power fed into said DC bus by said other one 

of said energy components. 

143. The method claimed in claim 139, further comprising the steps of: 
estimating phase voltage magnitudes and grid phase angle in a feedback process, 

in conjunction with measurements of actual phase voltages to improve estimated peak 
voltage magnitudes; 

estimating an instantaneous angle of each phase of an utility grid based on said 
estimated peak voltage magnitudes and measured phase voltages; 

utilizing most accurate angle estimate to calculate an estimate of an instantaneous 
phase angle of said grid; 
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9 differentiating and filtering to form an estimate of grid frequency; 

1 0 means for integrating said grid frequency to produce an estimated grid phase 

1 1 angle; and 

12 correcting said estimated grid phase angle to converge in phase with an estimate 

13 of an instantaneous phase angle of said grid. 
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ANALYSIS OF EMISSION CONTROL BY A CATALYTIC RECUPERATOR AND BY 

UPSTREAM MONOLITHS 

INTRODUCTION 



catalytic monoliths and catalytically active exhaust-side recuperator surfaces to remove 
hydrocarbon and carbon monoxide emissions from the Capstone tobogenerators-^BB 



calculations, results, and conclusions. Also attached are the MathSoft Mathcad™ 6+ documents 
comprising the analyses and an Excel 5.0 spreadsheet of the findings. 



TRANSPORT AND KINETIC ANALYSIS OF THE RECUPERATOR AND 
MONOLITHS 

The analysis of the conversion of heptene (C7H14), methane (CH4), and carbon monoxide 
(CO), were completed for several configurations of a recuperator and honeycomb catalytic 
converters. The method used was a MATHCAD program, and the procedures were similar to the 
document sent to Capstonej(^^(p(CAPSTONE.MCD document). 

The analyses consist of the following features (see printouts fer-CPSTONE5rMGD-aad 
ePSTONE6.MCD>: 

Gas Properties 

• The properties of the gases used in the combustion mixture are described on 
pages 1-8. 

• The gas compositions both for complete combustion air with fuel equivalence 
ratio 0.1 in moist (2 vol% water vapor) air, as specified in the problem 
statement, and also for an adiabatic temperature rise from 1 000° to 1 5 1 1 JF 
using heptene are defined on page 2. The difference in properties between 
these two compositions is insignificant. 

• Essential transport properties, such as the diffusivities of the added gases are 
determined on pages 3-4. 

• Mixture properties such as viscosity, thermal conductivity, heat capacity., and 
various dimensionless transport numbers such as the Prandtl and Schmidt 
numbers are defined and calculated on pages 5-8. 





iing post turbine 
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Recuperator and Converter Dimensions and Operating Conditions 

• The dimensions of the honeycomb sizes from 100 to 800 cells/in. 2 are defined 
on page 9 along with the size and dimensions of the recuperator and the 
upstream monolith converters, 

• The five operating combustion conditions representing the different power 
level cases for the turbogenerator are specified on pages 9-10. 

• The pressure drop for the 600 or 800 cells/in 2 recuperator is determined on 
page 1 1, and the analysis of heat transfer within the recuperator is described on 
page 12. 

Recuperator Performance 

• Mass transfer rates to the wall of the recuperator (assuming very fast 
conversion kinetics at the wall) are described on pages 13-14 for all three 
added components (C7H14, CH4, and CO). These calculations used average 
temperatures to determine an integrated formula for fuel conversion. 

• An alternative form of conversion was developed by looking at local Nusselt 
number (modified for the square channel configuration), TTiis analysis is 
described on pages 15 and 16. The two methods show good agreement for 
C7H1 conversion. 

• Catalytic properties and functions are defined on pages 17-18, and the specific 
rate expression for oxidation of dilute CH4 by supported palladium oxide 
(PdO) is described on page 19. 

• Slow (in CPSTONE5.MCD and related to the rate of CH4 oxidation) and 
moderate (in CPSTONE6.MCD and related to the rate of C3H6 oxidation) rate 
constants for C7H14 oxidation are described on p. 19. 

• The performance of the recuperator for conversion with temperature-dependent 
mass transfer properties and specific catalytic rates for C7H14 and CH4 are 
described on pages 20-2 1 . 

Performance of Upstream Monoliths 

• The operating conditions for 200, 400, and 600 cells/in 2 honeycombs is 
described on pages 22-25. 

• The pressure drop for three honeycombs (200, 400, and 600 cells/in. 2 ) is 
described on pages 22-23, and matched to the maximum specified conditions. 
From of these calculations, two honeycomb converter configurations, 400 
cells/in. 2 by 4 in. length and 600 cells/in 2 by 3.5 in. length, were chosen for 
mass transfer calculations. 

• The effluent concentrations of added components are calculated for the 
constant temperature monoliths (400 cells/in. 2 in CPSTONE5.MCD and 600 
cells/in. 2 in CPSTONE6.MCD) for the five operating conditions (cases) on 
pages 24-25, assuming very fast catalytic rates. 

• Finally, on pages 25-27, the conversion of C7H14 and CH4 using finite catalytic 
rates are described for the honeycomb converters (400 cells/in/ in 
CPSTONE5.MCD and 600 cells/in 2 in CPSTONE6.MCD). 



RESULTS OF THE ANALYSIS 



The principal results of the analysis of post turbine fuel component conversion are shown 
in the attached spreadsheet tables (CAPSTN5.XLS) and Mathcad documents (CPSTONE5.MCD 
and CPSTONE6.MCD). The main results of the analyses are summarized as follows: 

Catalytic Recuperator 

- ^TheAP/P calculations seem to ve ry closely fit your measured values. 
flHHHHHBpBHHP We c h° se as the best calculation of pressure 
drop an 800 ceils/in^ x 7-in. with 10% additional added length to describe the 
chevron geometry of the recuperator foils. 

• The estimated averaged wall temperatures, or ski n temperatures , within the 
recuperator seemed a little low (< SOT) based onSK^gg^gKKgtB 

discussions of recuperator thermal efficiency. As a consequence, the cell size 
was increased to 600 cells/in 2 keeping the length the same (7 in. x 1.1) for a 
conservative calculation of mass transfer rates. 

• For the mass transfer limited cases (assuming fast kinetics), virtually complete 
conversion CH4 and CO is predicted (see Table of Results 1). Except at the 
highest power level, conversion of CH4 and CO was essentially complete 
within the first inch or two of entry into the recuperator. 

• For the transport-limited cases with C7H14, maximum conversions are low 
only for Case 5 (full power). Even in this case, 50 ppm would be converted to 
< 1 ppm (0.2 ppm) while concentrations entering the recuperator above 250 
ppm would exit at levels above 1 ppm 

(500 ppm 2.1 ppm). The situation improves for the other cases (e.g., Case 4 
gives 500 ppm 0.02 ppms while 5000 ppm => 0.17 ppm) . 

• Conversions would be much higher for transport limited conversion in a the 
800 cells/in. 2 recuperator, as calculations with this cell size were found 
reproduce the observed pressure drop. In this case, the logarithm of the inlet 
concentration divided by the outlet concentration for all added components 
would be -35 % greater with the smaller cell size. 

• The step-wise calculation which includes the temperature on flow velocity and 
transport properties showed reasonably good agreement with the integrated 
calculation -using average properties. This method gave lower overall 
conversion than the averaged property formulas, but the results were close. 
Perhaps the higher property values at the inlet, particularly the greater velocity, 
more than offset the lower transport values at the outlet side. 

• CH4 conversion is catalytically rate limited under all conditions (see Table of 
Results 1). Conversions were calculated to be very low, with a maximum 
conversion of -66% under Case 1 and -32% conversion under Case 5. These 
conversions would increase somewhat, perhaps by a factor of 2, with the use of 
a thicker coating (e.g., 25 jam vs. 10 ^m). 

• CO conversion is likely to exceed 99% in all operating conditions for the 
recuperator. This is because the kinetics of CO conversion are probably much 
faster than CH4 and may well be close to the transport-limited case, at least 
through the (hotter) first several inches of the recuperator where most of the 
CO would be consumed. 



• The catalytic recuperator, C7H14 rate determines its conversion in the recuperator 
conditions (see Table of Results 1). For the slower rates (derived from the rates 
for CH 4 oxidation) using the 600 cells/in. 2 value for recuperator cell size, we 
find a range from 99.4% conversion under Case 1 to about 85,6% under Case 5. 
If twice the rate of combustion of propane is extrapolated to the higher 
temperatures of the recuperator, much higher C7H14 conversions result 99 98% 
under Case 1 , and 96.4% under Case 5. 

Catalytic Monoliths 

• The pressure drop reaches maximum specifications for the 400 cells/in 2 
honeycomb at 4-in. length. Maximum pressure drop would be reached in a 
shorter (-3.5 in. length) 600 cells/in. 2 monolith. Both configurations were 
chosen for analysis (CPSTONE5.MCI> for 400 cells/in 2 and 
CPSTONE6.MCD for 600 cells/in 2 ). 

• Transport-limited conversion shows excellent performance for methane and 
CO in all cases, including the full-power case, for both honeycomb 
configurations (see Table of Results 2). 

• As for the recuperator, methane conversion is again reaction rate limited (see 
Table of Results 3), with conversions range from 93.8% under Case 1 with the 
600 cells/in. 2 monolith to only 3 1 % conversion under Case 5 for the 400 
cells/in. 2 monolith. 

• C7H14 conversions under transport-limited conditions depends on the specific 
configuration (see Table of Results 2). For the 600 cells/in 2 monolith 
conversions are very high for Case 1, but fall to about 92% for Case 5. With 
the 400 cells/in 2 configuration, these fast kinetic transport-limited conversions 
fall to 87% conversion under Case 5. 

• The reaction rate of C7H14 strongly affects its conversion in either honeycomb 
m ?w- ? ^ See Table of ResuIts 3 )- For to extrapolated propane kinetics at 600 
cells/in/, 98% conversion is expected under Case 5, whereas the 400 cells/in 2 
using the more conservative extrapolation using PdO methane kinetics, 86% 
conversion is expected under Case 5. 

Effect of Washcoat Thermal Conductivity in a Recuperator 

• The application of a thin washcoat of porous alumina gives a very slight 
increase in the overall heat transfer coefficient^MM|MM| 

CPSNALT2.MCD). The thermal conductivityortn^^shcS, while 
significantly less than that of stainless 321, is greater than that of the gas phase 
boundary. Because the washcoat displaces the gas phase (creating a narrower 
channel) it causes a slight increase in overall heat transfer coefficient. 



OVERALL CONCLUSIONS AND RECOMMENDATIONS 



The conclusions and recommendations drawn from the results of the analysis is listed 

below: 

1) The actual C7H14 outlet concentration from the turbine are needed to assess the 
effluent emissions in the high-power case (Case 5) for either control strategy. 
Provided that fast catalytic rates can be provided in washcoats on either the 
recuperator or upstream monoliths, maximum conversions for C7H14 under full 
power range from 99,6% (for an effective cell size of 600 cells/in 2 ) to > 99 9% 
(if the effective cell size is 800 cells/in 2 ). 

2) The two methods (averaged properties and integrated formulas vs. local 
temperature dependent properties and finite difference integration) show good 
agreement for recuperator performance with infinitely fast catalytic reactions at 
the wall. 

3) Clearly, it will be necessary to apply and measure specific catalytic rates for 
C7H14 or gasoline to be sure that these estimated conversions could be 
achieved under actual operating conditions for either the catalytic recuperator 
or the upstream honeycomb monoliths. 

4) Thin wash coats (< 20 fun) should have minimal effect on the transport of heat 
through the recuperator given the much greater resistance of the gas phase 
boundary layers. 

5) For the current operating conditions, use of a standard 400 cells/in 2 x 

4-in length catalytic monolith upstream of the recuperator is recommended. 

6) If the recuperator temperature can be raised (~ 200°F) high conversions of CO, 
C?Hi4> and even CH4 can be expected under all flow conditions. In this case it 
is possible that only the first 3 to 4 inches of the recuperator be wash coated 
with catalyst because the downstream wall temperatures are too cool to be 
effective. 
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*****************************************************^^ 

Gas Properties and Constants 



cm 0.01 *m 
N : kg*m-sec 



R -8.3143- 



J-kg*m 2 *sec 2 
Pa N*m 2 



W-J-sec" 1 
atm-101300-Pa 



Pref = 1 -atm 
Mm- 0.00000 J *m 



Rc 8.3143-J-mole '-K' 1 



»»***»»*»****Ht»***;*«*»************»******** ******** ***** **************** ********************************** 



n Gas 



MolWt 



a(U) 



*LJ> 



0 CH4 MW o " l6 - 044 *g m °o r 3.785.10" ,0 -m e Q = 148.6K 

1 02 MW. = 3I.999-gm o = 3.464- 10" I0 -m e, = 106>7-K 



2 N2 

3C02 

4H20 



MW 2 = 28.018-gm o 2 = 3.798- 10" -m e 2 =71.4-K 
MW 3 = 44.010-gm o 3 =3.941-10" 10 -m 6 3 = 195.2-K 
MW 4 = 18.015-gm o 4 = 2.64M0" 10 <m £, = 809.1-K 



5C7H14 

methyf-cyclohexane 



Tc7 
1.2593 



Tc7 = 572.2-K Tb7 = 374.1 -K Pc7 =3.47-10 -Pa 

l 

a 5 = Vc7 3 *10" l0 -m 



MW 5 = 98,189-gm o $ -7.1661*10 *m s 5 = 454.3794* K 



6CO MW fi - 28.01 -gm ° 6 = 3.69*10" 10 -m e 6 = 91.7-K 

iHhrt***************^***lHMNrt***^^ 



'6 '-—6'" "6 

*************^**************^^ 

Define constants for transport integrals; 



<*a 0 


= 1.16145 


°l* 0 


= 0.14874 


aDa 0 


= 1.06036 


cDb 0 


= 0.1561 




= 0.52487 




= 0.7732 


oDa 


= 0.193 


aDbj 


= 0.47635 




= 2.16178 




"2.43787 


<?Da 2 


- 1.03587 


aDb 2 


= 1.52996 










GDa 3 


- 1.76474 


oDb 3 


= 3.89411 



Define transport integrals: 

Q«T) =<7iia 0 *expj-,o^b 0 -ln(T). 

QD(T) : oDa 0 -expi oDb 0 -In(T) 



2 

ip=l 

3 

V 1 



o^a -exp; opb -T 



aDa -exp oDb T 



0^(3) = 1.0394 



QD(3) =0.95002 



**************************************************************** 



Vc7 =368 
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**************** *************************************************^^ 



Set Inlet Conditions: 



yeqc7 " f 



0.205 



14 



yeqc7 =0.001464 



P 7 1 -atm 
yadc7 0.0021859 

, . yadc7 
<feq -<p-~ " 
yeqc7 

$eq =0.149281 



This valu e for yadc7 was determined via enthalpy 
functionsfHBHHHHHBHftf° r temperature 
rise from 1000F to 151 1F. 

The value for yeq was determine for a given equivalence 
ratio of 0.1 in moist air 



Define gas composition with C7H14 fuel using moist (2%) air: 



Jy7(y0,x) 



0 

{I - y0) 0.205- 10.5-y0*x 
(l-y0)-0.775 
7-(y0.x) 
: (1 -y0)-0.02+7«y0-x 
! yO-(l-x) 
0 



1 - 2.5-yO-x 



| yeq z y7(yeqc7 > 1) 
Jyad -y7(yadc7,l) 



yeq 



[0 

| 0.188634 
; 0.771043 
= j 0.010213 
; 0.03011 

I 0 
[0 



yad : 



! 0 
0.180613 
0.769103 
0.015218 
0.035066 
0 
0 



******** ********************** * * **** * * 



*********************** 



♦a******************************************** 



Define functions for evaluation of mixture properties (ideal gas) for complete fuel conversion and a given temperature. 
Density and average molecular weight: 



pmix(T) 



6 

V 
kk = 0 



yeq kk MW kk !■— 



MWmix 



kk-0 



yeq kk -MW 



kk 



pmix(300-K) = 1.162787'kg'm 3 



MWmix -28.631 "gm 



CAPSTNE5.MCD 



***************** 



Define diffusion coefficients for specific binary diffusion pairs: 



DifJ(n,m,T) r 1.8829- 10 



!;mw mw 

27 L \ n n 

,2 



1 -T 3 



0.5 



.10 8 -gm 05 cm 4 sec-'K' 15 



0_ +■ 0_ i i T 

v i_ 

n m/ 



P 

Pref 



Diff(0,2,300-K) -0.2225'sec 1 *cm 2 



Diff(2,6,300-K) -0.207'sec 1 'cm 2 



Define functions for use with various fuel-air mixtures: 

^Dch4(T) 



] Dch4(298-K) =2.1912*10 5 Tn 2 -sec 1 



jgas= 1 



Diff(OJgas,T) 



^Dc7hl4(T) 



^ ? Diff(5,jgas, 



T) 



Dc7hl4(298*K) =5.8372*10 6 "mSec 1 



^Dco(T) 



4 

jgas= 1 



yeq 



Jgas_ 



Diff(6Jgas,T) 



Dco(298-K) =2.0451*10 5 'mNec 1 



**************************************** 
Display gas mixture properties 

idsp 0..9 



Tin -1(1000 - 32)-- + 273.15 |-K Tin =810.9278*K 

L 9 J 



Tdsp =(750-50-idsp).K Tout =j(15ll 32)-| + 273.15 j-K Tout = 1094.8 1 67* K 



4 



2 5*10 4 



2*10 



'4 _ 



Dch4(Xdsp ida| J l 5 ., 0 -4 
^7hl4(Tdsp idsp i 
Dco ( Tdsp i£i j ,. l0 -4 



5*10 5 




1200 



yTdsp 



idsp_ 





7.5-10 2 -K 




8-10 2 -K 


r j 


8.5-10 2 -K 




9-10 2 -K 




9.5-10 2 -K 




M0 3 -K 




1.05- 10 3 K 




l.MOMt 




I.1510 3 *K 




l^lO^K 



1.082- lO^V-sec 1 


3.1097- 10" 5 -m 2 -sec -1 


12064- 10~ 4 m 2 -sec 1 


3.4776-10" 5 -m 2 -sec _1 


1.3359- 10" 4 m 2 sec _1 


3.8613-10" 5 m 2 -sec" 1 


1.4705-10" 4 *m 2 *sec" 1 


4.2603- 10* 5 -m 2 -sec -1 


1.6M0~ 4 -m 2 -sec~ l 


4.6743- lO'W-sec 1 


1.7543-10~ 4 -m 2 -sec _1 


5.1031- 10" 5 -m 2 -sec' ! 


1.9033- loW-sec 1 


5.5462- lO'W-sec 1 


2.057- 10" 4 m 2 sec" 1 


6.0035-10 "^m^sec" 1 


2.2152- I0" 4 *m 2 -sec" 1 


6.4746- 10" 5 -m 2 *sec" 1 


2.378- 10" 4 -m 2 -sec 1 


6.9593- loW-sec -1 



Dco(Tdsp idspZ 



9.8836-10" 5 - 


m 2 -sec _l 


1.10M0" 4 - 


m 2 -sec" 1 


1.2184-10" 4 * 


m 2 -sec _1 


1.3402- 10" 4 - 


m 2 -sec~ 


1.4665- 10" 4 


m 2 -sec" 1 


1.5972- 10" 4 


m 2 -sec" 1 


1.7322- 10' 4 


m 2 -sec~ 


1.8714- 10" 4 


m 2 sec _1 


2.0148- 10" 4 


m 2 -sec" 1 


2.1624 10" 4 


•m 2 -sec _1 



Define interpolation formulas for Dch4, Dc7h14, and Deo vs. T(K) via Jandel Science's Tablecurve program 



< n t ,'T\ 46.41736-K; cm 2 
Dch4T(T) : = exp^- 10.623294 t 1. 6259688- In^-J )•— 



Dch4(1000-K) =1.75426*sec 1 'cm 2 
Dch4T( 1 000-K) -1.75432 'sec 1 "cm 2 



. | T\ 82.65375-K; cm' 
Dc7hl4T(T)- = exp» 11.822189 ■ 1.62601 381n[ 



\K/ 



sec 



Dc7hl4(1000-K) =0.51031'sec ' 'cm 2 
Dc7hl4T(1000-K) =0.5103Tsec 1 xm 2 



24.53899-- )-cm 2 -sec 1 



Dco( 1 000* K) -1.59722 'sec 1 'cm 2 



DcoT(T) ^expf 10.831955 t 1 .6394196-lnj -- j 



DcoT(lOOO-K) = 1.59714'sec 1 *cm 2 
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Define viscosity functions for major constituents: 



KT): 



2.6693*10 



05 



t MW 2 -Tj 



0.5 



■ MW 3 -Tj a5 
1.116-iMW -T) 



26 iMW 0 -T. 05 10 5, cm , -0.5 
, .2 _ T sec \ K j 
£ o/ 



2.6693-10 



2.6693-10 



_26 /10jcm\ /gm 



sec 



0.5 



Ik 



sec 



I (3.115'I0" s -cm) 2 ^fi{-^— s 
1 n 514-K 



sec / 



0.5 



1.10281M0 
2.0603' 10" 5 
1.76981*10" 
1.51851-10" 
1.06784-10" 



kg m *sec 



Hmix(T) := 



4 

V 



yeq n -KT) n 



1 + 



;n(T) \°- s /mw ^ 2S 



[ . MW \ 

!s-ii + — «~\ 



0.5 



yeq 



0 

0.1886 

0.771 

0.0102 

0.0301 

0 

0 




700 



1200 



Tdsp. 
|750-K 



1800-K 



850 K 



900-K 



950-K 



1000-K 
1050* 
[ n00-K 
[ 1150-K 
1200* 



fimix^Tdsp 


idsp/' 




3.3875-10" 


5 <kgm" 


^sec" 1 


3.5348 10 


5 -kgm- 


-sec 


3.6787-10" 


5 -kg*m 


'•sec 1 


3.8194-10" 


5 -kg-m" 




3,9573-10" 


5 -kg*m 


1 -sec" 


4.0926-10 




1 -sec" 1 


4.2254-10 


_^kgm 


-sec 


4.356-10" 


5 ;kg*m 




4.4844-10 


5 -kg-m 


'•sec 1 


4.6108- 10 


5 -kg.m 


'-sec 1 
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MxT(T) ^expj 0.6984944 ^ 0.640549-lnj - ) f — ~- 5 - 5 — )-10' 6 -kg-m' '-sec 1 

\K;' T I 



Mx( 1000-K) =4.0926*10 5 'kg*m ! 'sec 1 



HmixT(l000-K) =4.09256*10 5 'kg"m Sec 1 



Define thermal conductivity functions for light molecules 

igas -0,1.. 4 



aX:= 



1.869 




" 0.08727 " 




1.179-10 4 




(3.61410 s ) 


- .3273 




0.09966 




-(3.743-10 5 ) 




9.732- 10" 9 


03919 


bX = 


0.09816 




-( 5.067- 10 -5 ) 


dX - 


1.504- 10"* 


-7.215 




0.08015 




5.477- 10" 6 




(1.053-10 8 ) 


7341 . 




.-0.01013 J 


1.801- 10' 4 




t (9.1-10 8 ) 



CH4 

02 
N2 

C02 

H20 



T / T < 2 I T \ 3 

aX. + bX..i + ci..-1-i- +- dX.-i i 
1 1 K ' \K/ ' IK 



Xmix(T) := Y, 



i = 1 



0.00001- 



W 
cm-K 



yeq-X(i,T) 




X(igas,1000K) 


0.16716-kgmsec 


3.K-I 


0.07163 -kg- m- sec 




0.06292-kgmsec" 


3.K-1 


0.06788*kgmsec 


3 «K l 


0.08631 -kg- msec 


3. K -i 



W/nrWK 



Xmix(1000-K) =0.000653 



W 
cm-K 



0.075 



0.065 - 




0 055 



1200 



Tds PKf<n 


XmixiTdsp idsp ; 


7.5-10 2 -K 




0.05323 kg m-sec 3 K 1 


8-10 2 -K 




0.055775 kg-m*sec 3 K ' 


8.510 2 -K 




0.058242 kg-m-sec~ 3 K 1 


9-10 2 -K 




0.06064 1 -kg-m-sec 3 K 1 


9.5-10 2 K 




0.06298 kg*m*sec 3 K 1 


M0 3 K 




OX^eS-kg-m-sec" 3 -^ 1 


L05-10 3 -K 




0.0675 l2kgm*sec" 3 *K 1 


I.IICPK 




0.06972 lkg-m-sec" 3 K~ ! 


1.15-ltf-K 1 




M? 1904 kg m-sec 3 -K 1 


L210 3 -Kj 

__ i 


0.074069-kg-m-sec 3 K 1 
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XmixT{T) -expi 7.073383 r 0.6397005-!n; T ! - - 65 "— — i-W-m l -K~ l 

K. T i 



lm\x( 1000-K) = 0.0653*m 1 *K~ ! -W 



A.mixT{1000-K) =0.0659'm n -K 1 *W 



Define heat capacities: 



Cpcoef 



19.25 .05213 0.0000 1197 0.00000001132 

28.11 -0.000368 0.00001746 0.00000001065 

31.15 -0.01357 0,0000268 -0.00000001168 

19.8 0.07344 -0,00005602 0.00000001715 

32.24 0.001924 0.00001055 -0,000000003596 

76. 1 9 0.7867 - 0.0004204 0,000000075 1 6 



Cp(igasJ) = 



3 

s 

11 = 0 



Cpcoef 



igas t j 



Aki k 



jjj = 0..3 



J/mol/K 



^Cp(0 > 300K)=35.6607*K 1 'J 
4 

Cpmix(T) - Cp(jk f T).yeq. k 
jk-1 



4 

Cpcfinix. .= £ Cpcoef jk . yeq jk ^ 



jk-1 



Cp(5,300-K) = 124.0133'K ' 
Cpmix(998-K) -33.5 144 ^K' ! 



Cpcfinix = 



30.493459 
-0.009725 
2.370305' 10~ 



1.094786-10 



*kg"m 'see 



CpmixT(T) - ]T Cpcfinix^ (|j JJ 

ii-o 



Cpmix( 1000-K) =33.52413'kg'm 2 -sec 2 *K 1 
CpmixT( 1000-K) = 33.524l3'kg-mNec~ 2 'K 1 



igas 0,1. .6 P = 1.013* 10 5 'Pa 
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************************************* 
Define momentum and thermal diffusivities: 



a(igas,T) 



X(igas /T);R-T 
Cp(igas,T)-P 



ct(0,300-K) =0.2344'sec 1 *cm 2 



amix(T) 



,Xmix T(T)-R-T 
CpmixT(T)-P 



amix( 1000-K) = 1.6132*sec 1 "cm 2 



umix(T) 



umixT(T) 



- IMhxTTT) 
pmix(T) 

,pmixT(T) 
pmix(T) 



omix(lOOO'K) - 1.1732'sec 1 *cm 2 



pmix( 1000-K) =3.48836" 10 4 *gm-cm~ 3 



pmixT(T) ^pmix(T) 



Dch4T( 1000-K) = 1.7543' sec 1 *cm 2 



Dc7hl4T( 1000-K) =0.5103*sec 1 'cm 2 



yPDefine Prandl, Schmidt, and Lewis numbers 

: umix ( T ) Pr( 1 000-K) = 0.7272 



rfr(T) 



amix(T) 



ftclT(T) := ^J1 SclT(lOOOK) =0.6688 Sc7T(T) = J™*^ 



S>ccoT(T) - 



Dch4T(T) 
omix(T) 



DcoT(T) 



SccoT( 1000-K) =0.7346 



Sc7T( 1000-K) =2.299 



SclT(Tdsp idsp ) 

Sc7T ( Tds Pidsp) 
2.5 



0.9 



0.8 



SccoTiTdsp 



idsp/ 



0.7 



0 6. 



700 



800 



900 



1000 



™sp 1( 



H00 



1200 



These relations have weak 
temperature dependence and 
(presumably) very weak 
dependence on composition 



This ends the extended section that defines gas properties for the combustion mixtures. 
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************************************************************************************************************* 

Honeycombs: Represent as square channel metal alloy waft 0.003 mil, uncoated; 0.004 mil coated 



dwalluc - 0.003 -in 



dwalluc =7.62* 10 3 *m 



dwallc 0.004 In 



dwailc = i. 016- 10 *m 



Cells 1 -100-in^ dwalll = dwallc 

dcellt -(Cellsl) 05 dwalll 

dcelil = 0.2438 *cm ksl ^0.002-cm 



Cells2 : 200-ui dwa!I2 = dwallc 

dcel!2 (Celis2) 05 dwaII2 
dceU2 =0.1694 'cm ks2 -0.002-cm 



ecelll - dcelil Cells I 



ecelll -0.9216 



ecell2 = dcell2 -Cells2 



ece!12- 0.89006 



Ceiis4 = 400-in 



dwa!14 -dwallc 



Cells6 600- in 



dwall6 -dwalluc 



dce!14 : = (Cells4)" 05 - dwa!14 dce!14 -0.1 168 'cm dceH6 ^(Cells6)"°- 5 - dwal!6 dce!16 =0,0961 'cm 

dce!14 =0.1168*cm ks4 -0.001-cm dce!16 =0.096 rem ks6 -0-crn 



ece!14 ~- dcelI4 -Cells4 ecell4 - 0,8464 



ecell6 = dcel!6 -Cells6 £ceI16 = 0.85843 



Cells8 : = 800-in' : 



dwa!18 -dwalluc 



\ dcellS ={Cells8)" 0 - 5 - dwallS dceH8 =0.0822 -cm 
1 dcell8 =0.0822 *cm ks8 = 0-cm 
^scel!8 = dceI18 2 Cells8 ece!18 = 0.83749 

~\ ***** *** *»*»***»»^*****»»*»»*fr**»** ****** * ******************'■ 

Define recuperator and monolith configurations: 



♦it************************************ ********** 



5 dreco ~17>in 
: dreci =6 -in 
6d ^ 1.375 in 
dhco = dreco - Sd 
dhci - dreci +■ $d 



She 



dhco* ~ dhci 



dreco =43.18 "cm Lrec =7-in 

dreci = 15.24 'cm Lrec = 1778 'cm 
Add space for mounting and insulation? 

dhco = 39.6875 'cm Lhc ^4-in 

dhci - 1 8.7325 'cm Lhc = 10.16 'cm 

She =0.096 I'm 2 



Srecup -0.0958*m 2 



***************************************** 
Define initial combustion conditions 

jcase ~0 t K.4 



^ -0.9964 
She 



0.2 1350 340 i 

0.26 1350 390 I 

case , 0.314 1350 440 i 

0.38 1350 490 : 

; 0.72 1190 675 I 



Mass flow (lb/sec) 
Recup T inlet (F) 
Recup T outlet (F) 



TKF(T) =[- - 273.15 -32 
\K 5 



TFK(t) (t 3 ^ } ' 5 * K . 273.15-K 
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Prec 1-atm 



Phc - 1-atm 



10 



mf . :: case. *0.454*kg-sec 



Trecin J _. -TFK( case j 

jcase-M * jcase, I, ' 



Trecout ^ . - TFK ( case. ^ , i 

jcase +- 1 * jcase , 2 / 



mf= 



0 

0.0908 

0.118 

0.1426 

0.1725 

03269 



kg* sec 



mf 



Trecin = 



Urec. , . 
jcase -t-l 



jcase + 1 



Trecin. 



jcase 



Srecup 



Spnly half of the recuperator cross section 
Is available for exhuast gas flow 



o 

1005,37 
1005,37 
100537 
100537 
916.48 



Trecout : 



0 

444.26 
472.04 
499.82 
527.59 
630.37 



Urec - 



! 0 

' 5.4633 
7.1023 
8.5774 
10.3803 

L 17.929 J 



m'sec 



"I 



Grec 



jcase 



mf , . 

_ jcase + 1 

+ 1 | Srecup | 



jcase +• 1 



mf , . 

jcase -h I 



pmix(Trecin case + 1 )-Shc-£cell2 



scel!2 



ecell2 



=Uhc2 - 



0 

3.058 
3.9753 
4.801 
5.8101 
l 10.0353 



nvsec 



Uhc4 = 



0 

3.2157 
4.1804 
5.0487 
6.1098 
10.553 



m'sec 



Tavg 



Trecin +■ Trecout 



I 



Tavg- 



724.8167 
738.7056 
752.5944 
766.4833 
773.4278 



Define average temperature etc/ for 
estimated AP and AT in the recuperator 



pmixrecm M Wmjx^Prec pmixrec^Tavg. ■ -0.4813'kg'm" 3 

R-T 
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Evaluate Reynolds numbers for use in calculating the momentum, heat, and mass transfer functions: 



Rerec6(G,T) 



. G ' d ? c !!l 
ymixT(T) 



Rerec6^Grec ( jTavgj | =54.9906 



Rerec8(G,T) 



G*dcelI8 
^imtxT(T) 



RerecSfGrec^Tavg,] =47.0389 



For laminar flow with entrance effect (adapted from Kays & Crawford, "Convective Heat and Mass Transfer"): 



cfricrec8(G,T,z) - — 1 1 + 0.046263- fRerec8(G,T)-^^ ) 

Rerec8(G,T) [ \ z ! 



10.45363 



14.227 



1 + 0.046263- Rerec6(G,T) 



Average coefficient of friction 
dcell6 \ 0 45363 for a smooth square channel 



Rerec6(G,T) 

^ Calculate monolith pressure drop in recuperator with entrance effect: 



(Urec. , A 

\ lease -f 1 y 



= &Prec8. ^. -4-cfricrec8iGrec, ^Tavg. e .,LrecV- 

-w jcase + l ^ jcase-j-t* °jcase + l' / 



pmixrec^Trecin^^J + pmixrec(Trecout Jcase + M 



I Treci 
Tree 



(Urec , . ) 
\ j case -h 1 / 



jAPrec6. , -4cfricrec6|Grec. , ,,Tavg. ,,Lrec)* 



pmixreciTrecin.^^ { ) + pmixrecj Trecout jcase ^ , 



/Treo 



jcase + 1/ Trec 



APrec8 
Prec 



0 ! 
0.0097 
0,013 | 
| 0.0162 
j 0.0203 
j. 0.0393 , 



APrec6 



0 

0.0071 

_ ( 0.0095 
Prec 10.0119 
: 0.0149 
I 0.0289 



APobs P 



! 0 I 
j 0.011 
i 0,014 

' 0.018 

i 

: 0.022 
0.042 



Chose 800 cells per in2 
with 10% length factor 
in Lrec for chevrons 
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Averaged Nu for 800 cpi square channels (fully developed 
laminar flow) taken from Rosner: 



Nuhavgrec8 



jcase +■ I 



= 2.976- 1 f 



7.6-(Lrec-1.1) 



Rerec8(Grec A ,,Tavg. M , . j-PrJavg, ^ t )*dcell8 : 



Nuhavgrec8 = 



fihrecavg 



0 

2.976 
2.976 
2.976 
2.976 
2.9762 



Entrance effect for heat transfer within the recuperator is only 
significant for the first half inch - average Nusseft No. = 2.976 for all cases 



. knix(Ta vg jcase + J >Nuhavgrec8 jcase ^ y 



>jcase + 1 



dcel!8 



hrecavg : 



0 

188.0022 
190.6356 

193.2448 

i 

195.8307 j 
197.1273 j 



j"kg*sec~ 3 *K~ l 



W/m2/K 



JATskin 



jcase -f- 1 



MWmixhrecavg t -4-Lrec-Ll 



ATskin = 



0 

-7.5657 

"9.274 

-10.5339 

"11.9549 

-13.635 



ATskin*2 = 

5*K 



0 1 
"27.2364 ' 
-33.3864 
-37.9221 
-43.0378 
-49.0859 



Approximate average AT (F) 
across wall in recuperator 



AT seems low - more conservative calculation uses dceil for 600 cells/in2 - in this case AT maximum 
increases from 50F to about 70F 



5-K \dcell8. 



0 i 
-37.2231 j 
45.628 j 
51.8268 j 
58.8182 j 
67.0839 



, dceii6 / 
dcelI8 



= 1.3667 
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Use 600 cells/in2 in recuperator dimension for conservative calculation of CH4, C7H14, and CO conversions 



Find mass transfer rates via Nusselt formulas: 



Numavgrec6cl jcase ^ x - 2.976- 



1 +■ 



7.6-Lrec*U 



. 3 



Rerec6|Grec. ,,Tavg ,j-SclT.Tavg e ir dcell6 



Numavgrec6cl 



Numavgrecec?^^^ -2.976- 



H-i 



7.6-Lrec*l.l 



,8 

\ 3 



1* 



Re re c6(G reCjciiseM) T a vg jcase + 1 )-Sc7T(Tavg jcase+i; .dcell6/ 



Numavgrec6c7 



I 

2.976 

2.976 

2.9761 

1 2.9761 

L 2.9765 
0 

2.9765 
2.977 
2.9775 
2.9784 
2.9891 . 



TNumavgrec6co. , . -2.976- 

JCESC T" 1 



1 



7.6-Lrec-U 



Rerec6(Grec ^Twg^ + ^SccoT^avg.^ ^dceltf 



Numavgrec6co - 



0 

2.976 
2.976 
2.9761 
2.9761 
2.9766 J 



Epnce again the entrance effects within the recuperator are negligible 



f|The average (integral) Nusselt numbers at 1-atm show very little contribution of the thermal entry length, i.e. the flow 
^velocities and heat and mass transfer coefficients have values close to those in long channels. 



Now determine the maximum possible CH4, C7H14, and CO conversions in a 600 ceils per in2 recuperator 



ych40 =0.005 



yc7hl40 '0.0005 



ycoO =0.005 



ych4rec6 



: ych40*exp 



- 4-Lrec- l,l«Numavgrec6cl 



jcasc + 1 



jcascti ' MRerec6(Grec jcaseM ,Tavg jcase + J-SclT^avg.^ M j-dceI16 



Start with maximum 
concentrations of added 
pollutants 



4- Lrec- 1 . 1 'Numavgrec6c i jca$e + , : 0 

Rer^Grec caK ,Tavg jcase ^S^^g""^ y'*^ 16 1 ° 

654083 ych4rec6^;° 

50.9771 y ch40 ' 0 

"427552: 0 

35J764 5.5932M0 9 
19.0017 
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yc7h!4rec6 



: yc7h!40-expi 



- 4-Lrec- 1 .1 •Numavgrec6c7 



j case i I 



-- ! 



yc7hl4rec6 > 



3-5287*10 
2.1834-10 
2.2859M0 
1.6805' 10" 
2.0597M0" 



12 



10 



~9 



yc7h!4rec6 
yc7h!40 



0 

7.057-10" 
4.367-10" 
4.572*10" 
3.361-10 
4.119-10" 



~5 



4-Lrec- 1 , 1 -Numavgrec6c7. 



jcase +• t 



Rerec6('Grec , ,Tavg. ^ , )-Sc7T, Tavg. ^ . )*dcell6 
jcase +■ I ' & jcase +- 1 / \ °jcase +- 1 / 



18.7692 



14.6441 



12.2956 



10.3007 



-5.492 



ycorec6 



: ycoO-exp 



4*Lrec*l.l Numavgrec6co. 



jcase +- 1 



*« 1 1 Rerec6(Gre W , ,Tavg jcase J -SccoT^Tavg.^ ^ , i-dcelI6 



U^corec6 = 



0 
0 
0 
0 
0 

1.457-10" 



10 



ycorec6 
ycoO 



0 
0 
0 
0 

6.4676-10 
2.914" I(P 



15 



4 • Lrec- 1 . 1 •Numavgrec6co 



jcase j- 1 



59.7854 



46.5803 



39.0561 



- 32.672 



17.3512 



Ljfl/ith SOOOppm C7H14 in turbine exhaust, exit concentration reaches 20 ppm at the highest power level 
for conservative 600 cefis/in2 calculation. For 800 cells/in2, the exit concentration in this case would reach 
about 2 ppm. With only SOOppm C7H14 in the inlet, should get <2 ppm in the outlet of the recuperator 
assuming favorable catalytic rates at the wall. Also this calculation used averaged temperature and flow 
conditions - will compare below a 1 -dimensional (axial) calculation that determines the effect of decreasing 
temperature on the kinetics and mass transport. 



Maximum possible (transport-limited) CO and CH4 conversions in a 7 in. 
600-800 cells/in2 recuperator are >99.99%, while cyclohexane conversion is 
>99.6%. The final CO and C7H14 levels are < 1 ppb and .2 ppm for post 
turbine concentrations of 1000 and 50 ppm, respectively. 
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********************************* 



An alternative method of calculating performance uses local Nussett numbers for heat and mass transfer to determine the 
wall temperature and combustion rates at the wall. First we define the local Nu function then check conversion assuming 
fast kinetics and compare with the previous averaged (integrated) Nu formula. 



Local Nu(x) taken from Kays & Crawford for hydro & thermal laminar flow entrance 
with circular channel (modified for square channel by a factor 2.96/3.66: 



ncase -5 



Numono(z) : = (! ♦ exp(- 1.09689 - 0381402-In(z) - 8.287774-z° 5 ))-2.976 Numono( 0.1 152163) =3.112 



zzchar6c7 



Rerec6(Grec, ca ^ Lrec 



jcase +- 1 



zzchar6c7 



= 0.9444 



zzchar6cl 



Rerec6( Grec. case ^ , .Trecin^ + ,)-SclT (Trecin.^ + ,) -dcell6 Lrec 



jcase -h 1 



zzchar6cl 



Case specific 



3.2574 Case specific 



Fnz-400 
Hz r 0.. nz 



s zrec6. r iz -f dwal!6 
3 ,z nz 



Lrec = 17.78 *cm dcell6 = 0.0961 *cm 



zrec6. 



zrec6. 



zz6c7 



IZ zzchar6c7_ 



zz6cl 



zzchar6cl 



zzchar6c7 = 



0 » 

4.896 

6.3647 

7.6867 

9.3023 

18.8265 



'cm 



Numono(zz6c7^ 
Numono (zz6cl ] 




no 



zzchar6cl = 



0 

1.4244 
1.8518 
2.2364 
2.7064 
5.4583 



cm 



Case specific ncase = 5 
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kgpbI6c7(T,z) 



Dc7hl4T(T)Numonoj 



dcell6 



zzchar6c7 



Case specific ncase = 5 



kgpbl6c7; Trecin^, Lrecj = 13.6159'sec *cm 



Dch4T(T)*Numonoi 



zzchar6cl 



kgpbl6cl(T t 2) 



dcell6 



kgpb!6cl fTrecin ncase ,Lrecj =46.9581 'sec "cm 



Trecin - Trecout 

rj, - ^ ncase ncase - 

Trecz. Trecin 12 

iz ncase ^ 



Trecz = 630.3722' K 

nz 



Trecout = 630.3722' K 

ncase 



Determine maximum conversion of fuels with very fast catalytic combustion rates (zero fuel concentration at the wall) 



sv yzc7hl4 z yc7hl40 



yzch4 =ych40 



ncase = 5 



Trecin \ 
tcellO ncase 



Hw7hl4. -II- kgpbl6c7(Trecz.,zrec6j.— 1 J-^^^ ni " )-yzc7hl4 

M izM \ or V iz iz/ ^ dce!16 Urec Trecz. / 



Case specific 



J"yzch4. : = 1 1 - kgpbi6cl (Trecz. zrec6. )•— • 1.1 4 



EceIK6 Trecin ncasei 



dce!16 Urec Trecz / 
ncase iz 1 



yzch4. 



Case specific 



4«c7M4 =3.385*10 ° yzc7hl 4 =0.0005 yZ c7hl4 - yzcfvt 

- 0 y "-6.77-10^ 



U yzch4 nz =2.1646*10 10 yzch4 Q =0.005 



yzc7h!4 



yzch4 



= 4.3292'10 



o 




Conversion of C7H14 under case 5 
(full power) is 99.3% [vs. 99.6% via 
averaged integrated formulas] with 
about 1 inch entry length. Conversion 
of CH4 (and also CO) is very high 
(>99.99%) within 3.5 inches. 



This result shows that the averaged 
integrated values for CH4, C7H14, and 
CO conversions are very good 
approximations and can be used for 
different values of fuel and CO 
Case specific ncase = 5 concentrations. 
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***^*****************************************************^********************* 



Add finite rate kinetics for PdO/AI203 methane combustion catalyst 
Begin by defining a few constants and catalyst functions 



R=82.076'K 1 -atmcm 3 



atm*cm3/mol*K 



Rem = 83143000-gmcm 2 -sec 2 K 1 



TOL ^0.0001 



Define Catalyst Functions: 



risph(+) 



3/ 1 1\ 



$ \tanh($) $ 



J 9 



Catalyst Physical Properties for Pd/AI203 

pcatpart =2.5-gm-cm 
dpcat =0.0002- cm 



OScat^Vgm 1 



-3 



^Vjipore -0.15-cm -gm~ 



-^Lcat -0,0020'cm 



pjlppore : - 6 



V^pore 
Scat 



ejipore =0.375 



tUpore 3.0 



4ipore = 1.0588' 10 *cm 



DeffCH4(T) =■■— - 



en pore 
iH pore 



I 



; Dch4T(T) 



|^R-T 



d^pore* 



dfipore- 



■R.Tin 



= 0.1 095 *sec 1 'cm 2 



DeffCH4(Tin) =4.43271" 10 3 'sec 1 'cm 2 



Dch4T(Tin) - 1.2343'sec 1 'cm 2 



1 
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Def!C7H14(T) 



enpore 
^pore 



1 



| Dc7hl4T(T) 



d^pore* 



M 



■R-Tin 



dupore- 



, 71 

' MW, 



= 0.0829* sec 1 *cm 2 



DeffC7H14(Tin) = 3.20528 *10~ 3 *sec 1 'cm 2 



Dc7hl4(Tin) =0.356'sec 1 "cm 2 



DmacCH4(T) ^*Dch4T(T) 
DmacC7H14(T) = ^-Dc7hl4T(T) 



^ kOT - l.O-gm^sec 1 



DmacCH4(800*K) =0.193-sec~ l -cm 2 



DmacC7H14(800-K) = 0.0556"sec 1 W 



Test value for kOT 



i Microscopic Thiele Modulus for CH4 



yiWtT) - [pcatpart 



mi 



\Pref/ DeffCH4(T). 



dpcat 
6 



<hi(kOT,Tin) =0.2042 



IB Macroscopic Thiele Modulus for CH4 



I— )• 

\Pnfl 



ly <t>mac(k,T) ' = J pcatpart- 1 

s:S5» 

Microscopic Thiele Modulus for C7H14 



nsph(fr(k,T)) Uat 
DmacCH4(T) 



+fi7(k,T) 



j R-T , 
pcatpart-^— fj . 



Def¥C7H14(T) 



dpcat 
6 



<*mac(k0T,Tin) = 1.8333 



♦(i7(kOT s Tin) =0.2402 



Macroscopic Thiele Modulus for C7H14 



W7(k,T) = r^^^I--kl- ^ 7 ^ T)) -Lc«t 
\f \Pref. j DmacC7H14(T) 



R-T 



RatekCH4(k,T) -T) S ph(Mk,T)) Tlfpl(*mac(k,T))-| pcatpart-) " \ k> 

' Pref i ' 



R-T 



tJ>rnac7(k0T, Tin) =3.4117 



RatekCH4(k0T,Tin) = 85999.49 14* sec 



RatekC7H14(k,T) = nsph(^7(k,T))-nfpl( < I>mac7(k > T))* ; pcatpart-: — - *k 

Pref- 



i 



RatekC7H 1 4( kOT, Tin) = 48479.296*sec 
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************************************************ 



Defines transport restricted rate constant function given cataiyst effective first order rate constant 
kOT, methane or cyciohexane concentration, x, and surface temperature, T 

PdO Catalyst Rate Constants 

ErpdO M2800-K Erpdl =4900-K Erpd3 = 15000-K 
RegtonO - iow T PdO to 625K Regionl - Intermediate T Region2 - Fiat 909K Region3 - 
ignores inhibition by H20 and C02 625K to 909K to 1048K decomposition Pd metal 



ApdO : = 2.4-10 3 -200-gm" 1 -sec" 1 



. i ErpdO Erpdl , 

Apdl =ApdO-exp(-^.— J 



Apd2 =A P dl.exp(^|) 



Apd3 -Apd2-exp 



Erpd3 
1048.1-K 



Interpolation formula for the four PdO+Pd rate constants 



itpd(T) 



, M 2 i 2-Erpd3\ 
Apd3 2 'expi + - 



1 



0.5 



A P d2" 2 , Apdl^exp^ j + ^0-2.^(2^0 



z si *** * * ** ***************************** 



* *********+**********************' 



Chx . = 0,1. .41 



T = 500-K 4- nx-25-K 



CATALYTIC RATE CONSTANT vs TEMPERATURE 
(mol/atmCH4/s/gcat) 




CAPSTNE5.MCD 
keflfrec6cl(T,z) =- 
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kgpb!6cl|Tavg ncase ,— j = 35.3023*sec 1 'cm 

\ 2 / 



\RatekCH4{kpd(T),T>Lcat/ kgpbKcl(T.z) 



keffrec6cl | Tav^^,— ! -0.7922* sec 1 'cm 



kefftec6c7(T,z) 



RatekCH4^kpd(Tavg ticase j,Tavg ncgse j*Lcat ^0.8104-sec *cm 



kgpb^lTavg^, — j ^10.1756-sec" 1 *cm 



\RatekC7H14(10-kpd(T),T)-Lcati ^ kgpb!6c7(T,z) 

RatekC7H14(l0-kpd(Tavg n ^),Tavg ncase )*Lcat : =7.3807'sec 1 'cm 



Lrec 
ticasc' 2 



keflfrec6c7^Tavg nf _,^— ) =4.2778* sec ' 'cm 



10 



kcf&ec6cl (Trecz^z^^ 

RatckCH4 (kpd (Trecz^ , Trecz^ Lcat 0 j 

kgpbl6ct (T^z^.^c 6 ^ 



001 



no 




Case specific ncase = 5 



keffirec6c7 .Trecz^, zrec6 ^ 

RatekC7H14 (W kpd(Trecz Jz/ ,Trecz |Z ; Lcat 0 j 

kgpbI6c7 (Trecz^, w«c6 lz ) 



001 



0.05 


0 

zrec 


1 


0.15 










































I - 1 






















**■" ■ — ■ 














*" * K 
















^^^^^ 





















0.2 



0.05 



0.1 
zrec6 , 



0.15 



0.2 
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Determine conversion of fuels with transport-limited finite catalytic combustion rates: 



Trecz. -Trecin 

iz ncase 



Trecin - Trecout 

_ ncase ncase 

nz 



Trecz = 630.3722* K 

nz 



Case specific ncase = 5 



Trecout = 630.3722" K 

ncase 



yzkc7M4 =yc7h!40 



yzkch4 -ych40 



f T rer 4 £reH6 Trecin- \ 

yzkc7h!4. . = 1 1 - keffrec6c7( Trecz. ,zrec6. V— ■ -yzkc7hl4. 

J a-H V "/ nz dce!I6 Urec Trecz.' lz 

\ ncase iz 



yzkch4 jz ^ } - j I - keffrec6c I ( Trecz. z , zrecd J 1 . 1 * 



Lrec t , 4 seell6 Trecln 5 j t UA 
11 . -yzkch4. 

nz dce!16 Urec Trecz. / a 



Case specific 

ncase - 5 

Case specific 



fyzkcVhH^ - 7.036* 10 yzkc7ht4 0 = 0.0005 yz kc7hl4, 
I yzkch4 nz s = 0.0034 yzkch4 0 =0.005 



yzkc7hl4 Q 



_ yzkch4 
m = 1.407' 10 1 — 55 -0.6838 



yzkch4 



yzlcc7hl4j 2 
yzfcc7ht4 0 0.1 

yzc7h!4 jz 
yzc7hl4 0 

yzkch4 iz 
yzkcM n 



0.05 



0.1 

zrec6, n 



0.15 



0.2 



Case specific ncase = 5 

Conversion of C7H14 under case 5 
(full power) is 99.4% [vs. 99.3% for 
the maximum BL transport-limited 
conversion] with gradually decreasing 
(rate limited) conversion. Conversion 
of CH4 is very poor (35%). CO 
conversion is probably very high given 
higher catalytic and transport rates 
than with C7H14. 

This result shows that the specific 
rates of catalytic oxidation of all added 
components needs to be measured as 
applied to a recuperator. The use of a 
recuperator for emission control would 
appear to be greatly enhanced with 
higher exhaust (recuperator) 
temperatures. 



In the following section, the conversion of C7H14, CH4 and CO is considered for various honeycomb monliths 
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ThC jcase t 1 = Tredl W'I 



Ghc2 



jcase t- 1 



mf. * 
jcasc * 

Shc-Ecell2 



mf 



Ghc4 



jcase +-1 



jcase ■ 1 Shc-scell4 



Ghc6 



mf f . 

jcase +■ 1 



Jcaseil Shc-ecell6 



Uhc2 = 



; o 

| 3.058 
3.9753 
4.80! 
j 5.8101 
! 10.0353 



"nrsec 



Uhc4 



0 

3.2157 
4.1804 
5.0487 
6.1098 
10.553 



"m'sec 



Uhc6 = 



0 

3.1706 
4.1218 
4.9779 
6.0242 
10.4051 



•nrsec 



lEvaluate Reynolds numbers for use in calculating the momentum, heat, and mass transfer functions: 



j Rehc2 



jcase +- 1 



Ghc2 i Cm c + r dce112 

mnixT(Thc jcase + 1 ) 



Rehc4 



, GW ,«t.' dw114 

\mixT(Th W + 1 ) 



Rehc6. 



Ghc6. -dcel!6 

jcase + 1 



Rehc2 : 



0 

43.776 

56.9088 

68.7283 

83.1744 

167.4548 



Rehc4 = 



0 

31.7427 

41.2655 

49.836 

60.3111 

121.4242 



Rehc6 = 



0 

25.7355 
33.4562 
40.4048 
48.8975 
98.4453 



cfrichc2(j,z) 



. 14.227 
Rehc2. 



i . . dcell2\ 
1 + 0.046263- |Rehc2 j - — -j 



0.45363 



cfrichc4(j,z) 



14.227 
Rehc4 



1 +■ 0.046263- I Retort 



dcel!4\ 



0,45363 



Average coefficients of friction 
for smooth square channels 



1 4 727 I . ^ dcel!6 \ 

cfnchc6(j,z) ~-7'l 1 + 0.046263* :;Rehc6- - - 
Rehc6. i_ J z y 



0.45363 
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APhc2. caseM " 4-cfrichc2(jcase t l,Lhc)- 



Uhc2 



j case 



■pmixhc 



Lhc 



J case4 ' dcel!2 



APhc4 



jcase 



Uhc4 .) 2 

tl =4-cfrichc4(jcase-H l.Lhc) Jcase " hi/ -pmixhc case ^ 



_Lhc_ 
'"dcelW 



APhc6. ^ M-cfrichc6(jcase-hl,Lhc)-^ JCMe + !/ -pmixhc. ■ Lnc 

jcase -H VJ 9 / 2 jcase -H 



dcel!6 



, APhc2 
: Phc 



0 

0.0013 

0.0017 

0.002 

0.0024 

0.0041 



APhc4 
Phc 



0 

0.0028 
0.0036 
0.0044 
0.0053 
0.0088 



APhc6 
Phc 



0 

0.004 

0.0053 

0.0064 

0.0077 

0.0127 



APhcmax P 



0 

0.003 
0.004 
0.005 
0.006 
0.011 



Chose 400 cells per in2 
honeycombs via estimated AP 



************* ********** * 



Use 400 ce!ls/in2 in honeycomb monolith for calculation of CH4, C7H14, and CO conversions 



Find mass transfer rates via Nusselt formulas: 



Numavghc4cl. 



jcase f I 



= 2.976- 



1 +■ 



7.6-Lhc 



Rehc4. -SclTfThc. 'i-dcell4/ 

jcase - 1 \ jcase + 1/ t 



Numavghc4cl 



0 

2.976 

2,9761 

2.9761 

2.9762 

2.9774 



Numavghc4co 



-2.976-i 1 + j _- 



7.6-Lhc 



jcase H — j - ; Rehc4 ^ SccoT / Thc , . dceU4 ; 

\ jcase -+- I \ jc - 



jcase +- 1 ' 



Numavghc4co = i 



0 

2.976 

2.9761 

2.9762 

2.9763 

2.9778 



*******************************^^ 



CAPSTNE5.MCD 



7.6Lhc 



■ 3 i 



Numavghc4c7. , " 2.976- M + | - 7— \ 

3 J case 1 ' Rehc4 -Sc7TfThc d . >-dccII4 



Once again the entrance effects within monoliths are generally small 



Numavghc4c7 - 



0 

2.977 

2.9781 

2.9795 

2.9817 

3.0123 



♦a*************************** 



Now determine the maximum possible CH4, C7H14 t and CO conversions in a 400 cells per in2 monolith 

ych40 -0.005 yc7hl40 =0.0005 ycoO =0.005 

: ych40expl 



ych4hc4. 



jcase -H 



4 * Lhc*Nuntavghc4c 1 



yc7hl40 =0.0005 
- 4*Lhc-Numavghc4cl 



jcase -H 1 



jcase -h 1 



! Rehc4 ica S e,.' SclT i Th V S e^i 



}-dcell4 





-48.7678 




-37.5142 


w - 


-31.0633 




-25.6688 




- 12.7326 



ych4hc4 
ych40 



0 
0 
0 

3.2315'10 
7.115M0 
2.9532*10" 



14 



"12 



4-Lhc-Numavghc4co jcase + j \ 



ycohc4 
ycoO 



0 
0 

1.4726*10 
5.2301*10" 
7.0997-10" 
9.1398-10" 



-15 



4LhcNumavghc4co i 



jcase +• 1 



Rehc4 jcase + 1' SCC0T ( ThC jcase + l) ' dCeI14 



44.3966 



34.1518 



-28.2792 



-23.3684 



11.6029 



24 



CAPSTNE5.MCD 



25 



, 4-Lhc-Numavghc4c7 
yc7hl4hc4. yc7hf40-exp| — J case t.J 



6.8527' lo" 7 
yc7h!4hc4 _| 1.8058" 10 5 

y c7hl4( > ! 0.0001 

i 0.0006 



4 • Lhc- Numavghc4c7 



jcase t- 1 



[ 0.0239 



14.19351 



10.9219 



9.0478 



- 7.482 



3.7348 



Only for C7H14 and only in case 5 (maximum power) does the BL transport-only 
calculated conversion (97.6%) fall below 99.9% with a 400 cells/in2 monolith. With 
50 ppm C7H14 in the exhaust in this case, the exit concentration would be 1 .2 ppm. 



For the final analysis, consider the finite catalytic conversion of C7H14 and CH4 in a monolith: mease = 5 



zzcharhc4c7. 



Rehc4. r Sc7T(Thc. ^ , VdceI14 



jcase +- I 



Lhc 



zzcharhc4c7 



= 0.6199 



zzcharhc4cl 



Rehc4. .SclT(Trecin. ^ 
. jcase + 1 \ jcase -f- 



jcase -f- 1 



j j*dcel!4 



Lhc 



zzcharhc4cl 



= 2.1382 



nz -400 



Lhc = 10.16 -cm dce!14 =0.1 168 xm 
zhc4 



Lhc 

zhc4. -iz ^ dwall4 zzhc4c7. r 

12 nz ,z zzcharhc4c7 



zzhc4c 1 =- 



zhc4 



Numono ^zzhc4c7- z j 
Numono(zzhc4cl 



lz zzcharhc4cl 




zzcharhc4c7 : 



0 i 
4.2621 i 
5.5407 
6.6914 ' 
8.0979 
16.389 i 



i'em 



zhc4 
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Case specific mease = 5 

Dc7hl4T(T)-Numonot Z 



zzcharhc4c7 / 

kgpblhc4c7(T,z) x - kgpblhc4c7;Thc ,Lhc) - 1 1.201 4*sec _l 'cm 

Ar-aMA ^ mease / 



dce!14 



Dch4T(T)'Numono[ 



^zzcharhc4cl , 

kgpb!hc4cl(T,z) — ~ — ! mcase/ kgpblhc4cl(Thc ,Lhc) -38.6 127" sec 1 'cm 

dcell4 mcase ' 



keffhc4cl(T,z) 



/ 



1 



\RatekCH4( kpd(T) , T) -Leal / kgpblhc4cl (T , z) 



1 



kgpblh<>tel(Thc^ f ^J-38.615-sec 1 -cm 



ke»c4cl(Thc ) - 1 .8533' sec 1 'cm 



RatekCH4(kpd(Thc ),Thc j-Lcat = I.9467*sec 1 'cm 

V r V mcase/' mcase/ 



'keffhc4c7(T,z) = — 

! i 



i 



RatekC7H14(10-kpd(T) ) T)Lcat/ kgpblhc4c7(T,z) 



1 



kgpbIhc4c7|Thc ncMe ,^J = lL2528*sec 1 'cm 



keffhc4c7|Thc ) = 6.7062'sec 1 "cm 

i 11X886 2 / 



RatekC7H14(l0-kpd(Thc mcase ),Thc fficase )-Lcat = I6.5977"sec 1 'cm 



kefB»c4c7(Thc mcase ,2hc4 iz ) 

RalekCTHU; 10kpd(Th Cmcase ) .Thc^^Lcat 0 , 
kgpbI6c7{Thc mcase ,zhc4 iz ) 



0.01 



0 0.02 0.04 0.06 0.08 0.1 0.12 

zhc4,„ 
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Determine conversion of fuels with transport-limited finite catafytic combustion rates: 
Case specific mease = 5 



yzkhc4c7hI4 Q = yc7hl40 



yzkhc4ch4 Q -ych40 



yzkhc4c7hl4 = [ 1 - keffhc4c7(Thc ,zhc4.}^ |-yzkhc4c7hl4 

a + l \ ^ mease izj nz dcel | 4 Uhc4 \ > u 



yzkhc4ch4- 2 ^ { 



= 1 1 - keffhc4cl (The ,zhc4 ^^\. yzkhc4ch ^ 

\ mease' «/ flz dce ,j 4 uhc4 ^ c ™ lz 

\ mease/ 

i 



The- 



0 

1005.3722 
1005.3722 
1005.3722 
1005.3722 
916.4833 



Case specific 



mease = 5 



yzkhc4c7hl4 nz = 7.281' 10 



yzkhc4c7hl4 n =0.0005 



yzkhc4ch4 nz = 0.003 



yzkhc4ch4 = 0,005 



yzkhc4c7hl4 



yzkhc4c7h!4 



™ = 1.456* 10 1 



yzkhc4ch4 n 
yzkhc4ch4 Q 



= 0.5959 



yzkhc^hH^ 
yzkhc4c7bI4 0 

yzkbc4ch4j z 
yzlcbc4ch4Q 




0.1 



0.01 



no 



Case specific mease = 5 

Conversion of C7H14 under case 5 
(full power) is 99.4% [vs. 99.6% via 
averaged integrated formulas] with 
about 1-inch entry length. Conversion 
of CH4 (and also CO) is very high 
(>99.99%) within about 3.5-inch. 

This result shows that the averaged 
integrated values for CH4, C7H14, and 
CO conversions are very good 
approximations and can be used for 
different values of fuel and CO 
concentrations. 



0.02 



0.04 



0.06 



0.O8 



0.1 



0.12 



zhc4 t . 
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Defines transport restricted rate constant function given catalyst effective first order rate constant, 
kOT, methane or cyclohexane concentration, x, and surface temperature, T 

PdO Catalyst Rate Constants 

ErpdO M2800-K Erpdl -4900-K Erpd3 -" 15000-K 
RegionO - low T PdO to 625K Regionl - Intermediate T Region2 - Fiat 909K 



ignores inhibition by H20 and C02 625K to 909K 



Region^ - 
to 1 045K decomposition Pd meta! 



ApdO ^2.4-10 3 -200gm '-sec' 1 



Apdl , A pd0.exp( -J».^L] A P d2 ^Apdl-exp'^ 



625-K 625-K 



909.1 -K 



Interpolation formula for the four PdO+Pd rate constants 

kpd(T) 



Apd3 -Apd2-exp 



I Erpd3 
\ 1048.1 K 



A pd3 2 -expf - 2 -^l + 



1 



0.5 



Apd2- 2 + Apdl^expf^i) + ApdO^exp/^^V, 



T / ~~ r ~~ ~~ r \ T 

Add expression for R(.5%)Pd(3.0%) catalyst with C3H8 here, use x2 for estimated rate constant with C7H14: 
Rate constant via crushed powder results for C3H8 



kptpdC3(T) : = 6.9-10*-exp(- ^^H gm-'-sec' 1 



10kpd(900K) =49.2392-kg '-sec 1 



2-kptpdC3(900 K) =886.4358*kg Sec 1 



nx:-0,1..41 



^( T nx) 
2.kp*dC3(Tj 



T = 500K + nx-25K 

nx 



CATALYTIC RATE CONSTANT vs TEMPERATURE 
(mol/atmCHx/s/gcat) 




JOOOK 
T 
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keffrec6cl(T,z) =- 



\RatekCH4( kpd(T) , T)-Lcat/ kgpb 16c I (T, z) 



kgpbWclfTava .^l-353023-scc 1 'cm 



keffirec6cl j Tavg^ , ~ J = 0.7922- sec" ' "cm 



keffrec6c7(T,z) 



\RatekC7H14(2.kptpdC3(T),T)Lcat/ kgpbJ6c7(T,z) 



RatekCH4(kpd(Tavg ncase j .Tavg^j -Lcat = 0.8104'sec 1 'cm 

kgpbiecT^Tavg^.-Jpj = 10.1756-sec - ' -cm 



keflfretfcTjTavg^,—] =7.4067-sec 1 -cm 

\ 2 / 



RatekC7H14(2-kptpdC3(Tavg ncase ),Tavg ncasc ).Lcat=27.2196'sec" 



kefficc6cl (Trccz jzf zrec6 iz ) 
RatekCH4 (kpd (Trecz jz ) , Trecz J2 ) teat 0 { 
kgpbI6c ! (Ttscz ir zrec6 i2 ) 



0.0 1 



1-10 




Case specific ncase = 5 



keffiec6c7 (Ttccz^, zrccd^ 

RatekC7Hi4i 2kptpdC3 (Trccz jz ) .Trecz^Lcat 

kgpN6c7 ( Trecz. t zrec6 » 
V iz m 




0.01 
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Determine conversion of fuels with transport-limited finite catalytic combustion rates; 
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Trecin Trecout 

Trecz. Trecin - ncase - n «* j z 

nz 



Case specific ncase = 5 



Trecz - 630.3722" K Trecout = 630.3722* K 

ILC ncase 



yzkc7hl4 ^yc7hI40 



yzkch4 0 -ych40 



nz dcel!6 Urec Trecz. ' 

ncase tz / 



yzkch^^j := 1 - keffrec6cl(Trecz jz ,zrec6 iz )*— IA 



Lrec , t 4 eceI16 



Trecin, 



nz dceI16 Urec Trecz. 



— -yzkch4 



Case specific 

ncase = 5 

Case specific 



yzkc7hl4 nz = 1.802" 10 5 yzkc7h!4 0 = 0.0005 yz kc7hl4 



yzkch4 nz = 0.0034 



yzkch4 =0.005 



yzkc7h!4 ft 



= 3.603*10 



-7. 



yzkch4 M 
yzkch4 



-0.6838 



yzkc7hl4 iz 
yzkc7hl4 0 0-1 

yzc7hl4- z 
yzc7hH 0 

yzkch4 iz 
yzkcM^ 




0.01 



no" 3 



0.05 



0.1 
zrec6>. 



0.15 



0.2 



Case specific ncase = 5 

Conversion of C7H14 under case 5 (full 
power) is 97.4% with estimated catalytic 
rates and transport limitations {vs. 99.6% 
for the maximum BL transport-limited 
conversion]. Conversion of CH4 is very 
poor (32%). CO conversion is probably 
very high given higher catalytic and 
transport rates than with C7H14. 



This result shows that the specific rates 
of catalytic oxidation of all added 
components needs to be measured as 
applied to a recuperator. The use of a 
recuperator for emission control would 
appear to be greatly enhanced with 
higher exhaust (recuperator) 
temperatures. 



In the following section, the conversion of C7H14, CH4 and CO is considered for various honeycomb monliths 
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The. 



J case t I 



Trecin 



jcase • I 



pmixhc pmixrec[ Trecin . . 

jcase ! I \ jcase t \j 



Ghc2. 



'jcase t I 



mf 

jcase -_l 

Shc-EceiI2 



Ghc4 , = + -± Ghc6 . 

JcaserI Shc-£cell4 Jcase + ! 



mf 

jcase «-_[ 

Shc-eceli6 



Uhc2 = 



0 

3.058 
3.9753 
4.80 1 



'm'sec 



I 5.8101 
l 10.0353 



Uhc4 = 



0 

3.2157 
4.1804 
5.0487 
6.1098 
10.553 



m'sec 



Uhc6 = 



0 

3.1706 
4.1218 
4.9779 
6.0242 
10.4051 



"m'sec 



Evaluate Reynolds numbers for use in calculating the momentum, heat, and mass transfer functions: 



Reh^^W^ Rehc4 . Gh «W**"4 Rehc6 = Gfa «W +l **" 6 



^Rehc2 = 



0 i 
43.776 i 
56.9088 • 
68.7283 ! 
83.1744 i 
167.4548 j 



Rehc4 = 



0 

31.7427 

41.2655 

49.836 

60.3111 

121.4242 



Rehc6 : 



0 

25.7355 
33.4562 
40.4048 
48.8975 
98.4453 



cftichc2(j,z) =!^ 2 ? 7 .! l +■ 0.046263- (Rehc2,— 12 -) 
Rehc2 \ i z 1 

i 



0 45363 



cfrichc4(j,z) - ^^ 27 J | +. 0.046263- (Rehc4 -~~ 4 - i 
Rehc4 I \ J z / 



0.45363 



Average coefficients of friction 
for smooth square channels 



cfrichc6(j,z) » *? 2 - 7 . 1 * 0.046263- iRehc6- dc - 11 - 6 ■ ' 
Rehc6 J z 



0 45363 



CAPSTNE6.MCD 
Calculate monolith pressure drop in recuperator with entrance effect: 
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APhc2. 



Uhc2 



jcase h 1 ) 



Lhc 



jcasc 



--4*cfrichc2(jcase t l,Lhc)~ J ~~— -pmixhc. 



, Uhc4 , ) 2 T , 

jcasetl 2 jcase+ - 1 dcel!4 



APhc6 jcase +■ l : = 4 * cfri *chc6(jcase +■ 1 , Lhc) 



\_ jcase -H/ . 



pmixhc. 



Lhc 

^ ase ^r dcell6 



APhc2 
Phc 



0 

0.0013 

0.0017 

0.002 

0.0024 

0.0041 



; 0 

I 0.0028 
APhc4 I 0.0036 



Phc 



! 0.0044 
! 0.0053 
L 0.0088 



| 0 
0.004 

APhc6 0.0053 



Phc 



0.0064 
0.0077 
0.0127 



3.5 APhc6 
4 Phc 



0 

0.0035 
0.0046 
0.0056 
0.0068 
0.0111 



APhcmax P - 



0 

0.003 
0.004 
0.005 
0.006 
0.01 1 



Chose 600 cells per in2 honeycombs but 
reduce length to 3.5in via estimated AP 



Lhc -3,5-in 



* Use 600 cel!s/in2 in honeycomb monolith for calculation of CH4, C7H14, and CO conversions 
Find mass transfer rates via Nusseft formulas: 



Numavghc6cl. casc+ ( =2.976- 



i ^ \ 7.6-Lhc 

Rehc6_ -SclTfThc. .)-dceII6 

L \ jcase f 1 V jcase + 1 / 



Numavghc6cl = 



0 

2.976 

2.976 

2.9761 

2.9761 

2.9767 



Numavghc6co. u . - 2.976-; 1 ^ j 

JcaseM ! lRehc6 ■ 

* jcase f- 1 



7.6Lhc ■ 

•SccoT[Thc. , -dcell6 ' 

V jcase * 1 ' / 



Numavghc6co : 



0 

2.976 
2.976 
2.9761 
2.9761 
[ 2.9769 
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Numavghc6c7 jcjse ^ j = 2.976 



7.6Lhc 



' l Rehc6 jca Se ^ C7T ( T Sa^7H Ce116 

Once again the entrance effects within monoliths are generally small 



Numavghc6c7 : 



0 

2.9765 
2.977 
2.9777 
. 2.9788 
= 2.994 
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Now determine the maximum possible CH4, C7H14, and CO conversions in a 600 cells per in2 monolith 

ych40 =0,005 yc7h!40 =0.0005 yco0^0.005 

- 4*Lhc-Numavghc6cl . 



y ch4hc6 jcaseti = ych40-exp 



jcase + I 



Rehc6 jc a s.^- Sc,T ( Th «= jc , a e^)- dcel16 



'% - 4-Lhc*Numavghc6c 1 . 

yj Rehc6. amm -SclTfTha ^ , )dcell6 

jcase -f- 1 \ jcase -H/ 





-64.0074 




-49.2367 




-40.7696 


yj 


-33.689 




-16.7076 



!ycohc6 



jcase -h I 



: yco0-exp 



4 * LhcNum avghc6co. 



jcase +• 1 



Rehc6 J c«e + .- SCCOT ( ThC jca S e4- 1 )- dcel16 



ycohc6 
ycoO 



0 
0 
0 
0 

4.7906*10" 
2.445* lO^ 7 



14 



ych4hc6 
ych40 



0 
0 
0 
0 

2.339-10 
5.5463-10 



15 



4 -Lhc-Numavghc6co 



jcase -f I 



Rehc6 ^SccoTiThc 

jcase -HI 



tine . j-dceIJ6 



58.27 



44.8234 



37.1153 



30.6695 



-15.2241 



JAPSTNE6.MCD 
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■ r - 3.5*Lhc*Numavghc6c7 
yc7hl4hc6 jcasetI -yc7hl40-exp J - C *^I 



3.5 -Lhc-Numavghc6c7. 



\Rehc6 -Sc7T(Thc (J -dcell6/ 

\ jcasefl \ jcaseiT / R e hc6. -ScTT The. 



j case 1 I 



j case 



jj-dcelie 



yc7h!4hc6 
yc7hl40 



8.3588*10 
3.5856* 10~ 
3.0905-10 

j 0.0002 

|_ 0.0141 



5 



16.2974 



12.5386 



10.3846 



-8.5841 



4.2631 



Only for C7H14 and only in case 5 (maximum power) does the BL transport-only 
calculated conversion (98.6%) fall below 99.9% with a 600 cells/in2 monolith With 
SO ppm C7H14 in the exhaust in this case, the exit concentration would be 0.7 ppm. 



For the final analysis, consider the finite catalytic conversion of C7H14 and CH4 in a monolith: 



mease =5 



zzcharhc6c7. 



jcase-f- 1 



Lhc 



zzcharhc6c7 



= 0.8137 



zzcharhc6cl. 

jcase-H 



■ „ Rehc Wr SclT ^ 



Lhc 



zzcharhc6cl 



-2.8064 



nz = 400 



Lhc = 8 .89 -cm dceI16 = 0.096 1 'cm 
zhc6. 



iz 



zhc6 - \z^^- -t dwal!6 zzhc6c7. - 

nz lz zzcharhc6c7 



zzhc6cl . 



zhc6 



zzcharhc6cl 



Numono (zihefcl-j 



Numono (zzhc6c 1 ^ 




no 



zzcharhc6c7 = 



0 

2.8414 
3.6938 
4.461 
| 5.3986 
[ 10.926 



*cm 



zhc6- 
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26 



Dc7hl4T(T)Numono 



/ 



kgpbihc6c7(T,z) 



kgpblhc6cl(T,z) 



zzcharhc6c7 



dce!16 



Case specific mease = 5 

, — c j 

mease' / 



kgpb!hc6c7fThc ,Lhcj -13.6172-sec 1 "cm 



Dch4T(T>Numonoi 



zzcharhc6cl \ 

mease/ 



dcel!6 



kgpb^c^Thc^^Lhc) -46.9581 -sec 1 -cm 



kef&c6cl(T,z) 



1 I 



kgpblhc6cl|Th Cmcjise ,^j =46.9589-sec~ 1 -cm 



UatekCH4(kpd(T),T)-Lcat/ kgpblhc6cl(T,z) 

RatekCH4(kpd(Thc inca J,Thc mcase ).Lcat = 1.9467-sec~' "cm 



kef&c6cliThc ,^)=1.8693-sec ' -cm 

y lUCSSC 2 



keffhc6c7(T,z) 



1 



kgpbIhc6c7[Thc^, — ] - 13.6468'sec 1 'cm 



RatekC7H14(2.kptpdC3(T),T).Lcat/ kgpblhc6c7(T,z) 

RalekC7H14(2-kptpdC3(Thc J -ThcJ-Lcat = 120.0371 -sec" 1 'cm 



Lhc 



kefftc6c7 Thc^ =12.2537"sec 1 'cm 

^ mease 2 



10 



]^ hc6c7 ( Thc m C a*' zhc6 i z ; 

_R S tekC7Hi4(2.kpt P dC3(Thc nlcase ) .Thc^-Lcat , 
kgpNecTiTlK^^.Ehc^ 



0 1 1 



















































































\ 









































0.02 0.04 0.06 

zhc6 iz 



0.08 0.1 
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Determine conversion of fuels with transport-limited finite catalytic combustion rates: 
Case specific mease = 5 



yzkhc6c7h!4 yc7h!40 



yzkhc6ch4 -ych40 



* mease/ 
* mease/ 



Thc = 



1005.3722 
1005.3722 
1005.3722 

f 

j 1005.3722 
I 916.4833 



Case specific 



mease = 5 
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yzkhc6c7h!4 n2 = 9.964- 10 



yzkhc6c7hH 0 =0.0005 



yzkhc6ch4 nz - 0.0028 



yzkhc6ch4 Q = 0.005 



yzkhc6c7hl4 
yzkhc6c7h!4 



112 = 1.993' 10 2 



yzkhc6ch4 a 
yzkhc6ch4 



= 0.5693 




Case specific mease = 5 

Conversion of C7H14 under case 5 (full 
power) with 600 cellsrfn2 x 3.5-in length 
and extrapolated catalytic rates (via 
C3H8 at 400C) is 98.0% [vs. 99.6% via 
averaged integrated formulas, 99.3% 
with temperature dependent transport, 
and 98.6% with 400 cells/in2 and 4-in 
length. Conversion of CH4 is limited by 
slow kinetics and CO is expected to be 
very high (>99.99%). 



zhc6.„ 
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Gas Properties and Constants 



cm=0.01-m 
N = kg-m-sec~ 



R = 8.3143- 



J = kg-m 2 -sec~ 2 
Pa-N-m" 2 



1 



W=J-sec 
atm = 101300- Pa 



Rc =8.3143-J-mole- , .KT 1 



***********-************ 



Pref = latm 
^im= 0.00000 1-m 



n Gas 

1 CH4 

2 02 

3 N2 



MolWl 



G(LJ) 



e(U) 



MW Q ^16.044 gm a Q = 3.785-10' 10 m z Q = 148.6K 
MW, =3l!999-gm a, = 3.464- 10" I0 -m e, = 106.7K 
MW 2 = 28.018gm a 2 =3.798* 10" 10 m e 2 = 71.4-K 



4C02 MW 3 =44.010-gm a, = 3.94M0 _I -m 8 3 = 195.2-K 



5H20 



MW 4 = 18.015-gm a 4 - 2.641* 10" 1 -m e 4 = 809.1 K 
*********^******************^^ 

Define transport integrals: 
Constants: 

=1.16145 anb 0 = 0.14874 aDa 0 =1.06036 oDb Q =0, 

a\xa y =0.52487 



1561 



a^bj =0.7732 aDaj =0.193 aDb, =0.47635 

0^=2.16178 anb 2 =2.43787 00^= 1.03587 aDb 2 = 1.52996 

orDa3 =1.76474 aDb 3 =3.89411 



Functions: 

IWT) = a^-expf (^b 0 .In(T)) ] + £ H«V0 0394 

i|X = 1 



3 

QD(T) ^aD^exp^aDbolnCT))]^ J] 0Da iD -exp(-aDb iD T) 

iD= 1 



2 

S 

in- 1 

3 

E 

iD= 1 



QD(3) =0.95002 
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***** *******^***************iMr***^H******#******* ****************************** 

Set Inlet Conditions: 

P =1.0-atm ych40 =0.035 

Define gas composition w/r methane using moist (3%) air; 



y(x) 



( 1 - ych40)0.2 - 2-ych40 + 2x 
( 1 - ych40)0.77 
ych40 - x 
( 1 - ych40)-0.03 + j2-ych40 - 



y(0.02)« 



0.02 

0.163 

0.743 

0.015 

0.059 



♦♦♦a************************************************************************************************************ 



Define viscosity function 



H<T) = 



f26 (MWo-T)^ lO^cm /gmf 5 
sec K 



2.6693-10 











(MW 2 -T 


I 0 - 5 








0, 



2.6693- 10 



/10 s . 



cm 



\ sec 



2.6693- 10 



.acflO^cmVfgm' 0 - 5 
sec 



2.6693- 10 



. 26 /10 5 -cm\ /girO 05 



1.116-(MW 4 -Ty 



0.5 



(3.11510- 8 -cm 



^(300-K) - 



1.103-10" 
2.06* 10" 5 
177- 10" 5 
1.519-10" 
1.068*10" 



•kg-m l -sec 1 



CPSNALT2.MCD 



3 



Hmix(x,T) = ^] 

n = 0 



y(x)„(KT)) n 



m = 0 



1 t- 



M(T) n \ 0S p; oa 



8-1 + 



0.5 



y(0.02) = 



utriix(0.02,300 K) = 1.759M0 5 -kg*m '-sec 1 

Define thermal conductivity functions 
igas -0,1. .4 



0.02 

0.163 

0.743 

0.015 

0.059 



mwt***m»*t****** ************** ************************** 



aA. = 



1.869^ 




0.08727 


- .3273 




0.09966 


0.3919 


bk - 


0.09816 


-7.215 




0.08015 


7.341 




-0.01013 



L17910" 4 

■ (3.743* 10" 5 ) 
(5.067- 10* 5 ) 
5.477- 10~ 6 



1.801- 10 



r4 



dX = 



- (3.614- 10" 8 ) 
9.732- IO" 9 
1.504- 10~ 8 

-(1.053-10" 8 ) 
(9.110" 8 ) 



CH4 

02 
N2 

C02 

H20 



X(i,T) 



Amix(x,T) - 

i = 0 



0.00001 



W 
cm-K 



y(x) r X(i,T) 



W/cm/K 



4 

j = 0 


MTU 05 [MW^ 025 ' 

1 + [n(T)jJ [mw,] 


2" 




r / mw.\i 

8 r MW j)J 


0.5 



A<igas,1000*K) 


0.16716-kg*m«sec 


3. K -1 


0.07163-kgm*sec 




0.06292-kgmsec" 


3. K -I 


0.06788 -kg- msec 




0.0863 1- kg- m- sec" 





Xmix(.0,1000-K) =6.64*10 



cm-K 



Qt * M5N / ^*3^ ?l6fQ&* ********* 
***************************************^ 



Define thermal conductivity of alumina and SS321 



f 83.71 
! 224.898 
A!tc =i 255.842 
! 132.727 
I 25.9695 



The parameter set Altcji was determined from published values of the 
thermal conductivity. Ref: HOP 

>wAI203(x) is the thermal conductivity function (of x with unit K) while 
XI is its derivative w/r to temperature. 



n = 0 



XA1203(x) = £ A.tc, .-^ XI(x).= £ ^-n- U ._ 



n- 1 



S32ltc 



4.69312 
4.9812 
-0.33357 

0 

0 



The parameter set S321tc_n was determined from published values of the 
thermal conductivity. Ref: Handbook of Metals 8th Ed. Vol. 1 p. 423 



kss321(x).= ]T S321tc n | 



/ x 



n = 0 



W 



1000-Ki m-K 



W 

Ass321(1000-K) = 9.341 ♦— 
nvK 



XAI2O3(1000-K) -7.897- 



W 



X\ ( 1000-K). 1000-K = -7.517 •■ 



W 
rrvK 



i =30,31. .200 



Ta : = lOi-K 



XA!203(Ta^ 20 



M(Ta.) 1000-K 
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x := 1500-K 



xO root(M(x),x) 



xO = 1.588*I0 3 -K 



Xss32I(300 K) =6.157' 



W 
m K 



0.003- in =7.62' 10 '•cm 



XA!2O3(300-K) =0.359' 



W 
cnvK 



XA12O3(200 K) = 0.479* 



W 
cm*K 



AAI2O3<800«K) =0.102" 



W 
cm-K 



^AI2O3(x0) =0.054* 



W 



Minimum thermal 



cm-K conductivity 



Determine thermal conductivity of porous wall: 

Ref: RHG - Handbook Chpt 6 for effective thermal conductivity including 
convection and radiative terms 



erad^O.6 e =0.42 



dp = 0.00005 cm 



dp =5*10 7 -m 



Xwc<x,T) ^mix(x,T).|0.093 U 2.0.9069-f, A ,^ A T H T \ JU^ll 
i 1^AAI203(T)- Ajnix(x,T)j \ \Xmix(x,T) 



+ 0.229- 10 6 . -^-.g-dp srad T 3 - 
m 2 -K 4 



e + ■ 



1 -£ 



0.229* 10" 6 — ^— r-e dp-EradT 3 
m 2 -K 4 



*A1203(T) 



/ J 



A,wc(O.Ot,1000K) =0.473* 
Xmix(0.01,l000-K) -0.067 



Xwcnr(0.01,1000 K) =0.473* 



W 

m*K 

W W 
• AAI2O3(700-K) = 12.354- - 



W 
m-K 



Radiation component 
negligible with such small 
particles in the washcoat 



Xss32l(1000 K) =9.341 



W 
mK 
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MJ2O3(1000K) 
Amix(0.01,l000 K) 



= 117.74 



A.wc( 0.01, 1000 K) 
fcss321(1000K) 



= 0.051 



XwcCO.OUlOOO K) 
\mix(0.01,1000K) 



The thermal conductivity of the wash 
= 7.046 coat is about 7 times greater at 1 0OOK 
than the gas phase conductivity 



\wc(0.01 , 1000-K) Xss321( IOOOK) 



0.001 cm 



0.003- in 



Xss321(1000K) 
0.003- in 



= 1 .386 



The washcoat increases 
the thermal conductance of 
the wall by about 39% 



Xss321(1000K) 
Amix(0.01,1000*K) 



= 139.275 



For a Nusselt no. of 2.976 (dcel(8 = 0.0961 *cm) for developed flow in the recuperator (full power) the 
mean heat transfer boundary layer thickness is much greater than either the SS321 wall thickness or 
the washcoat thickness. 



Sgprecup = 



0.0961 cm 
2.976 



Sgprecup = 0.032 *cm 



hrecavg(T) 



\ 



1 



Amix(0.0],T) "~ Ass32l(1000«K) 



\ Sgprecup / 



0.003in 



Overall local heat transfer coefficient 
for the recuperator (assuming gas phase 
conductivities are about the same 



hrecavgwc(T) 



f 



Xmix(O.Ol.T) 



Sgprecup - 



2 0.001 cm 
2.976 



1 



\ 



Effect of wash coat on heat 



1 



transfer coefficient 



Xmix(0.01,T) r /Xss321(1000 K)\ r Xwc(0.01,T) 
Sgprecup / { 0.003- in J 0.001 -cm 



hrecavg( 1000- K) = 0.01038 



W 
cm 2 -K 



hrecavgwc( 1000-K) =0.01046 



W 



cm 2 -K 



The wash coat makes about 0.1% increase in heat transfer 
coefficient because of the slightly decreased cell size and greater 
thermal conductivity of the washcoat vs. gas phase. 




COMBUSTION CATALYSTS FOR 



Introduction 




)I described possible alternative combustion catalyst configurations and 
detailed analysis of relevant performance factors, i.e., mass-transfer limited hydrocarbon 
conversions, pressure drop, catalyst temperatures, post combustion conversion of hydrocarbons 
and CO, etc. This report represents a comprehensive collection of my findings and comments. 
Attached are the final MathSoft Mathcad 6 analyses of the suggested alternative combustion catalyst 
configurations. 



Alternative Combustion Catalyst Configurations 

Two approaches are suggested for the catalytic combustion of gasoline vapors in air. 
These approaches represent more conventional and perhaps less expensive alternatives to the PCI 
Microlith™ technology. They are summarized as follows: 

• Monolith with bypass air — A standard 3.5", corrugated FeCrAl alloy, 600 
cells/in 2 monolith should be able to reach « 86% (naphtha) to >99.8% (natural 
gas) fuel conversions. Using an air by-pass would allow longer residence time 
and higher adiabatic temperature rise (which is probably necessary for natural 
gas but may not be for diesel fuel or naphtha) within the combustor to ensure 
fast catalytic reaction rates. Higher temperatures within and at the exit of the 
catalyst monolith would also increase homogeneous combustion reactions, and 
help convert residual fuel vapors in the region immediately downstream of the 
catalyst but upstream of the turbine. A controller for bypass flow could help 
with turn-down, startup, and tuning steady-state conversions. 

• Radial flow catalyst bed « Provided that the flow rate can be reduced to 
< 3 m/s, a thin packed bed of catalyst beads or granules may also achieve 
conversions similar to those of the monolith and PCI Microlith™ catalyst screens 
at equivalent pressure drop. The bed might intrude upon your fuel injection and 
premixing zone, but this could actually benefit combustion (within limits, e.g., 
±30% fuel/air ratio) by providing locally higher adiabatic temperatures. A 
packed bed has fairly good thermal conductivity and is virtually immune to 
thermal shock problems. A screened frame and possibly metal plates (or 
alternatively for robustness, silica-coated silicon carbide) could be located within 
the bed to guide the flow. Downstream hydrocarbon and CO conversion also 
would be necessary for low emissions as in the case of the monolith. 



• Exhaust converter A 400 cells/in 2 post combustion monolith located just 
ahead of the regenerator exhaust stream inlet could provide a very significant 
reduction in hydrocarbon emissions regardless of the configuration and 
performance of the combustion catalyst. An emission control monolith should 
remove about 99% of the residual unburned fuel components and > 99.9% of 
the CO. 




Analysis of Catalyst Performance 

Mathsoffs Mathcad 6.0 Plus Professional program was used to perform detailed 
calculations of pressure drop and conversion of Cj, C 5 , C 7 , and C9 hydrocarbons (methane, 
isopentane, cycloheptane and nonene) and CO. The analyses generally assume catalytic rates are 
fast under Capstone's specified combustion conditions. The analyses use various catalyst 
parameters such as particle/channel size and bed length for the two upstream (monolith and fixed 
bed) configurations. Only monoliths were considered for the downstream 1-atm conversion of 
residual CO and hydrocarbons. 

Gas Properties 

Initially the Mathcad document describes and defines various gas component and mixture 
properties and functions (pp. 1-13) that are used in the subsequent calculations. The gas and vapor 
components are identified with their physical parameters (p. 1) and the composition of combustion 
mixtures as a function of fuel conversion are defined (pp. 2-3). Mixture densities, diffusivities of 
individual fuel component in the mixtures, viscosity, and thermal conductivity are then defined and 
plotted as functions of composition and temperature (pp. 4-9). The composition of various fuel 
mixtures is shown to minimally effect mixture physical properties (< 0.5% under such lean 
conditions). These properties are then redefined as empirical functions of temperature for 
simplification of subsequent analysis (pp. 6-9). Heat capacities, heats of combustion and enthalpy 
functions of the mixtures are similarly defined (pp. 9-11), and the initial composition necessary for 



an adiabatic temperature rise of 511 °F are determined for each of the four fuel components (p. 1 1 
and re-entered on p. 2). Finally, transport functions, e.g., Prandtl Number, Schmidt Number, 
thermal diffusivity, momentum diffusivity, etc, are defined (pp. 12-13) prior to definitions of the 
catalyst configurations, geometry, and flow conditions. 

Catalyst Configurations 

The physical parameters of several catalyst configurations (honeycomb monoliths and 
packed beds) and the appropriate flow conditions are described (p. 14) and used in the subsequent 
analyses. 

Combustion Monolith Performance 

The pressure drop values for (standard) 3.5 or 4.0 inch length 400 and 600 cells/in 2 
monoliths are determined as functions of flow rate (pp. 15-16, and see figure on p. 16). Heat and 
mass Nusselt functions are defined and used to calculate single-pass fuel conversion vs. monolith 
depth (p. 17-19). It is possible to determine laminar flow hydrothermal Nusselt functions for 
square or sinusoidal channels (which may be more representative of actual monolith channel 
shapes), but this was not included in the current analysis. Accordingly the Nusselt functions 
cannot be considered more accurate than about ± 30%. 

The results of the performance analysis for the monoliths are summarized in a table (p. 18) 
for the two cell densities, the two lengths, and 0, 20, and 40% by-pass air. The transport-limiting 
conversion of fuel components is plotted as a function of pressure drop for constant flow rate but 
extended monolith length (p. 19). High fuel conversions (« 99.5%) are possible with the 
600 cells/in 2 monolith at 1% pressure drop (but it requires 10 inches length). 

Post Combustion Converter 

Analysis of post-combustion downstream (1-atm vs. 3-atm upstream) conversion of CO 
and fuels and the pressure drop in 400 and 600 cpi monoliths shows excellent results for such a 
wide area monolith (pp. 20-25). After defining component diffusivity functions and describing 
monolith configuration (pp. 20-21), the pressure drop and conversions of CO and fuel vapors are 
calculated and plotted as functions of monolith depth (pp. 21-23). Finally, the results of an 
analysis of the overall conversions of individual hydrocarbons and CO including upstream 
combustor, homogeneous combustion, and the exhaust converter are tabulated (p. 23). 

Differentia] (Local) Nusselt Analysis 

A different analysis, where the local hydrothermal Nusselt functions are defined and 
integrated over the length of a monolith (pp. 24-25), are used to check the integral formulas. 



These values (p. 25), which are derived from different sources (Kays and Crawford, "Convection 
Heat and Mass Transfer," 3rd, 1993, vs. Rosner, 'Transport Processes in Chemically Reacting 
Systems", 1986) show higher conversions than the integrated values, e.g., 94% conversion for 
C 7 H l4 vs. 86% with the integral expression. This difference appears to be caused by the greater 
effect of hydrothermal entrance (as opposed to hydrodynamically developed flow with thermal 
entrance) and shows the integral formulas to be conservative by a factor of about two. 

Combustion with Packed Beds 

Pressure drop and fuel conversion functions are defined and used to analyze the 
performance of packed beds (pp. 26-3 1). Pressure drop is a strong function of the superficial fl ow 
rate, nearly second order (pressure drops are calculated and plotted against superficial linear flow 
velocity for various sizes of catalyst particles on p. 27). Gas to particle Nusselt numbers for 
packed beds are developed plotted for methane, isopentane, cycloheptane and nonene for three 
sizes of catalyst granules (pp. 27-28). Shape factors are not included in this analysis, so the 
results represent values for packed spheres with equivalent effect on AP, not with equivalent 
diameter. The conversion of fuel components (again assuming very fast catalytic rates) is 
determined (pp. 28-29) for the oblique flow geometry (which is defined on p. 16). 

For full perpendicular flow such as with an annular cylindrical packed bed rather than the 
conical bed , the pressure drop function is much greater for equivalent conversion of fuel 
components. This is caused by the higher power dependence of AP on flow velocities (form 
drag), while the conversion per unit depth of bed increases almost linearly with increasing flow 
rate. At flows higher than about 2 m/s the catalyst granules begin to act as bluff bodies which 
cause the flow to accelerate around them. While this increases backmixing which helps 
convection, it causes a greater effect on the pressure drop. Conversions for methane and 
cycloheptane are plotted as a function of pressure drop for beds of I -mm and 1/8-in. granules (pp. 
29-3 1). For oblique flow (2.5 m/s superficial velocity through the conical bed) with the small 
granules, CjK\4 conversions are calculated to be very high (> 99.99%) at 1% pressure drop, 
while for perpendicular flow (8.2 m/s) with the 1/8-in. pellets, C 7 Hi 4 conversions are calculated 
to be very low (= 62%) at 2% pressure drop. Clearly some radial configuration must be used to 
achieve high conversions with the packed beds without a pressure drop penalty. 

Required Catalytic Rates 

Finally, the specific catalytic rate constants required for conversions of methane and 
cycloheptane in a 600 cells/in 2 monolith are determined (pp. 32-35) for Capstone's conditions. 



The effect of gas diffusional restrictions within the macroporous (pores between agglomerated 
washcoat particles within a catalyst layer) and microporous (within the porous particles) is defined 
for CH4 and CyH^ (pp. 32-33). The effect of rate constant on surface temperature (for a point at 
the entry to the channel) is determined (p. 34) and shows the rates necessary to move fully into the 
transport limited regime, i.e., when the concentration of fuel near the surface approaches zero and 
temperature approaches a limit. This analysis could be turned around by calculating the conversion 
for a given catalytic rate constant. 

But, further analysis of the expected conversions requires measurements of the actual rate 
of the washcoated catalysts and catalyst beds. The analysis concludes with a complex rate 
expression for methane combustion by supported PdO that was based on measured rates (in the 
absence of transport limitations). This rate was based on PdO, an hexa-aluminate support 
(LaAl j |0|g) with a much lower surface area than the La-stabilized supported used in the Capstone 
monoliths and pellets. The interception of the curve showing the surface temperature as a function 
of catalytic rate constant (as on p. 35) and the curve showing the actual rate constant of the catalyst 
as a function of catalyst temperature (plotted on p. 37) given the actual performance of the catalyst. 




More detailed calculations would integrate the conversion through the monolith or packed 
bed and it would include the thermal conduction along the monolith channel or within the bed. 
This would represent a one-dimensional calculation (axial surface and gas phase temperatures) 
which is only an approximation to the actual two- or three-dimensional convective transport and 
reaction problem. Such calculations could provide estimates of extinction conditions for specific 
catalysts. Such calculations and further analyses are beyond the scope of the current comments, 
but could be performed by extending the current Mathcad 6 program. 

Conclusions 

Catalyst Performance 

Several conclusions can be drawn about the transport limiting performance of a catalysts 
used to combustion gasoline and methane fuels for Capstone's 24 kW microturbine: 

1 ) Commercially available small cell (600 cells/in 2 ) honeycomb monoliths 
should be seriously considered as an alternative to the PCI Microlith 
technology. A catalytic monolith configured as annular flow 600 cells/in 2 
honeycomb can achieve very high conversion of methane (99.99%) and 



modest conversion for gasoline (94% for cycloheptane) etc. with low 
pressure drop. 

2) If the fuel air mixing can be moved upstream sufficiently to allow use of a 
longer catalyst monoliths, significantly lower emission could result without 
unacceptable (> 1%) pressure drop ratio. For 1% AP/P drop (with greater 
monolith length), transport limited conversion of C7H14 should reach 99.5% 
in an active 600 cells/in 2 honeycomb monolith. 

3) Post combustion monolith catalysts should be considered for use in EZEV 
turbines, regardless of the method of fuel conversion during combustion. A 
short 400 cells/in 2 Alpha-IV metal monolith coated with a layer of supported 
Pt located at the inlet to the recuperator can give very significant emission 
reductions for CO and all hydrocarbons because of the lower linear velocities, 
high exhaust temperature, greater diffusivities at 1-atm. Conversions of CO 
and CH4 can be very high, while conversions with C7H14 and other gasoline 
fuels is less, it still can bring the gasoline range fuel components in the 
exhaust to a few ppm with the assumption that only a modest degree of 
homogeneous combustion occurs upstream of the turbine. 

4) Provided that significant homogeneous combustion occurs within the 
channels and downstream of the monolith (= 80% conversion giving 60 ppm 
C7H14 and 840 ppm CO leaving the turbine) and that a 400 cells/in 2 post 
combustion hydrocarbon oxidation catalyst was used downstream, gasoline 
emissions should fall below 0,5 ppm. 

5) The oblique flow/conical catalyst packed bed could also be considered as an 
option to the Microlith™ and honeycomb monolith catalysts. The 
performance of the packed beds can be very good in terms of fuel conversion 
with the smaller catalyst granules. Yet, these beds must have some element of 
radial flow as the pressure drop for a flow linear flow rate will cause 
excessive pressure drop. A word of caution here, packed beds are more 
susceptible to flash back caused by the higher degree of backmixing and 
regions where local residence times can be quite long relative to the screens 
and monoliths. 



Catalyst Durability, Robustness, and Costs 

1) ^MHBVHVHfc^ conc&ms were ex P resse ^ about the tolerance of 

combi^ high fiiel/air ratios and local 

unmixedness. While clearly high degrees of F/A control and uniformly 
premixed combustion gases are desired, the combustion catalysts should 
show robustness sufficient to survive short excursions of off -specification 
F/A ratio conditions as well as showing the durability necessary for long 
active life under normal operation at high temperatures. 

2) The greater the thermal mass of the catalyst, the greater its robustness in the 
face of locally or globally high F/A. In this regard the packed beds would be 



superior to the fine screens of the Microlith™ and the thin walled metal 
monoliths. 

3) A packed bed also has greater ability to disperse locally high F/A mixtures to a 
greater extent than the Microlith™ and to a much greater than the honeycomb 
channels. High thermal conductivity also improves robustness in locally high 
F/A (hot spot) conditions. Here the packed bed is perhaps superior to the 
metal monolith because of it greater radial conductivity, while the screens 
have very little ability (radiation) to disperse heat from locally hot zones. 

4) The Alpha-IV alloy honeycomb should be much more resistant to fast 
overheating to high temperatures than the PCI -microlith. The AIzeta/Grace 
electrically heated monolith previously examined by Capstone (Nomac) had 
entrance channels too large to give high conversion except with long length 
and smaller channels downstream, hence the large volume and weight. 

5) Packed beds represent the least cost catalyst for combustion because they use 
configurations and materials common with commercial VOC incineration 
catalysts. The metal monoliths have low costs because of their widespread 
acceptance and increasing use in automobile catalytic converters and as VOC 
catalysts. The premixed, gasoline-fueled, and recuperated nature of the 
combustion process for Capstone's microturbines is a good fit (with a slight 
push to higher pressure, higher throughput, low pressure drop, and greater 
robustness) with the commercial VOC catalysts available from several 
sources, including Prototech/United Catalysts. It would not appear necessary 
to develop new catalyst and catalytic engineered combustors for Capstone' 
microturbine combustion applications. 



The following table summarizes the findings of this analysis: 



Property 




Ranking 




Fuel Conversion 


Packed Beds > 


Screens > 


Monoliths 


Pressure Drop 


Monoliths > 


Screens > 


Packed Beds 


Durability 


Packed Beds « 


Monoliths ~ 


Screens 


Robustness 


Packed Beds > 


Monoliths > 


Screens 


Cost & Availability 


Packed Beds > 


Monoliths > 


Screens 
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Scanned copy is best available. 



Final Recommendations 

Use of Catalysts in Existing Engines 

The following recommendations are made on the basis of improving the existing 
microturbine catalytic combustion performance: 

1) Instaj[400 cells/in 2 1 -atm monolith with a Pt-based catalyst for post combustion 
CO an d hydrocarb oiumQYal^This unit should provide the greatest improvement 
in emissions wniTBttie effect on pressure drop at modest cost (« $200/unit 
installed). It could be used with all catalytic and even with non-catalytic burners 
to reduce CO and hydrocarbon emissions* However, it cannot be effective in 
reducing NOx emissions from non-catalytic burners. 

2) Use a 600 cells/in 2 Alpha-IV alloy metal monolith of 4-in. length with premixed 
fuel and air as a lower pressure drop alternative to the PQ Microlith™ system. 
Fuel conversion may be lower than the PCI unit, but with lower pressure drop; 
at equivalent AP/P, the monolith should show superior fuel conversion. A small 
amount of by-pass air could help with turn-down and could be used to improve 
emissions during normal operation by increasing residence time and allowing 
higher temperatures in the homogeneous zones within the channels and 
immediately downstream of the combustor. 

Use of Catalysts in Advanced Engines 

For advanced combustion systems, the fuel air mixing could be moved upstream of its 
current position to allow great premixing and to give greater volume for the catalyst. 

1) Nominally 1-mm to 1/16-in. packed beds in a conical configuration with about 
2-3 m/s linear superficial bed velocity could give good performance, especially 
with a small bypass air stream for post catalyst hydrocarbon reduction. An inconel 
screen could be located upstream to suppress flashback and ceramic foams or 
honeycomb monoliths could be used as mechanical supports to increase robustness. 

2) A longer 600 cells/in 2 Alpha-IV alloy metal monolith (from 3.5- to 8- or 10-in.) 
could also decrease hydrocarbon and CO emissions with acceptable AP/P by 
relocating the F/A mixing zone upstream. 



Additional Catalytic R&D Activities 

Additional Analysis Fuel combustion rates can be measured for specific catalysts 
under Capstone's combustion conditions. The rates can be incorporated in an extended analysis of 




p. 8 



t 



combustor performance for use in predicting fuel conversion, specifying required catalyst activity 
maintenance (e.g., Is it close to extinction at 1000°F inlet temperature?), and assisting in the design 
of current and next generation combustors (e.g., How much by-pass air is necessary? What are the 
limited on F/A ratios?, etc.). 

Advanced Materials 

I^Apaore robust units could be developed with ceramic structural components and state-of-the-art 
catalysts capable of resisting short (a few seconds) temperature excursions to 3000°F (1650°C). 
The attrition resistance and crush strength of the Norton S A 6576 catalyst pellets supplied to 
Capstone by SRI International also need improvement (their very large specific volume of 
macropores, while desirable, leads to poor mechanical strength), and other catalysts (e.g., 1.5-mm 
extruded, Pt-promoted, Mn-substituted hexa-aluminates) could be examined for activity and 
attrition resistance. 

the longer-terin.-^BHII^HHHIHIIHHHHHHHHHHHHHNB 

■■■■feSRI could coat steels with chromium and aluminum (in situ preparation of FeCrAl alloy 
films) to prepare inexpensive, high-performance recuperators and monoliths for use in advanced 
higher temperature engines. 

Systems — Additional engineering analysis and catalyst preparation and testing could 
provide Capstone with cost effective assistance in the design of EZEV turbogenerators with 
alternative fuels such as LPG, CNG, and diesel for hybrid vehicles. LPG could be used in the 
current designs with combustion and exhaust converter monoliths; CNG could be used with 
combustion and exhaust converter monoliths, but with by-pass air to increase F/A ratios and 
catalyst temperatures; and diesel fuel could require upstream mixing with long combustion 
monoliths or conical beds and longer exhaust converter monoliths. 






Gas Properties and Constants 
************************************ *************** * * *>**************************im«*+*********^ 



cm^O.OIm 



J=kg-m 2 sec 2 



N^kg-m-sec 2 Pa^Nm" 2 



W^J-sec 1 
atm = 101300 Pa 



Pref = 1-atm 
mn=0.000001-m 



R :=8.3 143— 
K 

********************** 



n Gas 



Rc=8.3143-J-mole l K l 



MolWt 



s(U) 



1 CH4 MW Q = 16.044- gm <* Q = 3.785- 10" 10 *m e Q = 148.6-K 

2 02 MW, -31.999-gm ^ = 3.464- 10" 10 -m e 1 = 106.7K 

3 N2 ^2 z28018 *g m c 2 = 3.798- 10" i0 m e 2 -71.4-K 



4 C02 MW 3 = 44.010-gm a 3 = 3.941- 10" 10 -m e 3 = 195.2K 



5H20 



MW = 18.015-gm c = 2.641- 10" I0 m s =809.1-K 



************** 



^ 6C5H12 

f 3 2-methyi-butane 



Tb5 = 301-K 
Tc5 



Tc5 = 460.4*K 



Pc5 =3.39- 10°- Pa Vc5 =306 



1.2593 



<r 5 = Vc5 3 -10"" I0 -m 



MW = 72.151-gm £ = 365.59994*K o =6.73866' 10 10 -m 



Tc7 = 572.2*K Tb7 ^374.1-K Pc7 .= 3.47-10 6 -Pa Vc7 =368 



p 7C7H14 

^ methyJ-cydohexane 



8C9H18 
1-nonene 



9 CO 



S 6 = 



Tc7 
1.2593 



<s 6 .= Vc7 3 .10" l0 -m 



MW 6 = 98.189-gm a g =7.1661* 10 _l ° 'm £ fi = 454.37942* K 
Tc9 '=592-K Tb9 =420-K Vc9 ^580 Pc9 = 2.34 10 6 Pa 

Tc9 1 



7 1.2593 



a ? ^Vc9 3 -10~ 10 -m 



MW ? = 126.243 gm a ? = 8.33955* 10 10 "m £ ? = 470. 10244 *K 



s c =91.7-K 



n 0 =3.6910" 10 m 



MW g = 28.01-gm 

**************************************T*****************^ 
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Define constants for transport integrals: 



qia 0 := 1.16145 



ojib Q - 0. 1 4874 aDa Q = 1 .06036 



aDb 0 =0.1561 



W x := 0.52487 := 0.7732 oDa 1 =0.193 cDbj =0.47635 



c*ia 2 =2.16178 o|* 2 =2.43787 00^=1.03587 oDb 2 := 1.52996 



Define transport integrals: 



aD^ =1.76474 aDb 3 =3.89411 



Qn(T):=cna 0 .exp[-(oMb 0 -lii(T))]+ J] a^-exp^^T) 

i|i=l 



QK3) = 1.0394 



QD(T) =aDa 0 -e?q)[- (aDb 0 -ln(T))] + J] oDa^ *cxp(-aDb^-T) OD(3) =0.95002 



Set Inlet Conditions: 

ych40 =0.01183 
y0c7hl4 =0.0021859 



P =3atm 

y0c5hl2 = 0.0028968 
y0c9h!8 =0.0016878 



These values for yOt have been determined via enthalpy 
functions on p. 1 1 for temperature rise from 1 000F to 1 51 1 F. 



Define gas composition w/r methane using moist (2%) air: 



y(yo,x) 



yO(l-x) 
(l-y0)-0.205-2ych40x 
(l~y0)<0.775 
yO-x "T 
(l-y0)-0.02+2y0-x 
0 
0 
0 



y(ych40,0.1) = 



0.0106 

0.2002 

0.7658 

0.0012 

0.0221 

0 

0 

0 



Zy(ych40,0.5) = I 



, Define function for calculating stoichiometric gas composition for combustion of 
2-methyi-butane using moist (2%) air: 



y5(y0,x) := 



(l-y0) 0.205- 8y0*x 
(l-y0)-0.775 
5(y0x) 
(l-yO)0.02+6yO-x 
yO(l-x) 
0 
0 



I + 2-yO-x 



y5(y0c5h!2,0.5) = 



Iy5(y0c51il2,0.5) = l 



Define function for calculating stoichiometric gas composition 
for combustion of methyl-cyctohexane using moist <2%) air: 



0 

0.192262 

0.770523 

0.007221 

0.02855 

0.001444 

0 

0 



-f y7(y0,x) 



0 

(l-y0)-0.205- 10.5y0x 
(l-y0)*0/775 
7-(y0-x) 
(I-yO)0.02 + 7*yOx 
0 

yO(l-x) 
0 



1 +■ 2.5-yO-x 
Iy7(y0c7hl4,0.5)-1 



^ Define function for calculating stoichiometric gas composition 
for combustion of 1-nonene using moist (2%) air: 



y7(y0c7hl4,0.5) = 



0 

0.19255 

0.771199 

0.00763 

0.027532 

0 

0.00109 
0 



0\ 



y9(y0,x) = 



0 

(l-y0)-0.205- 13.5-yOx 
(1 -y0)0.775 
9(y0x) 
(l-y0)-0.02+9*y0x 
0 
0 

yo-(i-x) 



1 + 3.5yOx 



y9(y0c9h!8,0.5) 



Zy9(y0c9hl8,0.5) = 1 



a****************************** ************^-*** ******* ************ 



0 

0.192692 

0.771413 

0.007573 

0,02748 

0 

0 



8.414148- 10 



~4 



-1 



i 



Define functions for evaluation of mixture properties (ideal gas) given fractional fuel conversion and temperature. 
Density and average molecular weight: 



pmix(x,T) 



E y( y ch40 ' x >kk' NW kk I ^ MWmix(x):= J] tfycMO.x^MW^ 



\ kk = 0 



kk = 0 



pmix( 1 , 300-K) = 3,470573*kgm 



MWmix(l) =28.485 -gm 



Define diffusion coefficients for specific binary diffusion pairs: 



Difl(n,ra,T) -1.8829-10 



r27 



MW MW 

V a a 



\2 

a -i-o \ 
5 -QD 



in 8 5 4 „ - 1 -r- 1.5 

— 10 *gm -cm -sec K 



0.5 



P 

Pref 



Difi(0 t 2,300K)=0.07417-sec 1 *cm 2 



Difl(2,7,300-K)=0.01597-sec 1 -cm 2 



Define functions for use with various fuel air mixtures: 
l-y(ych40,x) 0 



S Dch4{x,T) : = 



Dc5hl2(x,T) := 



' y(ych40,x) jgas 

. 2-* Difi(0,jgas,T) 
Jgas= 1 

l-y5{yOc5hl2,x) n 



s 



ySlyOcShU^ 
Difl(5jgas,T) 



Dc7M4(x,T) 



jgas= 1 

1 - y7(y0c7hl4,x) n 



y7(y0o7hl4,x) jg| 
Difl(6,jgas,T) 



Dc9hl8(x,T) 



jgas= I 

1 - y9{y0c9hl8,x) n 



£ 



jgas- 1 



y9(y0c9hl8,x). ga 
Diff(7Jgas,T) 



Dch4{0.5,298*K) =7.31327-10 6 'm 2 -sec 1 



Dc5hl2(0.5,298 K) =2.28466*10 6 *m 2 'sec 1 



Dc7hl4(0.5,298 K) = 1.94254*10 6 *m 2 *sec 1 



Dc9hl8(0.5,298 K) = 1.53123*10 6 -m 2 'sec 1 



V^y Display gas mixture properties 



idsp -0..9 



Tdsp.^ =(750+ 5<Mdsp)K 



Tin =|(I000- 32)-^ + 273.15 



Tin = 810.92778-K 



Tout 



(1511 - 32)~ + 273.15 



K Tout-1094.81667-K 



Plots (below) of Di for CH4, C5H12, etc. show large gap between CH4 and the other hydrocarbons 
and show very little effect of fuel conversion (even for methane) on Di because the fuels are so dilute. 

Therefore we can choose methyl cycfohexane as the molecule representative of naphtha and ignore 
the effect of conversion on gas properties 



Dc9hl8(.5,Tdsp id|p ) 








7.5*10 2 K 




3.61 119- 10" 5 -ra 2 sec" 


8 l& K 




4.02644- 10" 5 -m 2 sec 


8.5 lO^K 




4.45884 10 5 m 2 sec" 






4.90801-10 "W-sec" 


9.5- lO^K 




5.37357-10" 5 -m 2 -sec" 


MO 3 -*: 




5.85522- 10" 5 - m 2 -sec- 


1.05- lO^K 




6.35265•10- 5 -m 2 -sec" 


l.MO^K 




6M56- 10" 5 m 2 sec" 3 


1.15- IO^K 




7.39384* 10" 5 -m 2 sec" 1 






7.93712- 10" W-sec" 1 



1.001 



Dch4 ('- Td3 P idJ 



1.0005 




0.9995 - 



0.999 



Tds pid» 


Dc7hl4(.5,Tdsp idap ) 


7.5- lO^K 




1. 03475- 10" W-sec -1 


SIO^K 




l^rMOWsec" 1 


8.5 K^K 




1.28482- 10" 5 -m 2 -sec -1 


MO^K 




1.41758- 10" W^ec -1 


9.5- lO 2 ^ 




1.55533- 10" 5 m 2 -sec 1 


1 10 3 K 




1.69798- 10" 5 m 2 sec 1 


1.05- lO^K 




1.84541- 10" 5 -m 2 -setf 1 


l.MO^K 




1.99755- 10" 5 m 2 -sec l 


1.15- lO^K 




2.15429- 10" 5 -m 2 -sec l 


I^IO'-K 




2.31555- lOW-sec" 1 



^ Define interpolation formulas for Dch4 and Dc7h14 vs. T(K) at x=0.5 via Jandel Scienttic's Tablecurve program 

Dch4(.5, 1000K) = 0.585522' 
Dch4T( 1000K) = 0.585523-sec" 1 -cm 2 



Dch4T(T) :=«p(- 1 1.718527 + 1.6257167h(l) - DoM(-5f 1000 . K) =0 .585522-sec" 1 W 



Dc7hl4T(T) =exp(- 12.923061 + 1.6260558- In g) - 82 - 48 ^ 98 K j Dc7hl4(.5,1000 K) =0.169798-sec" 1 -cm 

Dc7hl4T( 1000-K) = 0. 169798' sec" 1 'cm 2 



************ * 



Define viscosity functions for major constituents: 



2.6693*10 



0.5 



cm / giti 



sec 



U 



(MW f T) 
(MW 2 -T) ai 
(MW 3 -T) 0 ' 5 
1.116.(MW 4 T) 



2.6693-10 



r26 



r26 



2.6693- 10 



2.6693- 10" 26 - 



10 5 *cmi 



10 5 - 



K 



cm | f gm 
K 



0.5 



0.5 



lO^cm 



sec 



f) 



(3.115- 10- 8 cm) 2 W 

V 



514 K 



2.6693- 10" 26 - 



10 5 . 



sec 



cm i ( ' gm 



0.5 



K300-K) = 



1.10281*10 J 
2.0603* 10~ 5 
176981 -10" 3 



1.51851-10 
1. 06784*10 



rs 



kg'm 'see 



fimix(x,T):= ^ 



y(ych40,x) -n(T) 



n-0 



J] y<ych40,x) a 



m = 0 



i + 



/«T) \°- 5 /MW ^i 0 " 25 



/ MW \ 



8- 1 + - 



MW 



0.5 



jimix(0.5,300K)=0.00002'kg*m l 'sec 1 



y(ych40,0.5) = 



0.00592 

0.19074 

0.76583 

0.00592 

0.03159 

0 

0 

0 



ych40 =0.01183 



5*10 



4.5*10 



"5 _ 



^f(°' Tds Pi<bp) 



3.5* 10 



3*10 



r* - 




700 800 900 1000 1100 1200 



The viscosity of the CH4 - air mixtures vary 
little with composition (< 0.5%), again 
because the fuel is so dilute. We can 
therefore neglect the effect of composition by 
chosing x=O.S in all cases. 



1.004 



1.002 - 



^(• 5 ' Tda Ptdsp) 




0.998 - 



0.996 









7.5-10 2 K 




3.38338 10~ 5 -kgm" I -sec~ 1 


sio^k 




3.53056 10' 5 kgm" 


^sec" 1 


8.5- lO^K 




3.67429- 10" 5 kgm~ 


•sec" 


910 2 -K 




3.81491 -10" 5 -kg-m~ 


-sec 


9.5- lO^K 




3.95267- 10" 5 -kgm" 


l -sec" 1 


MO^K 




4.08779- 10" 5 kg nf ] 


-sec" 


1.05- 10*K 




4.22047 10" 5 kgm~ ] 


-sec" 


l.MO^K 




4.35088 10" 5 kg m ] 


-sec" 1 






4.47917- 10" 5 kgm" ] 


•sec" 1 


1.2-lO^K 




4.60547-10" 5 -kgm~ ] 


-sec" 1 



jimixT(T) :=exp[-0.70216045 + 0.640874-lnffl +- " ^ ^^ ^ j-lO ^kg-m ^sec" 1 



mnix(.5,1000 K) =4.08779-10 5 kgm ^sec 1 



jimixT(1000 K) =4.08782-10 5 -kg*m ^sec 1 



Define thermal conductivity functions for light molecules 



igas =0,1. .4 



aX: 



-1.869 




0.08727 " 


-.3273 




0.09966 


0.3919 


U = 


0.09816 


-7.215 




0.08015 


. 7.341 . 




.0.01013. 



cX = 



1.179-10* 

- (3.743- 10" 5 ) 

- (5.067- 10* 5 ) 
5.477- 10" 6 
1.801 10 -4 







CH4 




(3.614- 10" 8 ) 
9.732- 10" 9 


02 
N2 


<ft = 


1.504 10" 8 
-( 1.053- 10~ 8 ) 


C02 




-(9.1-10"') 


H20 



*(i,T) 



al. + bX,- + c X.- - +*.- - 
' 1 K '\K MK, 



0.00001- 



w 

cm-K 



^ Xmix(x,T) '— ^ 
i = 0 



y(ych40,x).X(i,T) 



Y, y(ych40,x).- 
j = 0 




Higas.lOOOK) 



W/cm/K 



0.16716 kg msec 




0.07163-kg-m-sec 


'-KT 1 


0.06292-kg-m-sec 


3.KT 1 


0.06788-kg-msec" 


3 -KT l 


0.0863 lkg-m-sec" 


3 -KT l 



Xmix(.0, 1000-K) =0.00066 * 



W 
cm-K 



0.075 



0.07 - 



0.065 - 



0.06 



0.055 



0.05 




1200 



Tds p^ 


lmix(0.5,Tdsp idsp ) 


7.510 2 K 




0.05362 kg- m-sec ^KT 1 






0.05621 kg-msec 3 ^ 1 


8.5- lO 2 * 




0.05872 kg m-sec 3 K { 


^lO 2 ^ 




0.061 17-kgmsec^K" 1 


9.5- lO^K 




OX^SSkg-msecf 3 -*: 1 


l-lO^K 




0.06589 kg m sec -3 KT l 


1.0510*-K 




0.06818kgmsec 3 *KT l 


Ll-IO^K 




0.07043 kg* m sec" 3 K: ! 


1.15 10*-K 




0.07266-kg*m sec 3 -KT 1 


1.2 10 s K 




O^SSkgmsec^KT 1 



1.01 



^nix(0 t Tdsp i<fap ) 
^(- 5 ' Td *idsp) 



1.005 - 



Xnux^i,T< 




0.995 " 



0.99 



H00 1200 



The thermal conductivity varies > 1% with 
conversion. Selection of 50% conversion 
gives and error range of +.5-.5% at 1000K. 
This is likely to be smaller than the accuracy 
of measured thermal conductivity. Therefore 
we again define a composition independent 
formula to increase speed in calculating 
derived transport properties 



XmixT(T) := exp|- 7.073383 + 0.6397005 -ln^j + ~ 65 ' 2 ^ 87>K ) .W-nT l -KT l 



Xmix(.5,1000*K) = 0.06589'm l 'K 1 *W 



A,mixT(1000-K) = 0.06589'm K K 1 *W 



Determine the heat capacity and enthalpy functions of the fuel-air mixtures for the evaluation of temperature rise 
Define heat capacities and heat of combustion: 



Cpcoef - 



19.25 .05213 

28.11 -0.000368 

31.15 -0.01357 

19.8 0.07344 

32.24 0.001924 
-9.52 0.5066 
-76.19 0.7867 
-37.18 0.8122 



0.00001197 
0.00001746 
0.0000268 
-0.00005602 
0.00001055 
-0.0002729 
-0.0004204 
-0.0004509 



-0.00000001132 
-0.00000001065 
-0.00000001168 
0.00000001715 
-0.000000003596 
0.00000005723 
0.00000007516 
0.00000009705 



(W*..T):= £ Cpcoef^- (| 1* ± 

jjj = 0 



J/moi/K 



Cp(0,300K)=35.66066-K 1 \T 
Cp(6,300*K) = 124.01332-K~ 1 • 



Cp(5,300K) = 119.44421"K 1 *J 



Cp{7 ? 300-K) = 168.51935"K *J 



Cpmix(x,T) := ^] CpGk,T)7(ych40,x) jk 
jk-0 



y(ych40,0.5) = 



Cpmix(0.5,998 K) =33.66993 \F-KT 



0.00592 

0.19074 

0.76583 

0.00592 

0.03159 

0 

0 

0 



Define heats of combustion 



ACpcomb(T) :=(Cp(3,T) + 2-Cp(4,T))- (Cp(0,T) + 2-Cp(l,T)) 
ACpcomb5(T) = (5Cp(3,T) + 6-Cp(4,T)) - (Cp(5,T)-h 8-Cp(l t T» 
ACpcomb7(T) -(7 *Cp(3,T) + 7Cp(4,T))- (Cp(6,T)+ 10.5Cp(l,T)) 
ACpcomb9(T) =(9-Cp(3 > T) + 9 Cp(4,T)) - (Cp(7,T) + 13.5 Cp(l,T)) 



CH4 + 202 => C02 + 2H20 
CSH12 + 802 «=> 5C02 + 6H20 
C7H14 + 10.502 ==> 7C02 + 7H20 
C9H18* 13.502 «> 9C02 + 9H20 



AHf := 



-74900- J 
OJ 
0-J 

- 393800* J 
-242000 J 

- 154600- J 
-119400- J 

- 103600-J 



ACpcomb(298K) = 10.38812'K "J 



ACpcomb5(298 K) £anAAa ^ t 2. ~2 tV r 
— ^ - = 6.97449*kg*m 'sec *K 



ACpco m b7(298-K) =9>35823 . K - 1 mJ 



ACpcomb9(298-K) _ fi3177 . K -i ^ 
9 



\ Define enthalpy function for individual components and AH for fuels: 

•T 

Cp(nn,Tx) dTx 



:Hf(mi,T):=AHf + 



•»298K 



UJ AHcomb(T) := [ (AH^ + 2-AHfJ - (AHf Q 4- 2- AHf J ] + 



ACpcomb(Tx) dTx 



298-K 

rT 



Jjt AHcomb5(T) :=[[(5-AHf) 3 + 6-AHfJ - (AHf Q -h 8-AHf I )] + 



298•K 
T 



H AHcomb7(T) :=[[(7-AHf) 3 + 7-AHfJ - (AHf 0 + 10.5-AHfJ] -f 
AHcomb9(T) :=[[(9-«K) 3 +9-AHfJ- (AHf 0 + 13.5-AHfJ]* 



298-K 

rT 
298-K 



ACpcomb5(Tx)dTx 
ACpcomb7(Tx) dTx 
ACpcomb9(Tx) dTx 



Hf(igas, 1000-K) 



36559.44819 kg m 2 sec~ 2 



22590.0387 1 kg- m 2 - sec" 2 



21484.78794 kg-m 2 -sec 2 



- 3.60367- 10 5 kg-m 2 sec 2 
-2.15959 10 5 kg m 2 -sec 2 



Evaluate a few values of AH for CH4, C5H12, C7H14, and C9H18; 

AHcomb(800-K)-«.00306-10 5 'J AH ^ (S0 °' K) - -6.69221 ■ 10 5 -J 



AHcomb7(800K) = ^ 2472M() 5 



AHcomb9(800K) = ^ 260()2 . I0 5 



******* ************** ******** ******* * *************************************** ********************************** *************************** 

Define the mixture enthalpies and determine the fuel composition needed for 51 1F temperature rise from 1000F. 



igas=0,l..7 



Hmix(x,T) := J] y(ych40,x)^-Hf(kgas,T) 



ych40 =0.01183 



kgas = 0 



Hmix(0,810.93K) = 1.006514* 10 4 *J Hmix( 1 , 1094.82K) = 1.006518* 10 4 *J 



Determines CH4 inlet 
concentration 



Hmix5(x,T) := y5(y0c5hl2,x)^-Hf(kgas,T) 
kgas = 0 



y0c5hl2 = 0.0028968 



Hmix5(0, 810.93 K) =1.058499' 10 4 'J Hmix5( 1 , 1094.82 K) = 1.058492" 10 4 *J 



Determines C5H12 inlet 
concentration 



|f Hmix7(x,T) := ^ y 7(y0c7hl4,x)^Hf(kgas,T) y0c71il4 =0.0021859 

;""= = ty* 5 = 0 

O Hmix7(0,810.93 K) = 1.075807" 10 4 -J Hmix7( 1, 1094.82 -K) = 1.075812" 10 4 *J 



Determines C7H14 inlet 
concentration 



Hmix9(x,T) := J] y9(y0c9hl8,x)^ffi(kgas,T) 
kgas = 0 



y0c9h!8 =0.0016878 



Determines C9H18 inlet 
concentration 



Hmix9(0,810.93-K) = 1.083024" 10 4 -J Hmix9( 1 , 1094.82-K) = 1.083028* 10 4 *J 



***********************************^***************^ 



********************* * ****1Mk **********+** 



Define momentum and thermal difushrities: 



«<igas,T) 



A(igas 5 T)RT 
Cp(igas,T).P 



P = 3.039' Hf 'Pa 



<x(0,300>K)=0.07813"sec 1 'cm 2 



amix(x,T) : 



,_ Xmix(x,T) R T 



Cpmix(x,T) P 



omix(.5 7 1000K) =0.5352-sec 1 -cm 2 



omix(x,T) := 



mnix(x>T) 
proix(x,T) 



pmix(.5,1000-K) = 1.04117-10 * *gm cm 



omix{.5 > 1000K)=0J926-sec 1 W 



Dch4(.5J0()O*K)=0.5855"sec 1 "cm 2 
Dc7hl4{0.5 ? 1000K)= 0.1698" sec" 1 "cm 2 



Dc5fal2(.5,1000K) -0.1957*sec 1 "cm 2 
Dc9hI8(.5J000*K)=0.13431"sec~ 1 'cm 2 



Define Prandl, Schmidt, and Lewis numbers 



Pr(x,T) = 



umix(x,T) 
amix(x,T) 



Pr( . 5, 1000K) =0.73359 



Sc(x,T) = 0mix < x ' T > Sc{J, 1000K) -0.67054 Sc5(x,T) - umix ^ x ' T > Sc5(,5, 1000 K) -2.00639 
E>ch4(x ) T) Dc5hI2(x,T) 

Sc7(x,T) - umix < x ' T > Sc7(.5, 1000 K) =2.31224 Sc9(x,T) := < * ng { x » T > Sc9(.5, 1000-K) =2.92325 
Dc7hl4(x,T) Dc9hl8(x,T) 



2.5 



K 5 > Tds Pi<fap) 

Sc5(.5,Td^ 



0.5 



1 - 



J_ 



700 



800 



900 1000 



1100 1200 



These relations have weak 
temperature dependence and 
(presumably) very weak 
dependence on composition 



LC(X,T) amixfx T) 1°°°^) = 1-09404 Le5(x,T) := 00511 12 < X > T > Le5(.5,1000*K) =036563 

' amix(x,T) 

Le7<x,T) ^ Pc7hl4(X>T) Le7(.5 T l000K) =0.31726 Le9(x,T) • = Pc9hl8 C x » T ) Le9(.5,10O0K) -0.25095 
amix(x,T) amix(x,T) 



1.2 



Le7(.5,Td^ idap } 
^(- 5 > T ^id3p) 0 .6 



0.4 



0.2 



700 



800 



900 1000 
Tds Pidsp 



1100 



1200 



The Lewis numbers show 
approximately the temperature 
rise at the waf! of an active 
catalyst relative to the adiabatic 
temperature rise at that point 



Clearly temperature (and composition) 
have little influence on the AT at the 
wall, while the nature of the fuel 
molecule has a strong effect (from 
0.25*ATadiab for C9 to 1 ,1*ATadiab 
fbrCH4) 



Can remove compositional dependence for most property fundctions 



pmixT(T):-MWmix(0.5)^ »mixT(T) = MmixT(T) 

pmixT(T) 

ScT(T) := MmixT(T) . = MxTffi 

Dch4T(T)pmix(0.5 , T) Dc7hl4T(T> pmix(0.5 , T) 



These functions speed calculations of Nu numbers 
and heat and mass transfer coefficients 



This ends the extended section that defines gas properties for the combustion of several fuel air mixtures. 



******************************** ******** ******************************************************************* ** 

Define honeycomb configurations 

Diao =9 in Diao= 22.86 -cm Diai =4 in Diai = 10. 16 "cm Lhc=3.5in Lhc = 8.89 'cm 

With bypass change Diao to 8.5 in Bypass = 0 Bypass - 0, 0.2, 0.4 Lhc - 3.5, 4 in 

Honeycombs: Represent as square channel metal alloy wall 0.05 mm, uncoated 

CelM-lOOm 2 dwalll := 0.009 cm Ceils2 - 200 m 2 dwall2 - 0.009- cm Describe 1 00 and 200 cpi 

dcelll -(Cells!)" 05 - dwaill dccll2 =(Cells2) 0 5 - dwa!12 dimensions, although not 



dceill =0.245 'cm ksl -0.002 cm dcel!2 = 0.17061 -cm ks2 -0.002- cm 

ecelll -dcelll 2 Cdlsl ecelll =0.93039 ece!12 := dce!12 2 Celis2 «cell2 = 0.90229 



expected to be used 



Cells4 -400 m 2 dwaII4 - 0.007 cm Cells6 -600 in 2 dwall6 - 0.007 cm 

dcelM -(Cells4) ° 5 - dwa!14 dcell4 =0.1 2 'cm dcell6 -(Cells6)'° 5 - dwa!16 dcel!6 = 0.0967 'cm 

dceltt =0.12'cm ks4 -0.001 cm dcelI6 = 0.0967 *cm ks6 -0.001 cm 

ecelM -dcell4 2 Cells4 £celI4 = 0.8928 £ce!16 -dcell6 2 -Cells6 sce!16 =0.86955 



Packed Beds: Cylindrical (conical) geometry 



Dinr -6 jn .__ 
Dinr = 15.24 *cm 



Doutr -7 in, 
Doutr = 17.78 'cm 



Wbed-7-in 



Lbed -0.75 in 



Wbed = 17.78 *cm 



Smono = ^-(Diao 2 - Diai 2 ) 



Sbed:=> D8nr ^ PoBlr -Wbod 



Smono 
Sbed 



-0.35714 



*************** **************** 



Define initial combustion conditions 



********* ******* ******************** 



Adjust U0 to fit U(Tin)=7 m/sec w/o bypass flow 



Uref-7 



m 
sec 



Tin -(537.778 -h 273.15)K Tin = 810.928'K 



P= 3.039- 10 5 'kg-m ^kxT 2 



M0=Urefpmix(0,Tin> 



UQ . MO (1- Bypass) 



7'((9-m) 2 -(4-m) 2 ) 
4 

1 



M0 =0.29601'kg'sec 



Tout- (821.667 -i-273.15)-K 
Tout = 1094.817-K 



U(T) -U0-— 
Tin 

Ubed(T) - 



M0 



pmix(0 ? Tin)-Sbed 



U(Tin) =7*m-sec 



Ubed(Tin) =2.5-m*sec 



GO =U(Tm)pmix(0,Tm) 



GO = 8.987 49* kg* m^ 2 * sec 1 



********#******#***#»*****#************« 

Evaluate Reynolds numbers for use in calculating the momentum, heat, and mass transfer functions: 



Rehc4(T) := 



GOdceM 
8cell4-MmixT(T) 



Refac4(Tm) -339.11451 



Rehc6(T) 



ecell6maixT(T) 



Rdw6( Tin) =280.56407 



For laminar flow with entrance effect (adapted from Kays & Crawford, "Convective Heat and Mass Transfer"): 



cfhc4(T,z) - 



14.227 
Reho4(T) 



1 + 0.046263 |Rchc4{T) ^i 



0.45363 



cfric6(T,z) 



14.227 
Refcc6(T) 



dce!16 



1 + 0.046263 *^Reiic6(T) 



0.45363 



Average coefficient of friction 




The entrance effect for AP of a channel in laminar 
flow is significant even to 3.5 in. length. 



cficx(x,z) =^^(l+0.0462363( X -z)) 0 - 45363 



. dce!16 Uds Pjd^, 



Rex4. 



_ dce!14 Uds Pjd ap 



,d * tcelK «mixT(Tavg) * w "kisp ^j^^ 

dceU4\ pmixT(T«vg) ( Uds PjdJ 2 Lbc 



APx4.. =4-cficx Rex4.. , 

,ds P 1 •** Lhc 



ecel!4 / <fcell4 



0.006 



0.005 



0.004 — 



APx6, 



0.003 - 



P 



0.002 ~ 



0.001 




AP/P < 0.5% to about 1 1 m/sec 



pmixT(T> 



APhc6(T) ^cfiiceCT^Lhc)- 



f U(T) \ 2 
ieccll6/ Lhc 



dceU6 



a*********************************** ********* 



Tavg = 952.8725*K U(Tavg) = 8.22528 nr sec 1 



^OM =0.00271 
P 



APfac6(Tout) =at00441 



APfac6(T g vg) ==000352 
P 



pmixT(T> 



APhc4(T) -4-cfric4(T,Lhc)- 



U(T)\ 2 
ecel!4/ Lhc 



dceU4 



APhc4(Tavg) 



-0.00228 



r***************** ******** ************************^ 



Find 1-D heat and mass transfer rates via Nussett formulas 



Averaged Nu for 600 cpi square channels (fully developed 
laminar flow) taken from Rosner: 



7.6-(2cm) 



NuhavgR6(z) = 2.976- 



1 + 



7.6-z 



Rehc6( Tavg) -Pr( .5 , Tavg) -dceU6 



Rehc6(Tavg)Pr( .5 , Tavg) dceU6 



Nnhavg^6(IJic) =2.98649 



-0.85369 



NumavgR6(z) -2.976- 



1 + 



7.6-z 



Rehc6(Tavg) ScT(Tavg) doell6 



NumavgR6(Lhc) =2.98441 



i Num5avgR6(z) = 2.976- 



1 

81* 



1+1 



7.6z 



\Rehc6(Tavg)-Sc5(.5,Tavg)dcell6 



Num5avgR6(Lhc) =3.11102 



' Num7avgR6(z) =2.976 



1 + 



/ — - 

\Rehc6(Tavg).S 



6z 



Sc7T(Tavg)dcell6 



Num7avgR6(Lhc) =3.16167 



Num7avgR4(z) - 2.976- 



i 

81* 



1 + 



/ 7.6z 

\Rehc4(Tavg)'Sc7T(Tai 



ivg)dce!14 



Num7avgR4(Lhc) =3.3902 



Num9avgR6(z) = 2.976* 



1+1 



7.6*z 



\Rehc6(Tavg)- Sc9( .5, Tavg)*dcell6/ 



Num9avgR6(Lho) =3.2775 



The Nusselt numbers for mass transfer show small (< 5%) compensation for the low diffusivity of the higher hydrocarbons. 

***************** ********************* ************************************************** **^^ 



Use exponential functions implied by the definition of average (integrated) Num functions to calculate exit mote fractions. 



yc7hl4(T,z) :y0o7hl4-exp( -**W**V> \ 
\Rdic6(T)-So7T(T)doeU6/ 



yc7hl4(Tavg,Lhc) =0 13705 
y0c7h!4 



ych4(T,z) -ych40exp^ 



I -4zNnmaygR6(z) 
Rehc6{T)ScT(T)dcell6 



ycl^Tavg.Lhc) =Qm533 
ych40 



yc5h!2(T,z) =y0c5hl2exp 



/ -4z-Nmn7«vgR6(z) \ 
\Rdic6{T)Sc5(0.5,T)-dceB6j 



yc5h!2(T«vg > Lho) _ 0101H 
y0c5hl2 



J F \Rfihc6{T)Sc9(0.5,T)*dceU6/ 



yc9h!8(Tavg,Lhc) 
y0c9hI8 



= 0.19604 



o.oi 



yc7hl4(T*vg t z idj9i ) 

,0.001 
yc3hl8(Tavg,z idj!p ) 



no 



0.02 



0.04 



0.06 



idjp 




0.08 



0.1 



z 9 =0.095'm 



z 9 =3.74016 'in 



Summary If 600 cpi Metal Monolith Performance 

%Bypass Length(in) Diao(in) Yout/Yin C7H14 YoutfYinCH4 AP/P 



0 


3.5 


9.0 


0.13705 


0.001533 


0.00352 


0 


4 


9.0 


0.10679 


0.000611 


0.00400 


20 


3.5 


8.5 


0.11888 


0.000906 


0.00324 


20 


4 


8.5 


0.09055 


0.000335 


0.00368 


40 


3.5 


8.5 


0.06266 


0.000089 


0.00240 


40 


4 


8.5 


0.04317 


0.000023 


0.00273 



Bypass = 0 
Lfac=3.5 -in 
Diao = 9 in 

These values are conservative 
and do not account for the 
higher average temperature 
with the higher F/A ratios and 
do not account for 
hydrothermal entry. 



A********************** 
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Conversions at AP/P=0.01 
CH4 = 1.06*10 A -8 




—in 




\ 
* 






C7H14=4.8*10 A -3 



0.01 0.02 0.03 0.04 

PLhc 



2=2525 cm 



zx r 25.25*cm 



APho6(Tavg)zx_ 
P-Lhc 



0.01000 yc7hl4(Tavg,zx) 
y0c7h!4 



Conversions at AP/P=0.02 
CH4-1.14M0M6 
C7H14=2.41*10 A -5 
C9H18=2.21*10M 
z=50.5 cm 



ych4(Tavg >ZX ) = ^ yc9h!8(Tavg,z X ) = lMm . 1(f 



ych40 



y0c9hl8 



The 600 cpi monolith has low pressure drop at 8.89 cm length. If if s length were extended to give 1 .0% and 2.0% 
fractional pressure drop, then C7H14 conversions would increase to 99.5% and 99,998% respectively 



Set U P properties for CO to calculate the conversion of dilute CO and hydrocarbons downstream of the turbine 



Dcol(x,T):=- 



P 
Pref 



\^ y(yoh40,x) jgas 



jgas = 0 



Diff(8jgas,T) 



Dcol(.5,Tin) ==U2854'sec 1 "cm 2 Pref = l'atm 
Dch4{.5,Tiii)=0.41195'sec n 'cm 2 P - 3 atm 



Dco7(x,T) : = - 



P 
Pref 



jgas=0 



Diff(8,jgas,T) 



Dco7(.5,Tm) = 1.12396'sec 1 'cm 2 Pref = Tata 
Dc7hl4(.5,Tin) = 0.1 1846* sec" 1 'cm 2 P=3*atm 



r=» Dch4(.5 



0.0003 



.0002 



Dcol(.5,T<fap i(J ^) 
Dco7(.5,Tdv^.0001 




600 800 1000 1200 



TdS Pid^ 

750K 



800K 



850K 



900-K 



950-K 



1000.K 



1050K 



HOOK 



1150K 



1200K 



9.90331 



10' 5 -m 2 -sec 



1.10324 10 4 m 2 -sec 



1.2208- 10" 4 m 2 -sec _ 



L34291 10" 4 m 2 sec 



1.46949-10 ~ 4 -m 2 sec 



1.60044-10- 4 -m 2 -sec' 



1.73572 



i0~ 4 -m 2 -sec" 



1.87523- 10" W-sec" 



2.01893 



10~ 4 -m 2 *sec~ 



2.16675- 10" 4 m 2 -sec 



DcoT(T) -exp 



- 10.80086 -h 1.6357427 m 



- 28.157906 - 



10" 4 m 2 - sec 1 Dcol(.5,Tavg) = 1.47689 sec 1 "cm 2 

DcoT(Tavg)^ 1.47689* sec 1 'cm 2 



dreco ~ 17 - in 
dreci -6 in 



dreco -43.18 'cm 
dreci - 15.24 *cm 



Lrec -3.5 in 
Lrec =8.89 'cm 



M0 =0.29601 -kg" 
M0 



GrecO 



Arecup 



Prec -l atm 



Arecup 



dreco 2 - dreci 2 



GrecO =2.30904-kg*m~ 2 ' sec 1 



pmixrec(T) 



MWmix(0.5) Prec 
RT 



Arecup = 0.1 282* m 



GO = 8.98749 -kg-m" 2 - sec 1 



pmixrec(T«vg) = 0.36422*kg*m 



Urec(T) :=- 



GrecO 
pmixrec(T) 



Urec(Tavg)=6.33964-m"sec 1 



SccoT(T) 



mnixT(T) 



DcoT(T)-pmixrec(T) 



Rerec4(T) := 



GrecOdceM 
MxT(T) 



Rerec4(Tavg) =69.96214 



Schmidt numbers are 
independent of pressure 



Rerec6(T) := 



GrecO- dce!16 
HmixT(T) 



Rerec6(Tavg) =5637495 



Tavg = 952.8725' K 



cfricrec4(T,z) = 



14.227 
Rerec4(T) 



1+0.046263- Rerec4(T) 



dce!14 



0.45363 



cfiicrec6(T,z) 



, 14.227 
Rercc6(T) 



1 + 0.046263- Rerec6(T) 



0.45363 



Calculate monolith pressure drop in several monoliths with entrance effect: 



k 4Prec4(T) ^^cfricreoK^Lrec) 



\ ecell4 / Lrec 



dcelM 



APreo4(Tavg) =0.00557 
Prec 



IB APrec6(T) =4 cfricrec6(T,Lrec) 



^T)JM)| 2 



\ ecel!6 / Lrec 



dcei!6 



APrec6(Tav g) _ 

Prec 



0.00898 



Numavgrec6(z) = 2.976- 



1 + 



7.6*z 



Rerec6( Tavg) - ScT( Tavg) dcelI6 



Numavgrec6(Lrec) =2.97616 



Numavgrec4(z) -2.976 



1 + 



7.6-z 



Rerec4(Tavg) ScT(Tavg) dcell4 



Numavgrec4(Lrec) =2.97649 



Numavgrec7R6(z) := 2.976- 



l 

*1« 



1 + 



7.6-z 



Rerec6(Tavg) Sc7T(Tavg) dce!16 



Numavgrec7R6(Lrec) =2.98024 



Nmnavgrec7R4(z) =2.976 



I 

8]* 



7.6z 



Raec4(Tavg)Sc7T(Tavg)-dceU4 



Numavgrec7R6(Lrec) =2.98024 



NwnavgreccoR6(z) -2.976 



i 

8l« 



/ 



7.6z 



Rerec6(Tavg)- SccoT(Tavg)dcell6 



NumavgR6(Liec> =2.98441 



NumavgreccoR4( z ) =2.976 



l 

«1 8 



1 + 1 



7.6z 



\Rereo4(Tavg)SccoT(Tavg)dcell4 



NumavgreccoR4(Lrec) =2.97663 



The average (integral) Nusselt numbers at 1-atm with the larger flwo cross section show very little contribution of the therma 
entry length, i.e. the flow velocities and heat and mass transfer coefficients have values close to those in long channels. 

Now determine CH4, C7H14, and CO conversions in 400 and 600 cpi honeycombs, assuming 80% C7H14 conversion 
downstream of the combustion monolith and upstream of the turbine 



Ly ych4rec6(T,z) :: ych4(Tavg > Lhc) exp 



- 4z-Numavgrec6(z) 
Rfirec6(T)-ScT(T)-dceU6 

5 



ych4rec6(Tavg,Lrec) =273301g . 1( p3 
ych4(Tavg,Lhc) 



ych4rec6(Tin,Lrec) =0 ych4(Tavg,Lhc) 1.81356' 10 



Bypass =0 



Lhc =3.5 'in 



Diao = 9 'in 



yc7h!4rec6(Tin,Lrec) =3.42819*10 



- 4-z Nqmavgrec7R6(z) 
Rerec6(T) Sc7T(T)-dcell6 



0.2-yc7hl4(Tavg,Lhc) 



0.2-yc7M4(Tavg,IJie) =5.99172" 10 



ycMrec4(T,z) :=ych4(Tavg,Lhc)exp 



ych4rec4(Tin,Lrec) =8.84152*10 



- 4zNmnavgrso4(z) 
Rerec4<T)«ScT(T)-dceM 



yc^rec^Tavg^Lrec) =6 945Q16 . 
ych4(Tavg,Lhc) 



10 



yc7h!4rec4(T,z) = 0.2-yc7hl4(Tavg,Lbc>exp 



- 4z*Nmnavgrec7R4(z) 



yc7M4n»4(Tin,Lrec) = 0 00773 
0.2yc7hl4{Tavg,Lhc) 



lRcrec4(T)* Sc7T( T)dcell4 ) 

^ yc7hl4rec4(TiD,Lrec)= 4.63263* 10™ 7 0.2yc7hl4(Tavg,Lhc) =5.99172' 10~ 5 0.4-7yc7hl4(Tavg,Lhc) ~8.3884'10~ 4 



- 4zNumavgreccoK4{z) 



1 - 4*z- 

ycon»4(T,z) =0.4-7 yc7hl4(Tavg,lJi C ) exp 

\Rerec4(T)-SccoT(T)-dcell4 

ycor^T^) :=04^1^Up> J - 4 z NnmavgreccoR6(2) 
j vwv j -ry 6 , F \Rerec6(T)SccoT(T).dceU6 



ycarec4(Tavg,Lrec) =3.06761M0 11 



ycorec6(Tavg,Lrec) =2.9626" 10 



15 



zx.^ = (idsp -I- 0.05)cm 



iQ yw~*(T«vg.aij (|l( J 




ruf* 



Sc7T(Tavg) =2.31842 

SccoT(Tavg) =0.73627 
ScT(Tavg) =0.67113 



ycorec4(Tavg,Lrec) = 3 656% . 1Q -| 
(0.4*7*yc71il4(Tavg,LlK>)) 



ycorec6(Tavg,Lrec) =3 53I7g . 
(0.4-7 yc7hl4(Tavg t IJic)) 



10 



"12 



This analysis assumes 80% homogeneous conversion of fuel exiting 
a 3.5in. 600 cpt monolith without bypass flow and assuming 40% 
production of CO in the post combustion monolith region. 

CO conversion in a 400 cpi recuperator monolith is >99.99%, 
while cyclohexane conversion is 99.2%. 

The final CO and C7H14 levels are < 1 ppb and .5 ppm while post 
turbine concentrations are taken as 840 and 60 ppm respectively. 
CH4 conversion also should be very high and use of a 600 cpi or 
longer monolith would give higher HC conversions (x20 for 600 cpi). 

Post turbine pre recuperation monolith pressure drops are modest 0.56% 
for 400 cpi and 0.90% for 600 cpi monoliths. 



********* *********** ****«*************************************1lr1^^ 

v - An alternative method of calculating performance uses local Nusselt numbers for heat and mass transfer to determine the 
wall temperature and combustion rates at the wall First we define the local Nu function then check conversion assuming 
fast kinetics and compare with the previous averaged (integrated) Nu formula. 

Local Nu(x) taken from Kays & Crawford for hydro & thermal entrance laminar with circular channel: 

Numoao(z) :=(I+exp(- 1.09689 - 0.381402*ln(z) - 8.287774-z 05 )) -3.6<Nwnoao(0.I I52I63) =3.82723 

zzchart : ,^hc6(Tavg)^Pr(.5 > Tavg).doell6 

2 zzchai6 = 8,90257 'cm =0.99859 

zzchar6 Y 
Rehc6(Tavg) =252.34858 * 

itest : = 0 . . 50 dce06 = 0.0967 *cm 

ztest + dwalft 

ztest = itest-0.05cm + 0.05- cm zz. : = — 

ltcat zzchar6 




Rehc6{Tavg) =252.34858 dceI16 =0.0967 -cm zz6ch4 := ^hc6(Tavg)^T(Tavg)-dceU6 

Pr(.5,Tavg) =0.72969 

Rehc6{Tavg) Pr(.5,Tavg) dcel!6 - 17.80513 cm 

^^,4 - Reho6(Tavg) Sc7T(Tavg) doeH6 

XmixT( T) -Numonof 



dcel!6 

Dch4T(T)Numono(»-^--i 
\zz6ch4/ 



^ kgpb(T,z) := \zz6ch4/ kgpb(Tm,Lfac) = 15.5934*sec 1 -cm 

dcel!6 



Dc7hl4T{T)-Numono 



kgpb7(T,2) 



\zz6c7hl4 



dcd!6 



kgpb7(Tin,Lhc)-4.50626*sec 1 'cm 



Prater Temperature for CH4 at the 
inlet with hilly developed flow 



kgpb(Ti n ,Lhc)-^- ycll40 - AHcomb(Tin > =309.73474'K 
RTin h(Tm,Lhc) 



Prater Temperature for C7H14 at 
the inlet with fully developed flow 



kgpb7(Tm,IA 0 )--^ y0o7hl4 " AHc ° mb7(Tin) =90.38489* 
* y RTin KTm.lOcm) 



Ikt tk kk t **********************************************^^ 

Determine maximum conversion of fuels with very fast catalytic combustion rates (zero fuel concentration at the wall) 

Tin+ Tout 



iz = 0.. 100 ych46 Q -ych40 



yc7h!46 0 =y0c7hl4 



Tavg 



, . j . Lhc 7 

z6, - dwali6 + xl 

42 100 7 



ych46. .Hl-kgpb(Tavg 5 z6.V— 1 ~ ^~)ych46. 

y iz+i \ & \ 55 a) ioo 7 aceitf U{Tavg)/ 



Lhc— =8.89 *cm 

7 



3.5 & 4-in. monolith length 



^-=1.2447.10- 

ych40 

/ 7^ 
ych4 Tavg, Lhc— 1 

* -^-=0.00153 



ych40 



:.(, - lgpb , (T ^,^.^.l^.^). J .™ 1 «„ 



yc7hi46 

!2i = 0.06145 

y0c7hl4 



By-pass length yC7H14out/yin yCH4outfyin 

0 3.5'in 0.06145 1.2447- 10' 4 

0 4-in 0.04256 3.288- 10" 5 

0.2 3.5-in 0.0502 5.96324- 10" 

0.2 4-in 0.03374 1.3974- 10" 5 

0.4 3.5-in 0.02054 2.17587-10"' 

0.4 3.5-in 0.01207 2.92928*10 



f7 



Tavg=952.8725-K 
Lhc =3.5 -in 
Diao = 9'in 
Bypass -0 

Nu function with circular channel hydothermal entry length 
shows 93.9% to 98.8% conversion of C7H14 



************************ ******* **************** **i 



***************************************iMr*^^ 

Pressure Drop Calculations for a Packed Bed 
given: e - void fraction; 

dp - effective particle diameter 
U - mean superfidal velocity 
v^Snematic viscocity 

p-density A x . '^^ i <>*- 

»*»*»***»**»»»****» * **»***♦*★★*****»»»*♦*»»»*»*****♦******♦*★* 

At T = 900K and 1 0 atm pressure, 



fbed(x) -—+1.75 
x 



pmix(0.0,300K) =3.47058-kg-m 
pmix(.5,Tavg) = 1. 09267 'kg-m" 3 

j =0..3 



umix(.5,Tavg)=0.00004'm 2 -sec 1 



. _ 125 . 
dp. m 

i =1.. 10 



Utx. -i-m-sec 



Lbed=0.75'm 
i)0 =umix(.5,Tavg) 



p0 =pmix(.5,T*vg) 



0.3175 
0.15875 
0.10583 
0-07938 



'cm 



Rebed . = Utx. — 



1J 1 o0(1~e) 



AP. . -Lbed fbedfRebed V 



1-e 



pO*(Utx.) 2 



*m*sec 



-1 



Rebed = 



[9 
10 








\o 0 


0 0 


f 0.00067 0.00179 0.00335 0.00536! 


10.00226 0.0054 


0.00943 0.01434! 


i 0.00474 0.01081 


0.01822 0.02695^ 


1 0.00814 0.01804 


0.02972 0.04318 \ 


[0.01243 0.02708 


0.04395 0.063041 


1 0.01764 0.03794 


0.0609 0.08651: 


1 0.02375 0.0506 


0.08056 0.11362; 


1 0,03077 0.06508 


0.10294 0.14435- 


1 0.03869 0.08137 0.12804 0.1787 ; 



0.04752 0.09947 0.15586 0.21668.; 



0.06 



0.05 



P 



0.04 



0.03 



1,2 



0.02 



0.01 
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5 6 
Utx- 



10 



dp 2 = 0.04167 in l y 



d Pl =0.0625 'in ^ 



dp = 0.125 -in 



Lbed=0.75'in 



5| AP(e,T,z) = zfbed 



Ubed(T)-dp. 



L oO-(l-e) 



1 - s pQ-(Ubed(T)r 

S 3 * 



dp = 



03175 
0.15875 
0.10583 
0.07938 



AP(.45 > Tavg,Lbed) 
P 



'cm 



0.00173 



0.00407 



0.00703 



0.0106 



e =0.45 



AP(,5,Tavg,Lbed) 
P 



0.00113 



0.00263 



0.00449 



0.00673 



Ubed(Tavg) = 2.5*m*sec 1 



For larger particles use axial beds with greater length 

k =0 dp 0 = 0.125 *in 



APaxial(E,T,z) -z-fbed 



U(T).dp k 



OO(l-E) 



1 - s P0 (U(T)) 2 



dp fc 



U(Tavg)-8.22528-m*sec 1 
APaxial(e,Tavg, I in) =0.0657 *atm 



P=3*atm 



APaxial(s, Tavg, 1 in) _ 



0.0219 1/8th inch pellets give about 2% pressure drop per inch 



********M**M*********t *************** ft**4MHt*^^ 



A********************************* 



Determine Nu values for heat and mass transfer between particles and the gas phase in the bed 



Rcpart(T) := ™«' t >* 

umixT(T) (l - e) 



Repart(Tavg) = 



398.16156 
199.08078 
132.72052 
99.54039 



0.3175 
0.15875 
0.10583 
0.07938 



Pr(.5,Tavg) =0.72969 



k=0 



e=0.45 



Sc(.5,Tavg) =0.67112 



NuhpartR(T) = 1.34- 



Pr(.5,Ty 



0,4 



0.4 (Repart(T) t )° 5 + 0.2 (Repart(T) k ) 3 



NubpartR(Tavg) =22.21552 



fcavgpart(T) 



, XmixT(T)Nalq)artR(T) 



havgpart(Tavg) « 445.592 Waf^K 1 *W 



Effective boundary 
layer thickness: 



dce!16 



NuhpartR(Tavg) 



= 142.9181 Mm 



NuhavgR6(Lhc) 



= 323.77508 -Mm 



dp k =0.3175 cm 



dceU6 -0.0967 -cm 



Nunq>arfR(T) := 1.34^ ScT( T) 0 ' 4 [ 0.4 (ReyarttDjJ 2 +■ 0.2- (Repart(T) k ) 3 



NumpartR(Tavg) =21.48438 



kgpartd) .^(^Numpart^T) 
*<*P t 



kgpart(Tm)=72.76151*sec 1 -cm 



Nupart7R(T) := 1.34Sc7T(Ty 



A4 



0.4 (RepartCT)^ 05 +^0.2 (RepMtd)^ 3 



Nnpart7R(Tavg) =35.27582 



kgfb7(T) := 



Dc7hl4T(T)'Nppart7R(T) 



kgfb7(Tin)=34.44665"sec 1 'cm 



ych46 iz+r = 



l-kgpart(Tavg) 



Lbed6(l-e) 1 
50 dp Ubed(Tavg) 



ych46.. 



ych46 5J 
ych46 



= 0.02369771 



ych4p(T\z) =ych40'exp 



6(1- c ) NmnpartR(T) z 
ERepart(T) k ScT(T) dp k 



ych4 P (Tavg,Lbed) =0Q2908261 
ych40 



yc7hl4(T,z) r yOc7hl4exp 



-fr(l -«) Nnpart7R(T) ^ 

*Rfipart(T) k *Sc7T(T) <*P fc 



yc7h!4(Tavg ) Lbed) =018611 
yOc7h!4 




dpart(cm) youtfyfnC7H14 youtfyinCH4 AP/P 
0.3175 0.18611 0.02908261 0.00113 

0.15875 0.0116 0.0000846 0.00263 



dp t =0,3175 -cm 

Ubed(Tavg)=2.5*nrsec 1 
£=0.45 

Wbed = 17.78 -cm 
Lbed = 1.905 'cm 




M Determine hydrocarbon conversions in the packed beds vs pressure to estimate effective conversion 
*B at equivalent pressure drop 



yc7h!4dp2(T,z) "y0c7h!4 exp 



6(1-8)- 



1.34*Sc7T(T) j 



0.4 



0.4(Repart(T) 2 ) 2 + 0.2(Repart(T) 2 ) 3 j 2 
*Repait(T) 7 -Sc7T(T) <*P 2 



ych4dp2(T,z) = ych40exp 



-6.(l-e)- 



1.34ScT(Ty 



0.4 



0.4- (Repart(T) 2 ) 2 +- 0.2- (Repart(T) 2 ) 3 J z 
e-Repart(T) y ScT(T) <*P 2 



Ubed(Tavg) ^2.5'nrsec 



IE 



1 

o.i 

0.01 
0.001 



>^7hl4dp2(T»vg,z id J _ 4 
y0c7hl4 



ycMO 



1*10 6 



1*10 



no 



no 




y0c7hl4 

ych40 



0.005 0.01 

^ ^ *&L>/- 



p 



0.015 




dp -0.10583 -cm 



dp 2 =0.04167 in 



z 9 = 3.74016 -in 



Lbed=0.75'in 



dp Q =0.3175 'cm 
dp Q =0.125 *in 



z 9 - 3.74016 'in 



Lbed=0.75*in 



Ubed(Tavg) =2.5 'in" sec 



For oblique flow through cone-shaped beds, large particles give a very low AP (even with thick beds) and high conversion 



zxx =2.71-cm 
APdp2(Tavg,zxx) = 



0.01000 



yc7hl4dp2(Tavg,zxx) ^ 1>32434 . 10 "5 ych4dp2(Tavg,zxx) m 5 45?93 . 10 ~n 
y0c7hl4 ych40 



For AP/P = 1%, C7H14 conversion = 99.9987% with fine particles, dp - 1 .06 mm 
zy =2.0- cm 



Repaxial(T) 



U(T)dp 0 



umixT(T)(l-E) 



Numpaxial(T) - 1.34ScT(T)' 



0.4 



0.4-(Rq>axiaKT)) 2 + 0.2 (RepaxiaKT)) 



Nump7axial(T) := l.34Sc7T(Ty 



r * % 

°' 4 10.4*(Repaxial(T)) 2 + 0.2(Repaxiai<T)) 3 



ych4paxial(T,z) =ych40 expf-6-(I - e) NroyaniKT) je_ 



yc7hl4paxial(T,z) -y0c7hl4exp -6 (1 - e) 



Nmnp7axifll(T) z 
€RfipaxiAl(T) Sc7T(T) dp Q 



o.i - 



yc7hl4paiial(Tavg,z iaap ) 001 
y0c7hl4 

ycfa40 0.001 



no 



no 



r4 




0.02 



0.04 0.06 
AP««d(s,T«vg,z^) 
P 



0.08 



Pressure drop increases 
approximately with 
thesquare of flow velocity 
while conversion 
decreases exponentially 



U(Tavg) = 8.22528* m* sec 1 zy = 2.32 cm 

AP»xial(£,Tavg > zy) r=Q02000 ych4paxial(T«vg,zy) ^ 0Q6S88284 yg 7hl4paxiaKT«vg,zy) ==0 28038 

P y ch40 y0c7hl4 

For flow perpendicular to the fixed bed, even with large catalyst particles, the equivalent pressure drop for 
/ high C7H14 conversion is very high. Even with AP/P = 2%, the C7H14 conversion is only 72%. For high conversion with 
practical fixed beds, the superficial bed flow rate must lower than about 3 m/sec to give high conversion with low AP 



* * ************** *********** ************************ ******************************************************** ******************************** #1 



Define Catalyst Functions: 



: 3_/_l l\ 

♦ ltanh(4) *] 



TOL =0.0001 



R=82.07601-K 1 'atmcm 3 atm*cm3/mol*K 
Rem = 83 1 43000* gmcm 2 sec" 2 -K~ 1 



Physical Properties for Pd/AI203 Washcoated cCatalyst 



Scat = 85 m 2 gin 1 



pcatpait : =2.5gmcm~ 



Vfipore -0.15 cm 3 - gin 1 



Loat =0.0020cm 

Washcoat thickness - 20 microns 



dpcat -0.0002 cm 

Size of Y-AI203 particle within washcoat - 2 microns 



dppore -6 



Vppore 
Scat 



efipore =0.375 
xjipore -3.0 



dMpore = 1,05882' 10 ° 'cm 



******************************** 



DeffCH4(x,T) := 



g^jpore 



1 



Dch4(x,T) 



DeflC7H14(x,T) := 



djipore* 

£^pore 
ilipore 



-•RT 



1 



Dc7hl4(x T T) 



djipore- 



-RT 



dupore* 



^RTin 



= 0.i0953*sec 1 'cm 2 



Dch4(.5,Tin)=0.41195-sec 1 *cm 2 



DeflPCH4(.5,Tin)=0 i 00419"sec 1 'cm 2 



dypore* 



-RTin 



= 0.04428'sec 1 W 



Dc7hl4(.5,Tin)=0.11846-sec *cnr 



DefiC7H14(0.5,Tin)=0.00l64-sec 1 -cm 2 



DmacCH4(x,T) :=— -Dcli4(x,T) 
2.5 

DmacC7H14(x,T) =^~*Dc7hl4(x,T) 
kOT := O.Ol gm" ! *sec~ 1 Test value for kOT 



DmacCH4( .5 , 800K) = 0.06442*sec~ 1 -cm 2 
DmaoC7H14(.5,800K)=0.0I851'sec n *cm 2 
Macropore porosity - 0.4 and tortuousity - 2.5 



Microscopic Thiele Modulus for CH4 



♦KkOT,x,T) := 



_ ^ , R-T\ kOT 
pcatpart 

ipref/ DeffCH4(x,T) 



dpcat 
6 



♦Kk0T,.5,Tin) -0.021 



Macroscopic Thiele Modulus for CH4 

<I>mac(k0T,x,T) := J pca^art f^j kOT 



Usph(»KkQT,x,T)) 
DmacCH4(x,T) 



<taac(k0T,ycIi40,Tm) =0.31774 



Microscopic Thiele Modulus for C7H14 
4n7<kOT,x,T) = pcatpart- 



kOT 



DeffC7H14(x,T)J 6 



dpcat 



♦ji7(k0T, . 5, Tia) =0.03357 



Macroscopic Thiele Modulus for C7H14 

<W7(k0T,x,T) 



pcatpart(|^)k0T 



ijsph(fr7(k0T,x,T)) aLwt Omac7(kOT,ych40,Tm) -0.59265 



DmacC7H14(x,T) 



RatekCH4(k0T,x,T)=-n^ 



pcatpart-)— i-kOT 



RatekC7H14(k0T,x,T) = -7tsph(ft7(IOT,xJ)) ^ 



RatekCH4(k0T,0 7 Tro) =-1610.07113 sec 



^•[^l-kOT 



RatekC7H14<kOT,0,Tm) =-1492.981 86 sec 



Defines transport restricted rate constant function given catalyst effective first order rate constant - 
kOT, methane or cyclohexane conversion - x, and surface temperature - T 



Ts-Tin kOT=0.01*sec 1 -gin 1 Tavg = 952.8725" K 

kgpb(Tin,0.1'cm) =26.77306- sec" 1 -cm RatekCH4(k0T,0,Ts) ="1610.071 13 -sec' 1 Lcat =0.002 -cm 



Th+ _ kgpbdavg^O.l.cm) ^ = ^ 

kgpb(Tavg,.l>cm) \ ^ h(Tavg,.lcm) R-Ts 
T,0,Ts)-Lcat/ 



RatekCH4(k0T, 



xs=ych40 xs =0.01183 Tin =8 10.928 'K 

Rk:=6.33-m" 2 -sco" 1 Ratech4 = 0.036 m 2 - sec 1 kOT = lOicg^-sec" 1 

^ Set up subroutine/function to solve for the surface temperature given a rate constant 
03 TOL = 0.0001 ych40 =0.01183 x:=0.01 
fi Ts ^1673-K Ts7-=1473-K 
Given 

p 

™ = y( ych40 , x) A 

W TsPTin + - m«*MT S ) 

1 1 hTin,^ 



-RatekCH4(kOT,x,Ts) Lcat . ./_ Lhc\ \ 100 

kgpb l Tin 'ioo) 

Tshi(k0T) = Find(Ts) Tshi( 100k0T) = 1010.1063 1'K 

Given 

y7(y0c7hl4,x) i: *-^- 

rp _ T . 6 R Ts AHcomb7(Ts) 

Ts7— Tin +- - } T< v 

1 1 hTin,^ 



-Ratd£C7H14(kOT,x,Ts)-Lcat . ._/_. Lhc\ \ 100/ 
Tshi7(k0T) - Find(Ts7) Tshi7( 100k0T) = 873.32275*K 

Rk=-RatekCH4(k0T,xs,Ts>Lcat Rk7 =-RatekC7H14(k0T T xs,Ts)Lcat 

Rk=6.72632*sec" 1 "cm Rk7 =6.4 1773 "sec" 1 'cm kOT = 10*kg n -sec' 1 

P 

xs 

^ Ratech4 :=— kgpb(Tin,Lhc) = 15. 5934 -sec" 1 'cm Ratech4 =0.012 15 'nT 2 ' sec 1 

— + 

Rk kgpKTin,Lhc) kgpb7(Tin,Lhc) =4.50626*sec 1 'cm 



I Ts:=Tin + Ratecl,4 - AHcomb(Ts) ych40 _ . _ 001 „ n 

h(Tm,Lhc) Ratecn4 -0.012 15 in sec 

Rk = 672632- sec 1 -cm 

Ratech4 

xs 

Rk 

Ts = 811.4626*K xs = 0. 18057 'm 3 Tin =8I0.928'K 

kOT = lO-kg^'sec" 1 Tshi( 100-kOT) = 1010.10631'K Tshi(0.01-kg l sec l ) -810.96346-K 

Tshi7(10Ok0T) = 873.32275'K Tsai7(0.01-kg" 1 sec" 1 ) = 810.96393'K 

xn=0,L.30 

kOTcalc^ := 10 5 -kg" 1 sec" 1 Tsnrf^ = Ts^kOTcatoJ Tsurf7 xn = T^7(k0TcakJ 



Tswf 



Tsurf7 




This plot shows the surface 
temperature within a honeycomb 
channel as a function of the specific 
catalytic rate constant for methane 
and methytcyclohexane combustion 
at a location within 1-mm of the 
leading edge 

Takes into account the internal pore 
diffusion (micro- and macropores) and 
boundary layer heat and mass 
transport. Upper surface temperature 
limit approaches the adiabatic 
temperature for methane and a much 
lower (Prater) temperature for C7H14 

kO unit - secM per kg washcoat 
Tunit-K 



We can also determine the actual effective first order rate constants for methane, propane, and other fuels and 
estimate the temperatures required for catalytic light-off and extinction, with roughly estimated Nu numbers. 

Other factors, such as the upstream surface temperature and gas phase radial temperature profile and thermal 
conduction of the catalyst wall, can also effect the light-off and and the critical extinction points. These calculations 
are beyond the scope of the current analysis. 



************** ************hh**************i 



**************** 



PdO Catalyst Rate Constants 



EipdO = 12800-K Erpdl =4900 K Erpd3 = 15000K 



RegionO-lowTPdOto625K Regionl - Intermediate T Region2-Flat909K Regions- 
ignores inhibition by H20 and C02 625Kto909K to 1048K decomposition Pd metal 



ApdO^lt^WgmW 1 / Eipd o Erpdl \ . « . ... /-Erpdl 

ApdS-Apdlexpi-^^' 



Interpolation formula for the four PdO+Pd rate constants 



kpd(T) 



Apd3^(-^) + 



0.5 



Apd2- + Apdl + ApdO- W^>) 



T / -V T "j/J 



nx =0,1.. 41 T =500K+nx-25K 



CATALYTIC RATE CONSTANT NEEDED FOR HEAT BALANCE IN A 
MONOLITH WITH FULLY DEVELOPED FLOW 



Compared with extrapolated but measured 
rate constants (mol/atmCH4/s/gcat) 



Lcat =0.00002-iii 



kOTcalc. 




0.001 



1.2 1.4 
1000-K 100OK 
T _ 'Tsurf 



Compared with extrapolated and estimated 

rate constants for C7H14 (lOx mol/atmCH4/s/gcat) 



kOTcalc 




0.6 



0,8 



1.2 1.4 
1000-K 1000-K 

T ' Tsurf7 
nx xn 



1.6 



1 j Define functions to calculate surface temperature profiles given gas temperature and composition 

x = 0.01 



TOL =0.1 Ts=81L4626*K 
M Tsx := 1500-K Tsx7 := 1400-K 

Given 



1 



.+ 1 = 7.01896*m l *« 



~RatekCH4(kpd(Tsx),x,Tsx)*Lcat kgpb(Tm,Lhc) 



y(ych40,x) 



Tsx = Tin+ 



°RTsx 



-AHcomb(Tsx) 
h(Tin,z) 



- RatekCH4(kpd( Tsx) , x, Tsx) Lcat kgpb( Tin, z) 



Tsrf(Tin,x,z) -Find(Tsx) 
Tsch4x -Tsrf(Tin,x,0.1cm) 

Tsch4x = 823.23891*K Tin = 810.928'K 

kpd(Tsch4x) =3.5131I-kg~ I -sec~ 1 
kgpb(Tin, 0.1 cm) =26.77306* sec 1 "cm 



Z inlet = 0.1 cm 
x = 0.0 1 Lcat = 0.002 -cm y(ych40 , x) Q = 0.01 17 1 



-RatekCH4(kpd(Tsch4x),x,Tsch4x)Lcat = L17316*sec 1 "cm 



Given 



y7(y0c7hI4,x) 



P 



Tsx7=Tin+ 



-AHcomb7(Tsx7) 



1 1 



h(Tm,z) 



-RatekC7H14(10kpd(Tsx7),x,Tsx7).Lcat + kgpb7(Tin,z) 



Tsrf7(Tin,x,z) = Find(Tsx7) 



Tsc7hl4x : = Tsrf7(Tin,x,0.1*cm) 



Z inlet = 0.1 cm 



Tsc7hl4x = 871.51735'K 



Tm = 810.928'K 



x = 0.01 



Lcat = 0.002 -cm 



y(ych40,x) 0 =0.01171 



10kpd(Tsc7hl4x) = 44.23122*kg 



sec 



-RatekC7H14( 10*kpd(Tsch4x) t x,Tsc7hl4x)Lcat 



9.21644-sec 



kgpb7(Tin,0.1cm) = 12.3599-sec 1 -cm 



yn Assuming the rate constant for C7H14 combustion is ten times the rate constant for CH4 combustion, the surface 

yg temperature near the inlet (1-mm in) is 871 .5 K (1 1 0OF) and near the transport limit for C7H14 and 823.2K (1 050F) and 

Q surface rate controlled for CH4. If this calculation is repeated for a series of steps downstream with corrections to the 

;" : \ gas phase concentrations and temperatures, new surface temperatures and rates can be calculated throught the length 

!J? of the monolith. 




